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Abstract 


This work revolves around the Jaumann wave absorber, its solutions (either obtained 
through maximally flat solutions and through genetic algorithms), its possible im- 
plementation for designing wave-absorbing and radar-cloaking systems, and mea- 
surements of ambient RF energy in four Duke buildings to gauge the feasibility of 
harvesting energy from the RF band. 

The first part of the introduction provides a short recap of radar, electromag- 
netic wave absorbers, and their equivalence to transmission line systems. The next 
part describes a brief summary of the history of the Jaumann wave absorber, its 
development over the years, and the short comings of previous attempts of finding 
its solutions. Finally, we conclude the introduction with contributions of this work. 

The second chapter presents the mathematical formulation of Jaumann absorbers 
which is crucial in deriving their maximally flat solutions. This chapter lays the 
foundation to explore other solutions through cut-and-try and genetic algorithms. 

The third chapter expands the search for solutions through automated means. We 
first try the simple cut-and-try method before moving on to a more sophisticated 
genetic algorithms. The bulk of this chapter describes the details of genetic algorithm 
and how they are implemented into our Jaumann wave absorber problem. 

The fourth chapter provides the results and analysis obtained from the algorithms 
described in the previous chapter. However, only samples of such results are pre- 


sented in this chapter while the bulk of the results is moved to Appendix A.9 for the 


sake of brevity. 

The fifth chapter describes our experimental works in verifying a sample result 
of the Jaumann wave absorber and the measurement of ambient RF energy in four 
Duke buildings. We also include the measurement results and analysis in this chapter. 
A brief comparison for the ambient energy measurement to other studies in other 
locations and the environment is also included. 


Finally, the sixth chapter concludes this work. 
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Introduction 


After James Clerk Maxwell published his famous equations, scientists and engineers 
have predicted the existence of electromagnetic waves. Only a mere 8 years after 
Maxwell's departure, Heinrich Hertz made the first observation of microwave elec- 
tromagnetic radiation [1]. Soon after, Nikola Tesla envisioned a connected world 
through wireless technology. 

Development of wireless technology has brought us the invention of RADAR 
(RAdio Detection And Ranging) during World War 1. Radar is mainly used for 
detection of an object, specifically to determine the distance, speed, and/or size of 
the object. A radar transmitter shoots high frequency radio waves and a receiver 
detects the reflection from an object. A radar antenna can act as both transmitter 
and receiver, switching between the two roles alternately. As a transmitter, the 
antenna is usually designed to focus on one direction to achieve a high gain but 
rotated around to achieve a wide detection area. As a receiver, it merely waits for 
the reflected wave to travel back. The time delay between transmitted and reflected 
wave determines the distance between the transmitter and the detected object. A 


more sophisticated radar system also detects the Doppler shift between transmitted 
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and reflected wave to determine the speed and direction on a radial axis with respect 
to the transmitter. The detection of size, however, is tricky since a radar system 
relies on the strength of the reflected wave to figure out the "size" of the detected 
object which may not be equal to the object's actual physical size. In particular, the 
stronger the reflected wave, the bigger the "size" detected by the radar system. This 
”size” is called Radar Cross Section (RCS) which will be discussed briefly in Section 
1.1. 


1.1 RCS Fundamentals 


RCS is defined as the effective area of a detected object as if it is the only area which 
receives the transmitted radar signal and perfectly reflects the radar signal back to 
the receiver. In other words, if we substitute the detected object with a perfectly 
reflecting sphere, the cross section of the sphere is the RCS. Mathematically speaking, 


the RCS is expressed as follows: 


RCS = c = (1.1) 








where 


RCS (o) = Radar Cross Section (m?) 
P. = reflected power (W) 


P, = transmitted power (W) 


G, = gain of transmitter antenna () 
r = distance between transmitter and receiver (m) 
Aur = effective area of receiver antenna (m?) 


In most cases, especially in the development of stealth technology, we are in- 
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terested to reduce the RCS of detected objects. On the contrary, designs of the 
commercial airplane would seek to maximize the RCS of the airplane for example. 


Several methods of reducing RCS are: 


e designing the surface of the object in such a way that the radar wave is reflected 


as much as possible to any direction except the incoming wave [2], 


e guiding the incoming wave around the object so that the transmitter receives 


no reflected wave back [3, 4], 
e coating the object with Radar Absorbent Material (RAM), [5, 2] 


e transmitting a wave which destructively interferes with the incoming radar 


wave (active cancellation stealth) [6, 7] 


In this project, we primarily focus on reducing RCS by designing an electromagnetic 


wave absorbing system. 
1.1.1 One layer absorber 


Consider a wave absorber placed 2 away from a conducting sheet as shown in Fig. 1.1. 
Using transmission line equation (Eq. 1.2) and setting o = 0 for lossless condition 


in Eq. 1.3, we obtain Eq. 1.4. 


Zr + Zo tanh (Yat) 
Zi = Z 1.2 
Ë g Zo + ZL tanh (Yat) ( ) 





Ya Qa + jpa (1.3) 


Zi = Zotanh(jñ,t) = jZotan (Bat) = co (1.4) 





e 
g 
g 
— 
7 
Ce 
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g 
lossy M conducting 
material alr sheet 


FIGURE 1.1: Wave absorber placed quarter wavelength away from a conducting 
sheet. 


Similarly, for the case inside the lossy conductor (0 < x < t), which ABCD matrix 


is expressed in Eq. 1.5: 

















A B |. | cosh (yt) Zosinh yt (1.5) 
C D ee | Yosinhyt  coshyt I 
AZ, + B A Zi 
Zn = mea e = : (1.6) 
CZL +D C tanh (yt) 
If (yt) is small : 
Bh X 
tanh (yt) ^ yt 
jop 
= o+jwe 
jwp(o + jwe) t 
1 
= WA 1.7 
(c + jwe) t Man) 


If c > jwe: 


1 
Zin = — = 3770 (1.8) 
ot 


Thus, given a specific value for sigma, the thickness can be tailored to achieve 
maximum absorption. Similarly, for square resistance sheet (width (W) = length 
(L)) we obtain: 

1L L 1 


oA oWt ot ER 





where R, is sheet resistance. Therefore, we can model the wave absorber system 
with a transmission line circuit as shown in Fig. 1.2 where Ryn = Zp. 

However, doing it this way yields a narrow bandwidth of absorption, due to the 
requirement of 1/4 separation between the conducting sheet with the absorbing layer. 
In the interest to improve the bandwidth, we use Jaumann wave absorber system 


which is discussed more in section 1.2 [8, 9]. 


E sk W 
E 


Incident í L 
wave Rs 







Conducting 
sheet 

















Incident 
wave Short 


circuit 


FIGURE 1.2: Wave absorber layers (top) and transmission line equivalent (bottom). 


Finally, Fig. 1.3 shows our proposed circuit to one layer Jaumann wave absorber 
model. Incoming E-field induces voltage difference as given by V = EL. The diode 
bridge rectifies the induced current so that the resistor always have polarity as indi- 
cated in the figure. In addition, we can also tune the length of conducting connectors 
so that the amount of stray inductances cancels out the innate capacitance present 


in the diodes. 


Antenna 











FIGURE 1.3: Proposed circuit equivalent to one layer wave absorber model. 


1.2 Background on Jaumann Wave Absorber 


Electromagnetic wave absorbing materials were first developed during WWII era 
when both the axis and allies were competing to build radar technology to detect 
enemies’ submarines and prevent their own from being detected. Therefore, a wave 
absorbing layer which renders a metal surface undetectable by radar was considered 
crucial for winning the war. On the axis side, a German engineer named J. Jaumann 
designed a set of solutions of impedance values which absorbs at least 90% of radar 
signal power between 1 - 10 GHz [8]. Meanwhile on the allies side, an American 
inventor W. W. Salisbury patented his design of an electromagnetic absorbent body, 
an early idea from which wave absorber layers today are named as Salisbury screens 


[10]. 


Throughout the rest of the 20th century, development of radar absorbing ma- 
terials (RAM) continued [9, 11, 12, 13, 14, 15, 16]. However, efforts in improving 
multiple layers of wave absorbers so far has been limited either by lack of compu- 
tational power or the scope of its application. For instance, Jaumann produced his 
solutions through a trial-and-error process assisted by an American invention called 
the Smith chart. He started from the zero impedance point (short circuit) in the 
chart, rotated A/4 and added in parallel an arbitrary impedance value Zi, rotated 
A/4 and added in parallel another arbitrary impedance value Z», and repeated this 
process until he ended on the center point of the chart (no reflection). The set of 
solution Zi, Z2, ..., Zn was then tested for its reflection coefficient. If the reflection 
coefficient was unsatisfactory, then the whole process was repeated to search for a 
new solution set. In another example, Knott et al [12] presents a more systematic 
approach by using an S-matrix to solve for maximally-flat response of the Jaumann 
absorber. However, they only provided the solution for 2 layers. Meanwhile, Fante 
and McCormack [14] used the multiple wave reflections approach to calculate the 
admittance for 1, 2, and 3 layers of electric screens, but their approach became 
incredibly complicated for more than three layers. Another work by Nortier et al 
[13] provides a design table for up to 7 layers, but their paper did not explain how 
they arrived at their values. Moreover, their table entries produce 511 plots with 
significant ripples in the pass band (not maximally-flat). 

In this work, we have successfully expanded the maximally-flat solutions of Jau- 
mann wave absorbers up to 6 layers or more, an improvement over previous attempts 
recorded in the literature (9, 11, 12, 17, 13, 14, 15, 16]. Furthermore, our results not 
only produce the maximum absorption at center frequency, but the mathematical 
procedure we developed can also be expanded to solutions of any number of layers, 
given enough computational resources. More details on this is presented in chapter 


2; 


In chapter 3, we describe how we use the result from mathematical procedure 
from chapter 2 and feed it into genetic algorithm to search for more solutions. 
Moreover, we can fine tune the parameters in genetic algorithm to obtain solutions 
which yield minimum bandwidth, desired minimum reflection coefficient at center 


frequency, and/or center frequency of the bandwidth. 
1.3 Contributions 
Our contributions are summarized as follows: 


e We produced a theoretical framework to obtain maximally flat solution of Jau- 


mann wave problem. 


e Using this theoretical framework, we explore other solutions to Jaumann wave 
problem using Genetic Algorithm. We compile 240 unique solutions, (120 each 
for -10 and -20 dB reflection coefficient limit; for 1 to 6 layers with 20 solutions 


each for frequency range which corresponds to fo - 20 fo for each solution). 


e In order to test the viability of wave absorber system as energy harvester, we 
measured the RF ambient energy level on 4 Duke buildings (CIEMAS, Teer, 


French, and Physics). The report of which is under preparation. 
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Mathematical Formulation of Jaumann Absorbers 


2.1 Determining Maximally Flat Solutions Through ABCD Matrix 


2.1.1 One section 


Zs = Zo Zi 





ZT Zi 
FIGURE 2.1: One section transmission line circuit 


By using ABCD matrix, the circuit shown in Fig. 2.1 can be represented as: 


u 1 0 cosÜ ` Zo sind 
dee = H BP cos 0 | 


cos jZo sin 0 
G,cos0 + jYosinO JG, Zo sin 0 + cos 0 


Where G = Ea 


And the total impedance of the circuit is: 





AZr+B B 
Zr = — 
CZ,+D D 
jZosind 





= 2.2 
jG1Zosind + cos 0 22) 


Setting Zo = 1 and Zr = 1 (which means we normalize Zo and set Zr = Zs for 


maximum power transfer), we obtain: 


jG, sin 0 + cos? — jsin0 = 0 
J sin (Gi — 1) + cos 0 = 0 (2.3) 


À 


When the length of transmission line is 4, we can solve for Gj = 1, where G is 


the normalized conductance. 


2.1.2 Two sections 


Zs Zo Zo 





FIGURE 2.2: Two sections circuit layout 


The ABCD matrix for the circuit shown in Fig. 2.2 is: 


J 1 0 cosÜ =f Za sind 1 0 cosÜ | j Zy sin 8 
Aa P: AAR cos Ø el Mel r cos 0 


= cos 20 + jG1Zo cos Ø sin 0 Zo sin Ø (2j cos 0 — Gi Zo sin 0) 
zE (G1Zo + (Gi + 2G2) Zo cos 20 + j(2 + G1G2Zê) sin 20) cos? Ø + j(G1 +2G2)Zo cos 0 sin 0 — (1 + G1G»Z2) sin? 0 


(2.4) 
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The load Zr with Zz = 0 (due to short circuit at the end) is: 




















7 AZ B. Æ 
T TGN UD 
8 Zo sin Ø (2j cos 0 — G1 Zo sin 0) „ jsin?g 
cos? 0 + j(G, + G2)Zo cos 0 sin 0 — (1+ G,G,Z2)sin20 ` jsin?0 
8 229 cot 0 + jG, Zó (2.5) 
` Zoeot O(G | + 2G3) + j(1 + G1G3Z2 — cot? 0) ` 
Setting Zo = 1 and Zr = 1, we obtain: 
cot OO: + 2G2) + J(1 + G1G9 — cot? 0) = 2cot6 + jG, 
cot 0(G1 + 2G — 2) + J(1 + G1G3 = Gi — cot?) = 0 (2.6) 


Note, make sure that all present trigonometric functions don't turn into an inde- 
pendent constant when 0 — 7/2 (any independent trigonometric function must go 
to zero, but trigonometric function which attached to a set of other parts are fine). 

Since we have set Zr = 1, setting all the coefficients equal to zero (in addition to 
0 — 7/2) in Eq. 2.6 will give the G: and G> values which satisfies the impedance 


matching: 


| Gi +2Go—2=0 (2.7) 


1 + G1G» ki G4 = 0 
Solution is: Gj = 1.4142, G> = 0.29289. 


2.1.3 Three sections 
1 0 cosÜ | jZosin0 1 0 cosÜ | jZosin0 
“| jYosnO cos? | `` Gi 1| | jYosinO  cos0 


(2.8) 
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E 
FIGURE 2.3: Three sections circuit layout 


The result of matrix multiplication is: 


A = eos° 0 + j(2G1 + G3)Zo cos? Ø sin 0 — (3 + G1G2Zj) cos 0 sin? 0 — jG» Zo sin? 0 


B = (1/2)jZysin 0(2 — G1G2Z3 + (4 + G1G3Z$) cos 20 + 23(G + G3) Zo sin 20) 





1 
C = IZ. (26,20 + JC Zu = G1G2G3Z3) cos 0 + Zo(2G1 | 2G» | 4G3 | G1G3G3Z2) cos 30 + 
0 


D = cos*0 + j(G4 + 2G2 + 3G3)Zo cos? Ø sin 0 — (3 + 2G2G3Z2 + G1(G2 + 2G3)Z2) cos O sin? 0 — 
jZo(G1 + G3 + G1G2G3Z3) sin? 0 (2.9) 
The load Zr with short circuit load is: 


AZ,+B B 


Zp = = 
á 0715 D 





Zo sin 0(—3 cos? 0+sin 0(—2j(G1+G2)Zo cos 0+(1+G1G2Z2) sin 0)) 
j cos? 0— (G4 +2G2+3G3)Zo cos? 0 sin 6—j(3--2G53Ga ZG (G2 -2G3) ZE) cos 0 sin? Ø+Z0(G1+G3+G1G2G3Z2) sin? 0 





(2.10) 
Setting Zo = 1 and Zr = 1 and cross multiply, we obtain: 


cos? 0(—3 + G + 2G2 + 3G3 — j cot 0) + ¡(3 + 2G2(—1 + Gs) + G1(—2 + G> + 2G3)) cos 0 sin 0 
—(—1 + G4(1 + G,(—1 + G3)) + Ga) sin? 0 = 0 


(2.11) 


12 


Therefore, the system of equations that we need to solve are: 


—3 + G1 + 2G2 + 3Gs = 0 
3 + 2G3(—1 + Gs) + G4(—2 + G> + 2G3) = 0 
-1+G1(1+ G>(—1 + Gs)) + Gs = 0 





(2.12) 


And the solution is: G: = 1.644, G> = 0.5143, Gs = 0.1091. 
In the interest of brevity, list of system of equations for four, five, and six sections 
case are moved to Appendix A.1. Table 2.1 lists the summary of solutions and the 


resulting relative bandwidth (RBW). 


Summary of the solutions 























Table 2.1: Summary of the G values up to six sections. 
# of sections | Gi Ga G3 G4 Gs Gs | RBW 
1 1.000 - - - - - 0.253 
2 1.414 | 0.293 - - - - 0.673 
3 1.644 | 0.514 | 0.109 - - - 0.926 
4 1.782 | 0.670 | 0.234 | 0.044 - - 1.089 
5 1.868 | 0.778 | 0.343 | 0.114 | 0.019 - 1.203 
6 1.921 | 0.853 | 0.430 | 0.188 | 0.057 | 0.008 | 1.288 





















































Table 2.2: Impedance values for the case Zo = 50 Q up to six sections. 
# of sections | Zi (Q) | Za (Q) | Za (Q) | Za (Q) Z; (Q) Ze (Q) 

1 50.000 - - - - - 

2 35.356 | 170.713 - - - - 

3 30.413 | 97.220 | 458.295 - - - 

4 28.052 | 74.683 | 213.950 | 1128.668 - - 

5 26.768 | 64.292 | 145.900 | 440.141 | 2659.574 - 

6 26.029 | 58.644 | 116.198 | 266.382 | 883.392 | 6097.561 
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Table 2.3: Impedance values for the case Zo = 377 Q up to six sections. 
# of sections Zi (Q) Z2 (Q) Za (Q) Za (Q) Zs (Q) Z6 (Q) 
1 377.000 - - - - - 
2 266.582 | 1287.172 - - - - 
3 229.318 | 733.035 | 3455.545 - - - 
4 211.513 | 563.107 | 1613.179 | 8510.158 - - 
5 201.831 | 484.763 | 1100.088 | 3318.662 | 20053.191 - 
6 196.262 | 442.177 | 876.133 | 2008.524 | 6660.777 | 45975.610 
2.2 Bandwidth Enlargement 
2.2.1 Reflection coefficient for one section case 
Earlier we found that the Zr for one section is: 
AZ, +B B 
Zr = = 
CZ,+D D 
— eosin’ (2.13) 
jG1Zosind + cos 0 
Then, the reflection coefficient is: 
Zr— Z 
B. uoc 
Zr + Zg 
|. —(Zscos0 + jZo(—1 + Gi Zs) sind) (2.14) 
H Zs cos 0 + jZo(1 + G1Zs) sind ` 





Assuming the case where Zç = Z = 502. Substituting the de-normalized 


solution for one section (Gj = Sch we have: 


— cos 0 
cos Ø + 27 sin 0 





(2.15) 
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2.2.2 Reflection coefficient for two sections case 


AZ,+B B 


Z === — 
S CZ,-D D 








—2 Zo sin 0(—2j cos0 + G4Zosin 0) 
= | | (2.16) 
(2 + G1G2Z2) cos 20 + Zo(-G1G229 + (Gi + 2G>) sin 20) 





Then, the reflection coefficient is: 


ZT — Zs 


e ee 
E Zr + Zs 


(-2Z5+G1Z3(1-G2Z3)) cos 20+Z0(G1Z0(-1+G2Z5)-j(-2+G1Z5+2G2Zs5) sin 20) 
((2Z5+G1Z3(1+G2Z3)) cos 20+Z0((-G1)Zo(1+G2Z5)+j(2+G1Z5+2G2Z5) sin 20)) 








(2.17) 


Assuming the case where Zç = Z = 500. Substituting the de-normalized 





solution for two sections (e: = SC Ga = 223), we obtain: 


0.447 — 2.605 cos? Ø — 0.6327 cos Ø sin 0 
— 1.828 + 3.829 cos 20 + 8.0003 cos 0 sin 0 





I» = 0.0794 (2.18) 


2.2.8 Reflection coefficient for three sections case 


AZr+B B 
Zr — — 
CZ+D D 








Zo sin 0(—3 cos? 0+sin 0(—2j(Gi +G2)Zo cos 0-- (1--G1 G5 Zê) sin 0)) 
j cos? 0— (G1 +2G2+3G3)Zo cos? 0 sin Ø-5(3+2G2G3Z3+G1(G2+2G3)Z2) cos 0 sin? Ø+Z0(G1+G3+G1G2G3Z2) sin? 0 





(2.19) 


Then, the reflection coefficient is: 


, 


Bb = ros 
5 Zu. 


LI (—Zs cos? 0—jZo(—3--(G1--2G2--3G3)Zs) cos? 0 sin 0+(3Z5+Z3(-2(G1+G2)+(G1G2+2(G1+G2)G3) Zg)) cos 0 sin 0? -jZo(-1--GaZs-Gi1(Zs--Go Zà ( -1--GaZs))) sin 03) 
(Zs cos? 0+jZ0(3+(G1+2G2+3G3)Z5) cos? Ø sin 0—(3Zs--Zg (2(Gi+G2)+(GiG2+2(Gi+G2)G3)Zs)) cos 0 sin 0? — jZo(1--Ga Zs -G1(Zs--Ga Zà (1+G3Zs))) sin 03) 





(2.20) 
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Assuming the case where Zs = Zo = 500. Substituting the de-normalized 





solution for three sections (e: = mE Go = SE Gs = m), we obtain: 


L. = —0.5j cos 0(1 — cos 0) (2.21) 
SE —1.4087 cos 0 + 2.4087 cos 30 + 1.268 sin 0 — 2.422 sin 30 ` 





2.2.4 Reflection coefficient for four, five and six sections case 


The expressions of the reflection coefficients after de-normalization and substitution 


of the G values for each case: 


m= — cost 0 (2.29) 
^ 0.085 — 1.795 cos 20 + 2.880 cos 40 — 3.5327 cos Ø sin 0 + 2.8837 sin 40 ` ` 





0.6257 cos 0+0.3127 cos 304-0.063; cos 50 


I's 0.1127 cos 0+2.2437 cos 30—3.3567 cos 50—0.072 sin 0—2.237 sin 30+3.356 sin 50 


(2.23) 


T — 0.312+0.469 cos 20+0.187 cos 40+0.031 cos 60 
6 ^" 0.024+0.100 cos 20+2.714 cos 40—3.837 cos 60+-0.177j cos 0 sin 0+2.712j sin 40—3.838j sin 60 


(2.24) 


2.2.5 Plot of reflection coefficients 


Fig. 2.4 shows the plots of Eq. 2.15, 2.18, 2.21, 2.22, 2.23, and 2.24. The plots show 
that = 0 when 0 = 5 nr, n = 0,1,2, .... 


where: 


Or (2.25) 
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Reflection coefficients vs frequency 





IF(9)] (dB) 








— One segment | 
=== Two segments 
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f (MHZ) 


FIGURE 2.4: Bandwidth plot, y-axis in dB and x-axis in Hz 























f = frequency 


fo = center frequency 
2.2.0 Proof of maximally flat response 


Obtaining a solution which yields maximally flat response requires that the reflection 
coefficient and its derivatives be set to zero at the center frequency. Despite of much 
more algebraically tedious, we will see that this approach also leads us to the same 


solution as ABCD matrix approach. 
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One section 


Starting from Eq. 2.14: 


L, 

Setting Zs and Zo = 1: 
L, 

Setting 0 — E 


I 


—(Zs cos + jZo(—1 + G1Zs) sin 0) 
Zg cos Å + jZo(1 + G1Zs) sin 0 








_ cos0 + j(—1 + Gi)sin0 
cos 0 + (1 + Gy) sind 





tæ 
uS (2.26) 





Solving for I' = 0 gives Gj = 1. 


Two sections 


From Eq. 2.17 and setting Ze and Zo = 1: 














(2.27) 








r = (G4(—1 | G3) | (-2 G4 T G1G3) cos 20 — j(—2 + G4 + 2G2) sin 20) 
m (-G4(1 | G3) | (2 G4 | G1G3) cos 20 + SCH + Gi + 2G3) sin 20) 
dT, — (43(4 — G? + G2 cos 20 + 2G, sin 20)) 
dü ((—G4)(1 | G3) | (2 | G4 | G1G3) cos 20 + 7 + G4 + 2G3) sin 20)? 





Setting 0 > 5: 








(2.28) 
CEC 
ue (2.29) 
2j(—2 + G? 





(1+ Gi + G41G3)? 
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Solving for I's = 0 and do = 0 yields: 
G, = 1.414 
Gə = 0.293 
Three sections 


From Eq. 2.20: 


(—Zs cos? 0—jZo(—3+(G1+2G2+3G3)Z5) cos? O sin 0+(3Z5+Z3 (-2(G1+G2) +(G1G2+2(G1+G2)G3) Z5)) cos 0 sin 6? +jZo(-1+G3Z5+G1(Z5+G2Z3(-1+G3Z3))) sin? 0) 
(Zs cos? 0+jZo(3+(G1+2G2+3G3)Zs) cos? 0 sin 9-(3Zs +Z (2(Gi+G2)+(GiG242(Gi+G2)G3)Zs)) cos O sin? 0—j Zo 1+G3Z5+G1(Z5+G2Z3 (1+G3Zs))) sin? 0) 





Ts 





(2.31) 


Setting Zs and Zo = 1: 


— cos? Ø-j(-3+G1+2G2+3G3) cos? 0 sin Ø4(3+G1G2—2(G1+G2)+2(G1+G2)G3) cos 0 sin? 0+j(—1+G1(1+G2(—1+G3))+Gs3) sin? 0 








Ls = cos? 0+7(3+G1+2G2+3G3) cos? O sin 0— (3+2G2(1+G3)+G1(2+G2+2G3)) cos sin? 6—j(1+G3+Gi (1+G2+Go2G3)) sin? 0 
Fa (j(24—8G2—4G2--3G1 G3 —4G 1 (2G2+G 1 (—2--G2)) cos 20+G2(8G1+4G2+G7G2) cos 40+16jG) sin 20—4jG2G» sin 20—8jG1GÀ sin 20+8jG2 sin 40+2jG2G2 sin 40+-4jG1G3 sin 40)) 
dø (4(cos? 0+j(3+G1+2G2+3G3) cos? Ø sin 0-(3+2G2 (1+G3)+G1 (2+G2+2G3)) cos Ø sin? Ø-j(14G3+G1(14G2+G2G3)) sin? 0)?) 





(2.33) 


(The expression for ¿ls is moved to Appendix A.2 since it’s too long.) 





Setting 0 > 5: 


P. TL O+ Ga(-14 Gy) + Ga (2.34) 
? 1461616) T Gs + GiG5Gs | 





(2.35) 





ds — — (23426165 + Gi(-2 + G2))) 
d — (14+- Gg + Gi1(14- Go + G5G3)? 














drs — (4(9+4G2(1+G3)+G1(8+7G2+8G3)+G7 (—4+G2(1+3G2+G3))+GF (—4(1+G3)+G2(—3+G2(1+G2+G3))))) (2 36) 
d02 (1+G3+G1(1+G2+G2G3))3 i 
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Solving for T3 = 0, ae = 0 and ¿Te = ( gives: 
Gi = 1.644 
G» = 0.514 
G3 = 0.109 


Since the results from maximally flat approach and ABCD matrix approach are 
exactly the same, this proves that the results from section 2.1 gives maximally flat 


response. 
2.3 Jaumann Absorbers with Varying Separation 


In section 2.1 we have already discussed on deriving the mathematical expression of 
Jaumann absorber problem which we will use in Genetic Algorithm (GA) in chapter 
3. However, we can step further into the search of optimal solution by generalizing the 
length of transmission line separation between the absorbers (which was previously 
assumed to be 4/4). 

On the first thought, varying the length of transmission line simply means as- 
signing different 0, to each section. However, separating the 0 into 0,, while simul- 
taneously varying 6, to plot 51; would add complication in the algorithm. To avoid 


this problem, we introduce a new variable X,, defined as: 


0, = X,0 (2.37) 
L 
Xo dmm 2.38 


Where L is the length of the separation between absorbers. We define X„ and L 
in this way for the sake of consistency with our preceding calculations. We maintain 
Xn = 1 for the quarter wavelength sections as the default and let the rest scale with 
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respect to quarter wavelength. This idea allows us to simply add X, in front of every 
0,, inside the trigonometric functions to control the varying separation length. 

In the interest of maintaining the brevity of this chapter, we place the full ex- 
pression of the reflection coefficient I" for varying separation case in Appendix A.4 


along with the MATLAB code. 
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3 


Genetic Algorithms and Their Application to 
Jaumann Wave Absorber Problems 


3.1 Motivation 


Although maximally flat solutions yield minimum reflection at the center frequency, 
it does not necessarily yield the most desired result. From an engineering perspective, 
once a certain reflection threshold is satisfied then it is a good enough solution. It 
would make sense then, to search for other possible solutions which simultaneously 
satisfy the reflection threshold and yield the widest bandwidth possible. Our effort 
in deriving the full mathematical expression, however, would prove to be useful since 
it enables us to automate extensive searches for other solutions. In doing so, we first 
try a simple cut-and-try algorithm which flowchart is shown in Fig. 3.1. The search 
space mentioned in the chart consists of equally-spaced 20 points centered around 
the maximally flat solution. The lower and upper limit of the search, however, are 
manually set and changed with each iteration. We repeat the search until the results 
stay fixed within 4 decimal points (default numerical precision in Matlab). 


This algorithm is the result of collaboration with William Kim, a junior at the 
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time in Pratt School of Engineering, who worked for a year as research assistant 
under supervision of Wiwi Samsul. The algorithm, while successfully discovers better 
solutions, but also has limitations in scope of how far it can go. To address this issue, 


we develop Genetic Algorithm to further extend the searching quality and flexibility. 








Load the 
maximally flat 
solutions 


Generate search space around 
the maximally flat solution 


Try one set of G values from the 
search space and calculate 
bandwidth 


Better bandwidth found? 







End of search space? 


Update bandwidth record 


FIGURE 3.1: Cut-and-try algorithm flowchart 


3.2 Components of Genetic Algorithm 


Genetic algorithm (GA) is an optimization algorithm which mimics the way natural 


selection process sorts out individuals who survive or do not survive by their fitness 
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[18, 19]. Through iterations of natural selection, we hope to find individuals which 
genes produces a desired fitness. 


In this work, the key components of GA are: 
e Population and individual 
e Each individual’s genes 
e Fitness criteria and selection threshold 
e Probability of gene mutation 
While the key processes include: 
e Generating initial population 
e Evaluating each individual’s fitness and selecting survivors 
e Mixing of survivor’s genes to produce new generation 
e Introducing mutation into the population’s genes 


In the following sections we describe the translation of GA key components and 


the implementation of each components on solving Jaumann wave absorber problem. 


3.3 Translation of GA components into Jaumann Wave Absorber 
problem 


3.9.1 Population 


We define the population as the set of all possible solutions. In addition, we present 
two types of solutions: solving only for the G values and solving for both G and X 
values. An individual consists of either a set of G values (G4, Go, ..., Gn) or G and X 
values (Gi, Go, … Gn, X1, X2, ..., Xn). Each G or X value corresponds to one gene. 
In other words, each individual in n layer case has n genes (if we only solve for G) 
or 2n genes (if we solve for both G and X). 
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3.9.2 Fitness criteria 


The fitness of each individual is defined as the relative bandwidth of the resulting I 


plot of its genes below a predetermined dB level (either -10 or -20 dB): 


Fitness = bah (3.1) 


fo 





fo 


center frequency 


fi = frequency when the first time I" plot dips below the dB level 
fa = frequency when the first time I" plot rises above the dB level 


after f; has been recorded 


We only include the first time I" plot dips below and rises above the dB level in 
order to ensure every individual which T plot oscillates at the dB level has low fitness 
and gets eliminated from gene pool. 

For the case of fixed quarter wavelength separation case, the I" curve is calculated 
using Eq. 2.14, 2.17, 2.20, A.14, A.15, and A.16 (with Zo and Zs are set equal to 1) 
for 1, 2, 3, 4, 5, and 6 layers case respectively. 

For the case of varying separation case, the I" curve is calculated using equations 


presented in Appendix A.4. 
3.9.3 Initial population 


The initial population is generated with uniform random distribution between value 
0 to 2. The limit was set as such because our previous maximally flat result shown 
in Table 2.1 provides reasonable ground to expect the optimal values for 1 to 6 layers 


case are between 0 and 2. 
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3.3.4 Survivor selection 


After each individual's fitness is calculated, the population is then sorted in de- 
scending order. Then we select the better half of the population as survivors. This 


threshold is arbitrarily chosen. 
3.9.5 Spawning next generation 


We randomly choose 2 individuals from the survivor pool as the parents and ran- 
domly select the genes from each parent to spawn a child. The probability of se- 
lecting each gene is 50% regardless of the parents’ fitness. The gene selection pro- 
cess also follows the proper sequence. For example, suppose parent A has gene 
Gia, Goa, Gza,..., Gna, parent B has gene G 1b, Gob, G3), … Gb, and the resulting 
child has gene Gc, Gac, G3c, … Gne. Then gene Gyc is randomly picked between 
Gra and Gb for all k € [1,n]. 


3.3.6 Mutation process 


The mutation process happens after the next generation is spawned. The rules of 


mutation are: 


e Mutation happens with equal probability (p) for each gene and independent of 


each event. We set probability of mutation p = 0.2. 


e Value of mutation (v) is determined by random uniform distribution with pre- 


determined boundary (b). In mathematical notation: |v| < b. We set b = 0.5. 


e If mutation occurs, the value of mutation is added algebraically into the gene. 


Otherwise, v = 0. 


The mutation process may cause a gene to turn into a negative value. We do 
not check for negative genes during the search process. However, if the algorithm 
produces final result with negative G or X values then that result is discarded. 
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3.4 Genetic algorithm flowchart 
Figure 3.2 shows the algorithm flowchart with additional specifications as follows: 


e Initial population is arbitrarily set to be 100. 
e Number of genes is equal to number of variables being searched. 


e Number of generation limit is arbitrarily set to be 200. This condition prevents 


infinite loop. 


e The saved population data comprises of the genes, relative bandwidth, band- 
width start, bandwidth end, theta limit, and the generation number when the 


best solution is found. 


e The fitness is evaluated at the specified dB level (either -10 dB or -20 dB in 


this work, but the algorithm can work for an arbitrary dB level). 
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Generate initial 


First iteration? H 
population 








Evaluate fitness 







Load saved 
population data 


Sort by fitness and 
determine survivors 


Record the best 
individual 








Randomly select 2 
parents 


Crossover parents” 
genes to produce a 
new individual 





Does the number of population have been 
replenished? 


Mutation process 


No 
Does the generation limit have 









been reached? 


Save results in a file 


FIGURE 3.2: Genetic algorithm flowchart 
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4 


Results and Analysis 


We structure this chapter by first laying out all the results then comparing and 


analyzing them in later in section 4.3. 
4.1 Cut-and-try Algorithm 


Table 4.1 presents the result from cut-and-try algorithm up to 5 layers case. The 
RBW increase is compared to the maximally flat solution. The bandwidth increase 
peaks at 3 layers case and then declines. We never managed to find the solutions 
for 6 layers case with this algorithm because of lack of computational power. The 
computers we use keep crashing in the middle of the searches. We think the main 
problem in this way of searching is because the algorithm has to keep track of every 
permutation of previous loops. Therefore, each additional layer increases the search 
space exponentially. On the other hand, GA does not encounter the same problem 
because older genes are erased and replaced with new genes on each generation. 
For the sake of brevity, we only show a sample plot of reflection coefficient from 


3 layers case. All other plots are available on Appendix A.8. 
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Table 4.1: Summary of the G values from cut-and-try algorithm up to five sections. 




















# of sections | Gi Go G3 Ga Gs | RBW | RBW increase (%) 
1 1.000 - - - - 0.253 0 
2 1.253 | 0.423 - - - 0.815 21.034 
3 1.368 | 0.514 | 0.316 - - 1.214 31.106 
4 1.500 | 0.605 | 0.290 | 0.290 - 1.402 28.691 
5 1.211 | 0.653 | 0.321 | 0.226 | 0.226 | 1.523 26.589 
































S4 comparison between maximally flat and cut-and-try solutions for 3-section case 
0 T T T T T T T 
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I —— Maximally Flat Solution 








V — — Cut-and-try Solution 
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FIGURE 4.1: Comparison of maximally flat and cut-and-try solutions for 3 section 
case. 


4.2 Genetic Algorithm 


We separate this section into two subsections: first subsection is comprised of results 
in the standard format like the preceding section and second subsection is comprised 
of results in expanded format while introducing a new variable we call ” theta limit” t 
(where £ = k fo, or the upper frequency limit of search range). The reason we intro- 
duce this new variable is because we realize that GA is able to efficiently search for 
solutions in any multiples of bandwidth range with respect to fo. In other words, for 
any given fo, GA can search for any solution which produce any relative bandwidth 


(RBW) we want by specifying the theta limit. 
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4.2.1 Standard Results 


The standard results of GA solution is shown in Table 4.2. The °% increase’ column 
shows the RBW increase compared to the cut-and-try algorithm. The two right 
most columns are the relative bandwidth start and end. In addition, we provide the 
start and end of bandwidth to highlight that the standard solutions only searches 
for bandwidth around fo while the expanded search (described below) can search for 
any bandwidth even outside of fo. 


Table 4.2: Summary of standard GA solutions (51; cut-off: -20 dB). 


























# G Go G3 G4 Gs Ge | RBW | 96 increase | BW Start | BW End 
1 | 1.016 - - - - - 0.254 0.416 0.873 1.127 
2 | 1.256 | 0.426 - - - - 0.872 6.963 0.565 1.436 
3 | 1.295 | 0.532 | 0.305 - - - 1.218 0.284 0.391 1.609 
4 | 1.365 | 0.593 | 0.302 | 0.274 - - 1.408 0.432 0.296 1.704 
5 | 1.250 | 0.637 | 0.347 | 0.217 | 0.238 - 1.525 0.098 0.238 1.763 
6 | 1.221 | 0.610 | 0.399 | 0.256 | 0.164 | 0.222 | 1.604 N/A 0.198 1.802 






































4.2.2 Expanded Results 


We take results from 3 layers case as an example of the expanded solution and show 
them in Figure 4.2 and Table 4.3 below. In the last column, "T" refers to "Theta 
limit”. As usual, we put the rest of the results in Appendix for the sake of brevity. 
Note that the ability of GA to search for solutions outside of fo gives us a huge 


flexibility in designing a system with any bandwidth in any frequency range. 
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FIGURE 4.2: 51 plot of GA expanded result for 3 sections case (theta limit = 20). 


Table 4.3: Summary of expanded GA solutions for 3 sections case (511 cut-off: -20 
dB). 

































































Gy Ga G3 Xi X» X3 | RBW | BW Start | BW End | T 
1.377 | 0.533 | 0.321 | 0.784 | 0.791 | 0.813 | 1.506 0.495 2.000 1 
1.365 | 0.508 | 0.313 | 0.402 | 0.405 | 0.400 | 3.015 0.976 3.991 2 
1.410 | 0.576 | 0.317 | 0.253 | 0.265 | 0.262 | 4.467 1.533 6.000 3 
1.230 | 0.541 | 0.273 | 0.201 | 0.198 | 0.192 | 5.995 2.005 8.000 4 
1.376 | 0.573 | 0.273 | 0.155 | 0.151 | 0.164 | 7.359 2.641 10.000 5 
1.385 | 0.480 | 0.300 | 0.131 | 0.129 | 0.125 | 8.907 3.093 12.000 6 
1.522 | 0.489 | 0.261 | 0.107 | 0.114 | 0.109 | 10.199 3.802 14.000 7 
1.462 | 0.507 | 0.309 | 0.102 | 0.100 | 0.095 | 11.935 4.065 16.000 8 
1.181 | 0.518 | 0.252 | 0.088 | 0.087 | 0.087 | 13.444 4.556 18.000 9 
1.263 | 0.518 | 0.284 | 0.077 | 0.076 | 0.076 | 14.836 5.164 20.000 | 10 
1.280 | 0.496 | 0.273 | 0.065 | 0.066 | 0.064 | 15.906 6.095 22.000 | 11 
1.270 | 0.499 | 0.274 | 0.068 | 0.067 | 0.064 | 18.030 5.970 24.000 | 12 
1.281 | 0.523 | 0.273 | 0.061 | 0.062 | 0.063 | 19.485 6.476 25.961 | 13 
1.428 | 0.561 | 0.334 | 0.057 | 0.057 | 0.057 | 20.942 7.059 28.000 | 14 
1.313 | 0.513 | 0.261 | 0.051 | 0.052 | 0.057 | 22.272 7.690 29.962 |15 
1.445 | 0.497 | 0.266 | 0.049 | 0.047 | 0.053 | 23.475 8.525 32.000 | 16 
1.357 | 0.487 | 0.270 | 0.045 | 0.047 | 0.043 | 25.058 8.942 34.000 | 17 
1.410 | 0.548 | 0.341 | 0.042 | 0.042 | 0.043 | 26.718 9.283 36.000 | 18 
3.654 | 0.964 | 0.380 | 0.014 | 0.034 | 0.038 | 23.548 14.452 38.000 | 19 
0.541 | 1.229 | 0.417 | 0.007 | 0.028 | 0.036 | 23.780 16.006 39.786 20 
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4.3 Results Comparison and Analysis 


4.5.1  Gain-Bandwidth Product 


Gain-bandwidth product (GBWP) is the product of gain and bandwidth at which 
the gain is measured. In this work, we use this attribute to compare our result with 
previous works particularly from ref. [20] and [13]. Since our goal is to produce 
bandwidth as wide as possible at the desired gain, the closer the plot to upper left 
corner the better the result. All calculations of GBWP in this section is performed 
at fo = 250 MHz. 

First, we show the comparison of maximally flat solutions from 1 to 6 layers in 
Fig. 4.3. As the number of layers increases, the the plot moves closer to the upper 


left corner which indicates a bigger bandwidth for the same range of frequency. 


Gain-bandwidth product vs dB level (maximally flat response) 
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FIGURE 4.3: Gain-bandwidth product comparison between maximally flat solutions. 


Next, Fig. 4.4 shows the comparison between maximally flat solution with GA 


expanded result (theta limit set to 1) at 3-layer case. 
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FIGURE 4.4: Gain-bandwidth product comparison between 3 layer maximally flat 
solution and GA. 


Finally, Fig. 4.5 shows the sample comparison of GBWP between our result (GA 
expanded result and theta limit set to 1) and results from ref. [20] (Kraus-Fleisch in 
the plot legend) and [13] (Nortier in the plot legend). All comparisons are performed 


at 3-layer case. 
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FIGURE 4.5: Gain-bandwidth product comparison at 3 layers between our work 
(GA), ref. [20], and ref. [13]. 
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4.5.2 Comparison with other works. 


Other works on Jaumann absorbers can be found in ref [21, 17, 15, 13, 12, 22, 23, 24, 
25, 26, 27, 28] while optimization of radar absorbers with genetic algorithm includes 
ref [24, 29, 30, 31, 32, 26, 33]. However, only ref [24, 31, 26] involve both GA and 
Jaumann absorbers so we will compare our work with the latter three works. 

Table 4.4 lists our result with Chambers et al [24]'s results for one, two, and three 


layers case. "NBW" column stands for " Normalized Bandwidth", which is NBW = 





24u p The fact that both our results has the same NBW suggests that they are the 
maximum NBW that can be achieved in the respective number of layers. However, 
Chambers et al's GA still uses the default quarter wavelength spacings while our GA 
varies the spacings and finds solutions which give shorter design. A shorter design is 
considered better to save space and materials. 

Both ref [31] and [26] incorporate spacing variation into their GA. However, result 
from ref [31] achieves 95% absorption at maximum while all our results achieve 99% 
absorption (-20 dB) which is usually the standard absorption level. Furthermore, 
the way our GA is implemented makes the algorithm flexible and can be set to 
search for an even higher absorption level. Finally, ref [26] maintains to achieve 
absorption level of -35 dB with total thickness of 2.10 cm (or 7.10 A as indicated by 
their classic Jaumann comparison) in 5 layers case. We set our GA to search for the 
same absorption level and produce the following gene: [0.46, 0.89, 1.81, 1.24, 0.40, 
2.50, 1.69, 1.67, 1.64, 1.55]. The first five genes are the G values while the latter five 
genes are X values which is a multiple of quarter wavelength and gives us a total 
length of 2.26 A (or 68% shorter). 

Finally, we also compare our result with ref [34] in more detail. Two takeaways 
from this reference are that any results in normal incidence perform best compared 


to any other angle, and any result that is specifically optimized in specific incident 
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angle perform best when subjected to normal angle. Thus, we only need to compare 
our results in normal incidence scenario. While their work include spacing variation, 
all of their results state the spacing with only a single number which must mean all of 
the spacings have equal length and they can not mean the total length because their 
result would have become absurdly improbable while taking it as spacing between 
each layers yields results which makes sense. We will compare our result with their 
best result at 4-layers case as shown in Table 4.5. Since their report only provides 
figures instead of numerical value, we have to approximate (generously) their NBW 
based on their reflection coefficient plots. Their best 4-layers result has spacing of 
7.15 mm each with e, = 1.1. Assuming fy = 10G Hz (their center frequency plot), 
this results in a thickness of 1.0006 A with NBW (generously) approximated at 1.40. 
In contrast, our best 4-layers result is [2.40, 0.90, 0.44, 0.28, 0.40, 0.73, 0.77, 0.77] 
with total thickness of 0.67 A. It appears varying the spacings while limiting them all 
to be equal would yield results that are not far from the default quarter-wavelength 
ones. In contrast, our GA is able to vary all of them independently of each other 
which results in much bigger search space. 


In conclusion, our work produces results which are: 

e 44% /28%/20% shorter in design length than ref [24], 
e at least 4% higher absorption level than ref [31], 

e 68% shorter in design length than ref [26], and 


e 33% shorter in design length than ref [34]. 


36 


Table 4.4: Comparison with ref [24]. 


















































Ref [24] Our work % Shorter 
NBW | Length (A) | NBW | Length (A) 
l layer | 0.25 0.25 0.25 0.14 44 
2 layers | 0.87 0.50 0.87 0.36 28 
3 layers | 1.22 0.75 1:21 0.60 20 
Table 4.5: Comparison with ref [34]. 
Ref [34] Our work % Shorter 
NBW | Length (A) | NBW | Length (A) 
4 layers | 1.40 1.00 1.24 0.67 33 




















37 








5 


Measurements 


5.1 Realization in transmission line 


We experimentally verify our 3-layer result for the Maximally Flat Solution (MFS) 
in transmission line setting. First we build a 50 (2 transmission line as shown in Fig. 
5.1, which —for the sake of brevity—calculation details and dimensions are provided 
in Appendix A.10. Then we pick the available off-the-shelf resistors which values are 
the closest to the theoretical solutions as listed in Table 5.1 and build the Jaumann 
absorber circuit. We first perform the experiment using lumped resistors, but the 
result is distorted such that the Su plot fluctuates around the -20 dB line (dotted 
line in Fig. 5.2). We hypothesize that stray inductances is the main culprit for the 
distorted result, so we simulate the theoretical solution with added 8 nH to all the 
resistors and we obtain good agreement with the experiment (Fig. 5.2). In addition, 
we also redo the experiment by replacing the lumped elements with chip resistors in 
order to obtain an improved result as shown in Fig. 5.3. 

In conclusion, we have experimentally proven a sample solution of Jaumann wave 
absorber problem. Any difficulty in real life implementation mainly revolves around 
producing precise resistive value either in transmission line model or wave absorber 
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model. 


Table 5.1: List of resistor values used in 3 layers MFS experiments. 

















RBW 
SCHENTEN 
Theoretical | 30.414 | 97.220 | 458.295 1.356 0.926 
Lumped 30 100 470 1.523 N/A 
Chip 30.1 100 475 1.625 0.975 
































FIGURE 5.1: Experimental set up for three sections Jaumann wave absorber in 
transmission line. 





0 250 500 
f (MHz) 


FIGURE 5.2: Comparison between simulation with added stray inductances (solid 
line) and measured 511 (dotted line) for 3 section case. 
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Experiment Results of 3 Layers Maximally Flat Solution Jaumann Wave Absorber 
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FIGURE 5.3: Du comparison between using lumped vs chip resistors. 


5.2 Measurement of ambient RF energy 


As part of examining the viability of Jaumann wave absorber on energy harvesting 
technology, we measure the availability of ambient energy that can be harvested 
especially in RF frequency. We perform this measurement of ambient RF power 
in four Duke buildings: CIEMAS, Teer, Physics, and French. Energy detection 
is performed with two dipoles, bowtie, and spiral antennas connected to spectrum 
analyzer. We record the spectral energy density in various locations throughout the 
buildings. We categorize the locations as: corners, halls, rooms, near router/repeater, 
and near window. Note that some locations fall into more than one category (e.g: 
corners which have routers nearby). Experiment details and results are summarized 


in the following subsections. 
5.2.1 Locations Stats 


Table 5.2 lists the summary of measurement locations. More details on these loca- 


tions is provided in Appendix A.11. 
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Table 5.2: Summary of measurement locations. 




















a Category 
Bulling Corner | Hall | Room | Near Router | Near Window "sata Bol 
CIEMAS 3 2 2 4 2 13 8 
Teer 2 0 3 4 2 11 9 
Physics 4 6 1 10 1 22 12 
French 7 8 3 7 3 28 18 
































5.2.2 Antennas Profiles 


Fig. 5.4 shows the 511 profile of each antenna used in the experiment. 


S-parameter vs frequency 
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FIGURE 5.4: 84, graphs of antennas used in experiment. 


5.2.3 Spectrum Analyzer 


We measure the ambient spectrum with Spectrum Analyzer from Agilent Model 
E4446A (serial number MY46180443) which settings of the analyzer is detailed in 
Table 5.3. A sample of spectrum measurement taken at CIEMAS 3rd floor skyway 
is shown in Fig. 5.5. We pick this spectrum to show because it has the highest 
detected peak (at -21.22 dBm around 0.9 GHz) compared all other measurements 


even surpassing the ones taken near router/signal repeater. The spikes around 0.9 
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GHz indicate presence of cellphone signal while the spikes around 2.4 GHz indicate 
presence of WiFi. We hypothesize this is because of the skyway lacks concrete walls 


which is the main culprit of attenuation of cellphone signals. 


Table 5.3: Setting details of the spectrum analyzer. 





Min. Frequency 0.01 GHz 
Cent. Frequency 2.01 GHz 
Max. Frequency 4.01 GHz 
Resolution bandwidth | 3 MHz 
Video bandwidth 3 MHz 



































Sweep time 38 

Num of data points 601 
Ref level 0 dBm 

Data format „CSV 
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FIGURE 5.5: Sample of a spectrum taken in CIEMAS 3rd floor skyway. 


5.2.4 Ambient Energy Calculations 


Our goal is to obtain a figure of power-to-area of ambient field which Eq.5.1 and 
5.2 show the general calculation process. However, the details of the calculation is 


not as straightforward as it seems. Referring to table 5.3, we notice that while the 
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resolution bandwidth is 3 MHz, the spectrum analyzer records the measurement in 
601 data points over the span of 4 GHz which means there is a 6.6 MHz gap between 
data points. In order to counter this mismatch, we first interpolate the spectrum to 


match the resolution bandwidth of 3 MHz gap between data points. 





Feng 

Pk Te 1051/10 (5.1) 
Ind Puma 
` = ` 5.2 
unit area Ae 2 
where: 
Pamb. = ambient power (W) 
Peas. = measured power (W) 


Su = reflection coefficient (dB) 


A. = effective area of antenna (cm?) 


Furthermore, we need to obtain the effective area of the spiral antenna which does 
not have a straightforward formula to calculate. Hence, we opt to simulate the spiral 
antenna in CST and obtain its maximum gain then calculate the effective area using 
Eq. 5.3. Table 5.4 shows the parameters of the spiral antenna. Inner diameter refers 
to the spacing between spirals or in other words the difference of diameter between 
a spiral and its subsequent spiral. Outer diameter refers to the total diameter of the 


outermost spiral. 


M 
A EG (5.3) 
AT 
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Table 5.4: Parameters of spiral antenna. 


























Name Description Value 
D; Inner diameter 1.5 mm 
D, Outer diameter 143 mm 
N Number of turns 96 
h Handedness Right handed 








Putting all of them together, we obtain the average power spectrum across all 
buildings and all location categories as shown in Fig. 5.6. We also integrate the 
spectrum to obtain the total power between 0-4 GHz as 0.9263 wW/cm?. For com- 
parison, other studies conducted in Belgium [35], UK [36, 37], Netherlands [38], and 
Japan [39] on ambient RF spectrum maintain that the energy level varies between 
0.3 uW/cm? - 3.6 mW/cm? where the lower and upper limit were obtained from 
measurements in rural and urban area respectively. Those other studies are implied 
to be conducted in outdoor environment while our measurement is in indoor setting. 

Detailed breakdowns of average ambient energy by location categories and build- 
ings are shown in Table 5.5. Surprisingly, ‘corner’ sites has the highest concentration 
of RF energy, almost 5 times higher compared to "near router’ sites. In addition, 
‘hall’ also has higher energy than "near router’. Combined with the fact that ’room’ 
has the lowest energy, this phenomenon suggests that ambient energy tend to con- 
centrate near walls (despite the definition of ’room’ as space contained within 4 walls, 
none of the rooms can be considered claustrophobic; width of the rooms tend to be 
wider than halls). It is also worth noting that ‘near window’ has higher energy con- 
centration than "near router’. We suspect it is because of contribution from cellphone 
signals from base stations which takes attenuation only from a window compared to 
multiple walls in halls’ or rooms”. Finally, Table 5.6 shows the average by buildings. 
Since category of "corner" holds the most concentrated energy, we take a closer look 


at the proportion of corners in each building. CIEMAS has 2 out of 11 (18.18%), 
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Teer has 2 out of 11 (18.18%), Physics has 4 out of 21 (19.05%), and French has 7 out 
of 26 (26.92%). Physics and French have greater proportion of corners compared to 
Teer yet Teer still holds a higher ambient RF energy concentration. Therefore, it is 
not because of the corners that makes Teer and CIEMAS have higher concentration 
than Physics and French. We hypothesize the reason is due to the number of stu- 
dents around. Both Teer and CIEMAS have dedicated sites for students to work and 


the student's laptops and cellphones increase the ambient RF energy concentration. 
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FIGURE 5.6: Average power spectrum across all buildings and all location categories. 


Table 5.5: Summary of average ambient energy by location category. 




















Category Average ambient energy (uW Zem?) 
Corner 5.2031 
Hall 1.1915 
Room 0.2877 
Near router 1.1120 
Near window 1.4031 
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Table 5.6: Summary of average ambient energy by buildings. 

















Building | Average ambient energy (uW /cm2) 
CIEMAS 1.2234 
Teer 1.6808 
Physics 0.3576 
French 0.8480 
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6 


Conclusions and Future Work 


6.1 Conclusions 


We have briefly explored the origin of Jaumann wave absorber and provided a solid 
mathematical framework to solve for its maximally flat solution for any number of 
layers while past attempts to do so were merely blind trial-and-error [8], stopped 
at 2 layers [17], impractical to go beyond 3 layers [14], did not produce a reliable 
absorption spectrum [13], did not optimize for absorber spacings [21, 17, 15, 13, 12, 
22, 23, 24, 25, 27, 28], or still has room for improvement [21, 31, 26]. In comparison 
to the last three references, we surmise that our GA gives better results because the 
way we design our GA by starting from the framework of maximally flat solution 
and automate the search for other solutions seem to provide a greater flexibility 
in searching with any absorption level, varying both resistances and spacings, and 
finding a more, if not the most, optimum solution compared to other works. 

In order to support the feasibility of building a wave absorber system, we perform 
measurements in four Duke buildings and compare our result with other studies. We 


conclude that the average ambient RF energy in indoor environment ranges between 
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0.2877 - 5.2031 uW/cm?. We also conclude that the best spots to install wave 
absorbing device would be in corners since it is the category in which the highest 


ambient RF energy concentrates. 
6.2 Future Work 


This work sees possibility of usage for constructing energy harvesting devices through 
absorbing electromagnetic wave. By implementing the right values of resistances and 
their spacings, one can engineer the absorption level and bandwidth to suit one’s 
needs. However, given the measurement results of ambient RF energy in buildings, 
such devices must operate in uW level unless we artificially pump the ambient RF 
energy, effectively turning the whole system into wireless charging. Another possible 
usage comes in radar cloaking devices. By engineering the wave absorber to absorb in 
targeted bandwidth with certain absorption level, one could render radars ineffective 
in detecting the cloaked item. This idea offers another approach to the existing 


stealth systems. 
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Appendix A 


Appendix List 


A.1 System of equations to obtain solutions for four, five, and six 
sections case 


A.1.1 Four sections 
The ABCD matrix for four sections case is: 
1 
A=; ( — AG3G3Z2 cos 20 + 2(4 + G3G3Z2 + G (Ga + G3) Z2) cos 40 + 


Zo(—2((—G2)G3 + G. (G + G3))Zo ES 2j(2G3 + Gi(—2 + Catia Sol sin 20 + 





Boy, sin 9 (4j cos? 0 — (3G, + 4G2 + 3G3) Zo cos? 0 sin 0 — 





2j(2 + Catia Zä + Gi(G3 + G3)Z3) cos 0 sin? 0 + Zo(Gi + G3 + G1G2G3Z2) sin? 0) 





1 
C= 87, (D'Sich + G3 = G3G3G4Z2) cos 20 + Zo(AG» | 4G3 | SCH, | 2GG3G4Z2 + 
0 


G4 (4 + 2G3G4Z3 + G2(G3 + 2G4)Z2)) COS 40 — j((—9)(2G1G3G4Z + G2Z0(G1(G3 + 








D = cos? Ø | (Gi 2G» | 3G3 | 4G4) Zo cos? Ø sin — (6 + 2G2G3Z2 + 3G3G4Z; + 





Gi (Ga + 2G3 np 3G4) Zé) cos? 0 sin? 0— 7 Z0(G3 + 4G4 + 2Ga2(1 + G3G4Z2) + 
G1(3 + 2G3G4Z3 + G2(G3 + 2G4)Z3)) cos 0 sin? 0 + 


(1 + G3G4Z; + G1Z ¿(Ga + Ga + G2G3G4Z2)) sin* 0 — G2G4Z; sin? 20 (A.1) 
The load Zr with short circuit load is: 


AZ, LP B 
= = A2 
Zr CLA D (A2) 
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Setting Zo = 1 and Zr = 1 and cross multiply, we obtain: 





j cos? 0(—4 + Gi + 2G» + 3G3 + AG4 — J cot 0) — 





(6 + 3G3(—1 + G4) + 2G5(—2 + Gs + 2G4) + G1(—3 + G> + 2G3 + 3G4)) cos? 0 sin 0 — 


J(=4 + Gs + 2G2(1 + Gall + G4)) + 4G4 








G4(3 + 2G3(—1 + Ga) + Ga(—2 + G3 + 2G4))) cos 0 sin? 6+ 
(1+ G3(—14 G4) +G (—1 + Goll + G3(—1 + G4)) + G4)) sin? 0 = 0 


(A.3) 
Therefore, the system of equations that we need to solve are: 


—4 + G; + 2G, + 3G3 + 4G, = 0 

6 + 3G3(—1 + G4) + 2G2(-2 + Gs + 2G4) + Gi(—3 + Ga + 2G3 + 3G4) = 0 
—4 + G3 + 2G2(1 + G3(—1 + G4)) + 4G4 + G1(3 + 2G3(—1 + G4) +... 
G3(—2 + Ga + 2G4)) =0 
(1+G3(=-1+G4)+Gi(-1+ + G,(—1+G,)) = 0 














(A.4) 


And the solution is: G4 = 1.7824, G> = 0.6695, G3 = 0.2337, G4 = 0.0443. 
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A.1.2 Five sections 


The ABCD matrix for five sections case is: 


A 


cos? 0 + J(4G1 3G» | 2G3 | G4) Zo cos? 0 sin 0 — (10 + Ga G4Z; + 2G2(G3 + G4)Z2 4 














G1(3G» + 4G3 + 3G4) Z2) cos? 0 sin? 0— JZo(6(G3 | G4) | Ga(4 | G3G4 Z8) + 2G, (2 4 
G3GLZ2 + G2(G3 + G4)Z2)) cos? 0 sin? 0 + (5 + 3G3G4Z; + 2G2(G3 + G) Zi + 


G1Ze(Ga + Ga + G2G3G4Z3)) cos 0 sin* 0 + J ZÁ (G> + Ga + G2G3G4Z3) sin? 0 





Ži sin (5j cost Ø — 2(2G, + 3G, + 3G3 + 2G4) Zo cos? Ø sin 0 — 34 (G3G4 + 


Gi(G3 + Ga) Zé cos? 0 sin? 0 + 2Z0(Ga + G3 + 2G4 + G2G3G4Z; + 
G1(2 + G3G4Z; + G2(G3 + G4)Z2)) cos 6 sin? 0 + (1/2)5(2(1 + G3G4Z + 


G1 Z2(Gy + Ga + G3G3G422)) sin? 6 — (5 + 2G1G322 + 2G2(G3 + G4)Z2) sin 26°)) 
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1 
m (2Z0(G2(4 + G Zå(-2G4 — 2Gs + G3GaG5Z2)) — 2(G1G4G5Z2 + Ga(—2 — 
0 


G4GsZ$ + G (Ga + 2G5) Z2))) cos 0 — Z9(AGG3G 5 Zi + 2G4(—4 + 4GsG;Z + 
Go(G3 + 2G5) Z2) + G,(—2(4 + G3 (Ga + 2G5) Zó) + 

G3Z5(—2G4 + G3(2 + 3G1G543)))) cos 30 + 

8G Zo cos 50 + 8G2 Zo cos 50 + 8G3 Zo cos 50 + 8G4Z cos 50 + 16G5 Zo cos 50 + 


2G1G3G32Z3 cos 50 + 2G1G3G4Z3 cos 50 + 2G1GaG4Z3 cos 50 + 2G3G3G4Z3 cos 50 + 








4G1G3Gs5Z$3 cos 50 + AG 1 G3G5Z2 cos 50 + AG2G3G5Z3 cos 50 + AG 1 GAG5Z3 cos 50 + 








AG>G1G523 cos 50 + AG3G1G523 cos 50 + G1G3G3G4G5Z; cos 50 — Belle sin 0 + 
8jG3G3Z, sin 0 + 4jG1G4Z,; sin 0 — 44G3G4Z; sin 0 + 8jG3GsZ; sin O — 

SJG3G5 Zê sin 0 — 27G1G2G3G1Z2 sin O — 

4jG1G3G3GsZq sind — 43G1G3G4GsZq sin 0 — 43G41G3G4GsZg sin 0 + 
8jG3G3G4G5Z, sin 0 + 4jG1G3Z; sin 30 — 4jG3G3Z; sin 30 + 

8jG1G4Z, sin 30 + 4jG3G4Z, sin 30 + 8G1GsZ; sin 30 — 8jG4G5Z, sin 30 — 

jG, G2G3G4Z, sin 30 — 2jG1G3G3G5Z, sin 30 — 2;G1G3G4G5Z, sin 30 — 
2jG1G3G4G5Z, sin 30 — 6jG2G3G4G5 Zo sin 30 + 16j sin 50 + 4jG1G3Z, sin 50 + 
4jG1G3Z; sin 50 + 4jG3G3Z; sin 50 + 4JG1G 473 sin 50 + 


4jG3G4Z; sin 50 + 4jG3G4Z; sin 50 + 8jG1G5Z; sin 50 + 











8jG3G;Z, sin DÉI + 8jG3G5Z, sin 50 + 8jG4G5Z, sin 50 + ¡G1G2G3G4Z, sin 50 + 
2jG1G3G3G5Z, sin 50 + 2;G1G3G4G5Z, sin 50 + 2;G1G3G4G5Z, sin 50 + 
24G3G3G4Gs Zi sin 50) 


(A.5) 


53 


D = cos 0 + (G; +2G2+3G3 + 4G4 + 5G5) Zo cos" 0 sin 0 — 





kaa TG 


2G2(G3 + 2G4 + 3G5) Zé + G, (G | 2G3 | 3G4 | 4G) Z2) cos? 0 sin? 0 — 





jZo(4G3 + AG4 + 10G5 + 3G3G4G5Z; + 2G2(3 + 2G4G5 Zê + Ga(G4 + 2G5) Zo) + 
G, (6 + 3G4Gs5Z; + 2G3(G4 + 2G;)Z + G2(G3 + 2G4 + 3G5) Ze) cos? 0 sin? 0 + 


(5 + G3G1Z§ + 2G3G5Z6 + AG1G5Z3 + 2G2Z3 (G3 + Gs + G3G1G5Z2) + 





G1 Z5(2G3 + G4 + 4G5 + 2G3G1G5Z2 + 

G»(3 + 2G GsZ + Ga(G4 + 2G5) Sall cos 0 sin? 0 + 

jZo(G3 + Gs + G3G4G5Z; + G (1 + G ,GsZ + GoZ3 (G3 + Gs + G3G1G5 43))) sin? 0 
(A.6) 


The load Zr with short circuit load is: 


AZL+B B 
Zr = = 
CZ,+D D 
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Setting Zo = 1 and Zr = 1 and cross multiply, we obtain: 


cos? 6(—5 + G4 + 2G2 + 3G3 + 4G, + 5G; — j cot 0) + 








j(10 — 6G; — 4G, + 3G3G4 + 6G3G5 + 4G4G5 + 


2G2(-3 | G3 | 2G4 3G5) | Gi(—4 | Ga | 2G3 3G4 | 4G5)) cos? O sin 0 — 

















(UE (SiGe) EE E NE 


NNN SEN 








G3(—3 + G3 + 2G4 + 3G3))) cos? 0 sin? 0 — 3(5 + G3(-2 + G4) = 4G 4 + 


2(G3 + 2G4)Gs + 2G2(—1 + Gs — GaGa + Gs + G3G4G5) + Gi(—4 + G4 + 2G3(1 + 





G4(—1 | Gs)) | AC, | Go(3 | Gs(—2 + G4) — 2G4 + 2(G3 + G4)G5))) cos 0 sin? 0 + 








(-1 + G3 — Gill + Goll + G3(—1 + G4)) + G4) 
G3G4(—1 + Gs) + Gs + G1(G2 + Ga + G2G3G4)G5)sint9 = 0 


(A.7) 
Therefore, the system of equations that we need to solve are: 


—5 + G4 + 2G2 + 3G3 + 4G, + 5G5 = 0 

(10 m 6G; = 4G 4 + 3G3Ga4 Sp 6G3G5 + 4G4Gs; "pes 

2G5(—3 + G3 + 2G4 + 3G5) + G4(—4 + G> + 2G3 + 3G4 + 4G5)) = 0 
(G3(4 + 3G4(—1 + Gs)) + 2(—5 + 2G, + 5Gs) + 2G2(3 + 2G,(—1 + Gs) 
G3(—2 + G4 + 2G5)) + Gi(6 + 3G4(—1 + Gs) +... 

2G3(-2 | Ga | 2G5) | Ga(-3 + G3 + 2G4 + 3G5))) =0 

(5 + G3(—2 + G4) — 4G4 + 2(G3 + 2G4)G5 +... 

2G3(—1 + Gs — G3G1 + Gs + Ga3G4Gs) +... 

G4(—4 + G4 + 2G3(1+ G4(—1 + G5)) +4G5+... 

Go(3 + G3(—2 + Ga) — 2G4 + 2(G3 + G4)G5))) = 0 

(-1+G3 — Gill + Goll + G3(—1+ G4)) + G4) +... 

G3G4(—1 + Gs) + Gs + G1(G2 + Ga + G2G3G4)G5) = 0 












































(A.8) 


And the solution is: Gi = 1.8679, Gy = 0.7777, Gs = 0.3427, G4 = 0.1136, 
Gs = 0.0188. 
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A.1.8 Six sections 


The ABCD matrix for six sections case is: 


A = cosf0 + 3(5G1 + 4G> + 3G3 + 2G4 + G5)Zo cos? Ø sin 0 — (15 + 2G3G4Z2 + 








2G3G5Z3 + G4GsZe + 2G1(2G2 + 3G3 + 3G4 + 2Gs)Z$ + 
G>(3G3 + AG 4 + 3G5) Zó) cos? 0 sin2 0 — J(GsGuG5 + 


GONE 





3G5))) Za cos> 0 sin? 0 + (15 T 6G3G4Z3 + 6G3G5 Zó + 6G4G5Z; T 














G3Ze(4(G4 + Gs) + Ga(4 + G4GsZ2)) + 2G1Z2(Gs + G4 + 2Gs + GsG4GsZ; + 


Cat? + G4GsZ; + G3(G4 + G5)Z6))) cos? 0 sin“ 0 + j Z9(3G3 | 2G4 | DCH, | 





3G3G4GsZe + 2G2(2 + Gate Ee + Ga(G4 + G5) ZG) + G (1 + G4GsZ$ + GaZe (G 3 + 
Gs + G3G4G5Z6))) cos 0 sin? 0 + (1/4)(—4(1 + G, Gs Zê + G2Z3(G3 + Gs + 


G3G4G523)) sinf 0 — j(5G1 + 4G2 + 5G3 + 6G4 + 5G5)Zo sin? 20) 





B = Zysin 0(6j cos? 0 — (5G + 8G2 + 9G3 + 8G4 + 5G5) Zo cos“ 0 sin O — 








2j(10 + 3G3G4Z2 + 3G3G 528 + 2G4GsZ; + G1(2G2 + 3G3 + 3G4 + 2G5) Z3 + 


G>(3G3 + 4G4 + 3G's) Ze) cos? 0 sin2 0 + Zo(8G» | 6G3 | SCH, | 10G; | 





3G3G4GsZ, + G1(10 + 3G1G5 E + 3G2(G3 + G5) Z2)) cos? 0 sin? 0 + 
2j(3 + G3G4Z2 + G3G5Z2 + Mats Zé + G228 (G 3 + Gs + G3G4Gs Ze) + 
Gi Zo (G3 + G4 + 2Gs + GaG4GsZ2 + G3(2 + G4Gs Zë + 

G3(Ga + G5) Z2))) cos 0 sint 0 — Zosin 0((G3 + Gs + GaG4GsZg + 

Gi(1 + GiGsZ + GoZ3 (G3 + Gs + G3G4G5Z2))) sint 0 — 

(G (GoGs + G3(Ga + Gs)) + Go(GuGs + G3(Ga + Gs))) Z2 sin? 20)) 


(A.9) 
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(IC: | Go | G3 | G4 | Gs | Gel Zo cos? 0 + j(6 + G3G4Z + 2G3G5Zj + 


G1G5Z3 + 8G3G6Z6 + 2G1G5Z3 + GsGeZg + Go(G3 + 2Ga + 3Gs + 4G6)Z2 + 











G. (G | 2G3 | 3G 4 | 4G5 | Bel Zê) cos° 0 sin 0 = Zo(6G3 | 7G4 | 10G; | 
15G6 + G3G1G5Z3 + 2G3G1G56Z3 + 2G3G5G56Z3 + G1G5G523 + Ga(7 + 


3G5GsZ; + 2G (Gs + 2G5) Z? + Ga(G4 + 2Gs + 3G6) Zo) + G1(10 + 3G Gs Z2 + 





DC Gs Zo + 4G5G6Zo + 2Ga(G4 + 2G; + 3G6) 23 + 


Ga(G3 + 2G4 + 3G5 + 4G5) Z2)) cos* 0 sin? O — 





JZ) (4G3G4 + 6G3G's + 6G4Gs + G3G1G5G6Z? + Go(4Gs + 2G4(2 + G5G6Z?) + 
G3(4 + 2G5GeZ2 + G (G; + 2Gs) Z2)) + G1(AG4 + 4G; + 3G4G5G6Z? + 
2G3(3 + 2G5G6Z) + Ga(Gs + 2G6) Za) + Go(6 + 3G5G6Z3 + 2G4(Gs + 2Gs) Za + 


G3(Ga + 2G5 + 3G6) Z3))) cos? 0 sin? 0 + Zo(9G3 | 7G4 | 5Gs | 15Gg | 





3G3G4GsZ; + 6G3G4GoZ, + 6Ga3GsGeZg + 6G4G5G6Zo + Go(7 + 4G5G6Zó + 


2G4(Gs + 2G6) Zé + Gs Ze (2G4 + 2G'5 + 4G 4 G4G5G6Z6)) + G1(5 + GaG5Z3 + 





2G4G6Zó + 4G5G5Z3 + 2G3Z2(G4 + Ge + G4G5G6 Zo) + GaZ3 (2G4 + Gs + 





4Gg + 2G4GsGeZg + G3(3 + 2G5G6Z; ar G4(Gs + 2Gs)Z2)))) cos? 0 sin? 0 + 








4(6 + 3G3G 423 + GuG5Z3 + 3GaGsZ; + 2G4GeZa + 5GsG6 Z + 3GaG4GsGgZy + 
G5Z2(2G4 + Gs + AGs + 2G4GsGeZa + G3(3 + 2G5G6Zó + G4(Gs + 2Gs) Z2)) + 


Gy Zo (Ga + Gs + G4G5G6Zo + Gall + GsGe Zi + 














G3Zó(Ga Gs G4G5G6Z6)))) cos 0 sin? 0 + (1/4)(—4Z0 (Gu + 


Gs + Gatistie Zo + Goll + GsGs Z + G3Z ¿(Ga + G + GaGsGeZZ))) sinf 0 — 





3010 + 5G1G6Z; SIE 4G3GsZ; + 5G3G6 Zo + 6G4GsZ; + 5G'sG'gZp) sin? 20))/ Zo 
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D = cos? 0 | (Gy | 2G» | 3G3 | 4G, | 5G; | 6G) Zo cos? O sin 0 — (15 + 3G3G4Z; + 





6G3G5Z5 + AG4GsZ; + 9G3GeZ, + 8G1G5 Za + 5GsGeZg + 


Ses IGJEN AO 








G, (Ga + 2G3 + 3G4 + 4G5 + 5G6) Zo) cos! O sin? 0 — 


J(3G3G4G5 + 6G3G1G, + 6G3G5Gg + 4G4G5G6 + 





VONG LIEN Gr ESN 











2G3(G4 + 2G5 + 3G6) + Go(G3 + 2G4 + 3Gs + 4Gc))) ZË cos? O sin? 0 + 

(15 + AG3G4Z2 + 2G3GsZ2 + AG4G3Z2 + 6G4GeZ2 + 8G4G6 Z2 + 

10G5GsZ; + 3G3G4G5GsZq + 2G3 Z5 (Gs + AG + 2G4(1 + G5G543) + 

G3(3 + 2G5G5Z3 + G4(Gs + Ziel Sall + Gi ZG (4Gu + 4G5 + 10Gg + 3G4G5G6Z; + 


2G3(3 + 2G5G6Z2 + G4(Gs + 2G6) ZE) + Ga(6 + 3GgGeZe + 





2G4(Gs + Miel Zi + Ga(G4 + 2Gs + 3G6)Z¿))) cos? 0 sin” 0 + ¡Zo(3G3 + 4G4 + 
Gs + 6G + G.G Gs Za + 2G3G1G6Z + 2G3G5G6Zó + AG4GsGeZg + 

2G5(1 + G5G6 Zo + G3Z5 (Ga + Go + G4G5G6Z3)) + 

G1(5 + Gate Sa + Mate Z + 4G5G6Z) + 2Ga Ze (G4 + 

Go + G1G5G6Zó) + GaZ2 (2G1 + Gs + AGg + 2G4G5G6Z; + 

G3(3 + 2G5G5Z3 + G4(Gs + 2G6) Z2)))) cos O sin? 0 — 

(1 + G5G5Z3 + GaZ2(G4 + Gs + GaG5GeZ2) + G1Z8(G4 + Go + G4GsGgZg + 


Go(1 + G5G6Zó + GaZe(G4 + Ge + G4G5G6Z)))) sin? 0 — 





(1/4)3 Zo(BG1 + 5G3 + AG4 + 5Gs + 10Gó + Go(6 + 4G4GsZ;)) sin? 20 


(A.10) 


The load Zr with short circuit load is: 


AZr+B B 
Zr — — 
CZ,+D D 


58 








Setting Zo = 1 and Zr = 1 and cross multiply, we obtain: 





(—j) cos? 0(—6 + G1 + 2G> + 3G5 


4G4 | DC, | 6Gg — j cot Ø) + 


(15 = 9G3 = 8G4 + 3G3Ga4 = 5G5 + 6GsG5 + 4G4G5 + 9G3G6 + HCH Ce 7 





DCH Ge + 2G2(—4 +t Gs + 2G4 + 3G5 4 4G) 





Gis 











G2 + 2G3 + 3G4 


4Gs + 5G6)) cos? 0 sinó + j(—20 + 10G3 + 8G4 = 6G3G4 + 10G; = 


6G3G5 — 4G4G5 + 3G3G4Gs + 2(10 + 2G4Gs + 


3G3(G4 + G5))G6 + 2G2(6 + 3G'5(—1 + Gç) + 2G4(—2 + Gs + 2G6) 








(esae asc dco rasa r 








2G3(-3 | Ga | 2Gs | 3G) | Ga(—4 | G3 | 2G4 | 3Gs | 4G6))) cos? 0 sin? 0 — 


(15 — 6G3 — 8G4 + 4G3G1 — 10G; + 2G3Gs + 4G4G5 — 3G3G4Gs 4 





(8G4 + 10G; + 3G3(2 ENG + 2G2(—4 + Gs + 2G4(1 + Gs(—1 + Ge)) 4 


4Gg + G3(3 + 2Gs(—1 + Ge) + G4(—2 + Gs + 2G6))) d 





Gi (Ga(4 + 3G5(—1 + Go)) + 2(—5 + 2G5 + 5G6) 





2G3(3 + 2Gs(—1 + Ge) + G4(—2 + Gs + 2G6)) + Ga(6 + 3G5(—1 + Ge) 4 








2G4(—2 + Gs + 2Gs) + Ga( —3 + Ga + 2G5 + 3G6)))) cos? 0 sin? 0 — 


j(—6 + 3G3 + 4G4 — 2G3G4 + Gs — 2G3G5 — 4G4Gs5 4 





G3G1G5 — 2G2(—1 + G3(1 + G4(—1 + G5)) + Gs) + 2G2(G3 + Gs 4 


G3G1G5)Gó + 2(3 + 2G4G5 + G3(G4 + Gelle + G1(5 + G4(—2 + Gs) — 4Gs + 


2(G4 + 2G5)G6 + 2G3(—1 + G4 — G4Gs + Go + G4G5G6) 4 





Go(—4 + Gs + 2Ga(1 + Gs(—1 + Ge)) + 4G6 + G3(3 + G4(—2 + Gs) — 2G5 4 


2(G4 + Gs)Ge)))) cos 4 sin“ å + (1+ Gs(—1+ Gg) + 





EG issu =G 











Ge + G1G5G6) + Gi(—1 + G4 — Goll + G3(1 + Ga(-1 + Gs)) + Gs) 


G4Gs(—1 ER Ge) Ze Ge + G2(G3 + G5 + G3G1G5)G6)) sin? 0 
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0 


(A.11) 


Therefore, the system of equations that we need to solve are: 


—6 + G4 + 2Gy + 3G3 + 4G4 + 5G5 + 6Gg = 0 
(15 — 9G3 — 861 + 3G3G4 — 5G5 + 6G3G5 + 4G4G; +... 
9G3G6 + 8G4G + DG5G6 + 2G2(-4 + G3 + 2G4 +... 
3G; + 4G6) + Gi(—5 + G> + 2G3 + 3G4 + 4G; + 5G6)) = 
(—20 + 10G; + 8G4 — 6G3G4 + 10G5 — 6G3G5 — 4G4Gs "t 
3G3G4G'5 + 2(10 + 2G,G; + Jed + G5))Ge +. 
2G2(6 + 3Gs(—1 bag Gg) + 2G4(—2 + Gs 4 Í 2G6) 
Gs(—3 + G4 + 2G; + 3Gg)) + G,(10 = 6G4 — 4G5 De 
3G ¿Gs + 6G4Gg + 4G5G6 + 2G3(—3 + G4 + 2G; 4 sech peus 
G3(—4 + Gs + 2G4 + 3G + 4Gg))) = 
(15 — 6G3 — 8G4 + 4G3G4 — 10G5 + 2G3G5 + dee 
3G3G4G; + (8G4 + 10G5 + _ GOGH 
2G2(— 4+ G; + 2G4(1 + G5(-1 + Ge)) 4 4Gg 4 Fesk 
G3(3 + 2G5(—1 + Ge) + G4(—2 de G's + 2Gg))) +... 
Gi(G4(4 + 3Gg(—1 + Gg)) + 2(—5 + 2G; + 5G6) +... 
2G3(3 + 2G5(—1 + Gs) + G4(—2 + Gs + 2G5)) + 
G2(6 + 3G5(—1 + Ge) + 2G4(—2 + Gs + 2G) +... 
G3(—3 + G4 + 2G5 + 3G6)))) = 0 
(—6+3G, + 4G4 — 2G3G4 + G; — 2G,G; — 4G G; +... 
G3G1G5 — 2G2(—1 + G3(1 + Ga(-1 + Gs)) + Gs) +... 
2G2(G3 + Gs + G3G4G5)Gó + 2(3 + 2G4Gs +... 
G3(G4 + Gs))G + G1(5 + Ga(—2 + Gs) — 4G5 +... 
2(G4 + 2G5)Gg + 2G3(—1 + G4 — G4Gs + Gg +... 
G4G3Gs) + Go(—4 + Gs + 2G4(1 + Gs(—1 + Gg)) +... 
AG, + G3(3 + G4(—2 + Gs) — 2G; + 2(G4 + Gs)Ge)))) = 0 
(1+ Gs(—1 + Gg) + Ga(—1 + G4 — G4Gs + Ge + G4G5G6) + 
iere d Od Gy is OE EE. 
G4Gs(—1+ Gs) + Gs + Go(G3 + Gs + G3G4G5)G6)) = 






























































(A.12) 


And the solution is: G4 = 1.9209, G> = 0.8526, Gs = 0.4303, G4 = 0.1877, Gs = 
0.0566, Gg = 0.0082. 
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dT 


E for three sections case 


A.2 Expression of 


z = [(s6:(-8 + 3G3) + G?(64 + 104G2 — 30G3 + 64G3) + G$ (64 + G2(28(1 + Gs) — 5G2(4 + G2 + G2G3))) + 4G1(—88 + G2(—64(1 + G3) + 3G2(4 + 3G2(1 + G3))))) cos + 
3(-8(-12 + G3)(2 + G> + 2G3) + 4G1(24 — G2(—16 + 3G2)(1 + G3)) + 2G2(—32(1 + Gs) + 3G2(-8 + G2(2 + 3G2 + 2G3))) + G1(—32 + 3G,(—4(1 + Gs) + Gold + G> + G2G3)))) cos 30+ 
(—8G2(—2 + G, — 2G3) — 4G1(—4 + 362) (8 + G> + G2G3) + 2G2(64(1 + Gs) — G2(—8 + 15G2(2 + G> + 2G3))) + G3(32 + Ga(4(1 + Gs) — 5G2(4 + G> + G2G3)))) cos 50+ 

Go((2 + G1)3(2 + G2)? + (32G2 + 24G2G + 12G1(4 + G3) + G3(4 + G2))Gs) cos 70— 
¡(8G2(—8 + 3G2)(1 + Gs) + 2G2(32 — G2(—52 + 15G3)(1 + G3)) + 4G1(88(1 + G3) — G,(—80 + 3G,(4 + 5G2 + 4G3))) + GH(-128(1 + Gs) + G2(—132 + 5G2(4 + 7G» + 4G3)))) sin 0+ 
3j(—8(—12 + G3)(2 + Ga + G2G3) + 4G1(G2(8 + G2(—4 + Ga — 4G3)) + 24(1 + Gs) + G3(—32(1 + Gs) + G2(—20 + G2(8 + 7G2 + 8G3))) + 2G2(—32 + G2(—8(1 + Gs) + G2(14 + G2 + G2G3)))) sin 30— 
j(—8G3(—2 + Ga + G2G3) + 4G1(—32(1 + Gs) + Gal A + 3G2(4 + 3G2 + 4G3))) + GH32(1 + Gs) + G,(—4 + G2(16 + 7G2 + 16G3))) + 2G1(—64 + Ga(—8(1 + Gs) + 3G2(14 + 3Ga(1 + G3))))) sin 50+ 

jGa((2 + G1)*(2 + G3)? + 2(24G1G2 + 2G3G + 4(4 + G3) + 3G2 (4 + G2))G3) sin 70] / 


16(cos Ø + ¿(3 + Gi + 2G» + 3G3) cos? 0 sin Ø — (3 + 2G2(1 + Gs) + Gi(2 + G> + 2G3)) cos 0 sin? 0 — j(1 + Gs + Gi(1 + G> + G2G3)) sin? ey] 


(A.13) 
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A.3 General form of reflection coefficients for 4, 5, and 6 layers case 


Below are the general form of Eq. 2.22, 2.23, and 2.24 (before substitution of Zç = 


Zo = 1 and maximally flat solutions) respectively. 
































Lla ae Ó 

(A.14) 
Ts = 

(A.15) 
To = 

(A.16) 
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A.4 MATLAB code of GA for variable length case 


clear 
for t = 1:1 % The theta limit to adjust for spacings 
close all 


hh 
ht 
1. Define a single figure of merit which will be used to evaluate all the 
solutions in the same manner. 
2. In this case, we will be using S11 as the figure of merit (the 
lower the S11 the better. (how this figure of merit will be calculated can 
be changed to suit the search requirements). 
3. S11 = (use the S11 formulas for each layer). 
4. Start the population with random genes (random number or randomly 
manually set between 0-1) 
5. Test their fitness (find their figure of merit). 
6. Choose the best 
AY 


tic; 
^^ Generating initial population 


number of population = 1000; 4 Arbitrarily set 

number of layers - 4; 

number of genes = 2*number of layers; % This equals the number of 
Jaumann layers AND the variable length 


number of generation - 1000; 


^ Initialization 
population = zeros(number of population,number of genes); 


^ Generating population 
^ Genes vary from 0-2 
population = 0 + 2*rand(number of population,number of genes); 


^^ Manually defining initial population with best survivors from 
previous search 

^4 The ones marked with ?*? are updated with latest iteration 

hh 1 layer 
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Z population(1:50,1) = 1.0157; 


^^ 2 layers 


% population(1:50,1) = 1.486353; 
% population(1:50,2) = 0.511298; 
% population(1:50,3) = 0.040800; 
% population(1:50,4) = 0.042000; 


hh 3 layers * 


% population(1:25,1) = 1.6798; 
% population(1:25,2) = 1.8442; 
% population(1:25,3) = 0.7360; 


% population(1:25,4) = 0.0830; 
% population(1:25,5) = 0.0135; 
% population(1:25,6) = 0.0284; 


population(1:25,1) = 1 
population(1:25,2) = 0.533; 


population(1:25,3) = 0.321; 
population(1:25,4) = 0.784; 
population(1:25,5) = 0.791; 
population(1:25,6) = 0.813; 


population(1,1:6) = [1.345680 0.565396 0.353227 


0.801603 0.820233 


^^ 4 layers * 

^ population(1:25,1) = 0.5636; 
% population(1:25,2) = 3.4802; 
^ population(1:25,3) = 1.1403; 
^ population(1:25,4) = 0.4061; 
H 


% population(1:25,5) = 0.9573; 


1; 


0 
% population(1:25,6) = 1.1320; 
% population(1:25,7) = 0.2257; 
% population(1:25,8) = 0.2014; 
population(1,1:8) = [1.238353 0.600010 0.330545 
8.589016 8.448836 8.666715 8.253772 


population(1,1:8) = [2.397 
0.765 0.765]; 


0.899 0.437 0.277 0.400 
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0.773810 


0.237590 
1; 


0.733 


4^ 5 layers * 


% population(1:25,1) = 0.6800; 
% population(1:25,2) = 2.3352; 
% population(1:25,3) = 1.8294; 
% population(1:25,4) = 0.7995; 
Z population(1:25,5) = 0.2149; 
^ population(1:25,6) = 0.1212; 
^ population(1:25,7) = 1.4906; 
population(1:25,8) = 0.1198; 
0.1355; 


^ 
% population(1:25,9) = 
^4 population(1:25,10) = 0.1806; 


hh 6 layers 


^ population(1:25,1) = 1.160500; 
^ population(1:25,2) = 1.405500; 
^ population(1:25,3) = 1.981900; 
^ population(1:25,4) = 2.343700; 
^ population(1:25,5) = 0.978500; 
^ population(1:25,6) = 0.350083; 


^ population(1:25,7) = 0.302400; 
^ population(1:25,8) = 0.773300; 
^ population(1:25,9) = 0.025700; 
^ population(1:25,10) = 0.428900; 
^ population(1:25,11) = 0.040500; 
^ population(1:25,12) = 0.042700; 


^^ Loading the entire table 

At 

fID = fopen(’C:\Users\Wiwi Samsul\Documents\MATLAB\GA_results.txt’); 
B = textscan(fID,’%s’,’ HeaderLines’ ,1,’delimiter’,’\t’); 
fclose(fID); 


TableSize = length(B{1}) ; 

NumberOfLayers = zeros(1,round(TableSize/360)) ; 
MaxNumberOfLayers = 6; 

NumberOfEntriesPerLayer = 20; 

NumberOfGenes = 6; 


% test = cell2mat(Bi1) (10) 
for n = 1:round(TableSize/360) 

NumberOfLayers(1,n) = str2double(cell2mat (B{1}(1 + (n-1)*360))); 
end 
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GValues = zeros(NumberÜfEntriesPerLayer,NumberÜfGenes,length( 
NumberOfLayers) ) ; 

XValues = zeros (Number0fEntriesPerLayer , NumberOfGenes, length ( 
NumberOfLayers) ) ; 

RelativeBandwidthValues = zeros(Number0fEntriesPerLayer,1,length( 
NumberOfLayers) ) ; 

StartBandwidthValues = zeros (Number0fEntriesPerLayer ,1,length( 
NumberOfLayers) ) ; 

EndBandwidthValues = zeros(NumberÜfEntriesPerLayer,1,length( 





NumberOfLayers) ) ; 


for m = 1:NumberOfEntriesPerLayer 
for n = 1:NumberDfGenes 
for p = 1:length(NumberOfLayers) 

GNaNCheck = isnan(str2double(cel12mat (B{1}(18*(m-1) + (2+(n 
-1)) + 360*(p-1))))); 

if (”GNaNCheck) 
GValues(m,n,p) = str2double(cell2mat (B{1}(18*(m-1) + 

(2+(n-1)) + 360*(p-1)))); 

end 

XNaNCheck = isnan(str2double(cel12mat (B{1}(18*(m-1) + (8+(n 
-1)) + 360*(p-1))))); 

if (^XNaNCheck) 
XValues(m,n,p) = str2double(cell2mat (B{1}(18*(m-1) + 

(8+(n-1)) + 360*(p-1)))); 

end 

RelativeBandwidthValues(m,1,p) = str2double(cell2mat (B 
{1}(18*(m-1) + (14) + 360*(p-1)))); 

StartBandwidthValues(m,1,p) = str2double(cell2mat (B{1}(18*(m 
-1) + (15) + 360*(p-1)))); 

EndBandwidthValues(m,1,p) = str2double(cell2mat (B{1}(18*(m 
-1) + (16) + 360*(p-1)))); 

end 
end 
end 


% Automatically populate with the newest result taken from the table 
for m = 1:number_of_layers 
population(1:25,m) = GValues(t,m,number_of_layers) ; 
population(1:25,(m + number_of_layers)) = XValues(t,m, 
number_of_layers) ; 
end 


AY 
hh 
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fitness_collection = zeros(1,number_of_population); % Initialization 
all history fitness collection = zeros(number of generation, 
number of population); % For collecting all history of fitness over 
generations 
all history generation collection = zeros(number of generation, 
number of population,number of genes); ^ For collecting all history 
of population 
best fitness gene = zeros(1,number of genes); % To store the best 
gene 
best fitness record = 0; % To store the best recorded fitness 
generation number record = 0; % To record at which generation the value 


was found 


n = 0; 
for n = 1:number_of_generation 
if (mod(n,100)==0) 
fprintf(*%d th gen finished. \n’,n); 
end 
% while ((max(fitness_collection) < 1.4018) && (n < 
number of generation)) 
% n=n+ 1; 
for m = 1:number_of_population 


fitness_collection(1,m) = SllFitness_with_var_L_var_theta( 
population(m,:),t); 


end 


all_history_fitness_collection(n,:) = fitness_collection(1,:); 
all_history_generation_collection(n,:,:) = population(:,:); 


[fitness_sorted, original_index_in_sorted_order] = sort( 
fitness collection,'descend?); 


^4 Take the better half/quarter of population based on its fitness 
survivors = zeros(ceil(0.5*number of population) ,number of genes); 


^ Record the best fitness 
if(best fitness record < fitness sorted(1)) 
best fitness record - fitness sorted(1); 
best fitness gene - population(original index in sorted order 
(1),:); 
generation_number_record = n; 
end 
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for m = 1:length(survivors) 
survivors(m,:) = population(original_index_in_sorted_order (m) 


DE 


end 


new gen = zeros(number of population,number of genes); 


ht 
^ COMMENT BLOCK SPECIAL CASE FOR 1 GENE POPULATIONS 


% Crossover process (by averaging) 

^ Select two surviving parents randomly and average their values 
% This only works for population with 1 gene 

for m = 1:number_of_population 


^4 Create a random number between 1-length(surviving gen) 


% To randomly take 2 parents 
parenti index = round(1 + (length(survivors) - 1).*rand(1,1)); 
parent2_index = round(1 + (length(survivors) - 1).*rand(1,1)); 


new gen(m) = 0.5*(survivors(parent1_index)+survivors (parent2_index) 
); 
end 


hy 


^4 Cross over genes 
for m = 1:number_of_population 


^4 Randomly take 2 parents 
parenti index = round(1 + (length(survivors) - 1).*rand(1,1)); 
parent2 index = round(1 + (length(survivors) - 1).*rand(1,1)); 


^4 Send for cross over 
new gen(m,:) = crossover(survivors(parenti index,:),survivors( 


parent2 index,:)); 
end 


^ Mutation 

mutation magnitude - 0.5; 

mutation, chance = 0.2; 

population = abs(mutation(new gen,mutation magnitude, 
mutation chance)); 


end 
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x_axis = 1:number_of_population; 


figl = figure(1); 
for m = 1:number_of_generation 
hold on 
plot(x_axis,all_history_fitness_collection(m,:)); 
hold on 
end 
grid on 
% axis([1 number of population floor(min(min( 
all history fitness collection))) ceil (max (max ( 
all history fitness collection)))]) 
axis([1 number of population O ceil(max(max( 
all history fitness, collection)))]) 
xlabel(?Individual number / Solutions?) 
ylabel(’Fitness / Relative bandwidth at -10 dB?) 
title(C Evolution of solutions?) 
^4 legend(?Gen 1’,’Gen 2’,’Gen 3’,’Gen 4’,’Gen 5’,’Gen 6’,’Gen 7’,’Gen 
8’,’Gen 9’,’Gen 10’) 
image filename = sprintf(?Genetic algorithm result % 
d layers with var L population %d theta.png?,number of layers,t); 


saveas(gcf,image filename); 


% Line markers 


lineMarkers = [’-k ’;’:k ?;?--k ?;?-.k ?;?-sk ?;”:sk ’;’--sk’;’-.sk 
?;?-+k Pe kk YEAR KIS = Ok > sok >;?--ok’;’?-.0k’;’-xk Bae 
xk ?;?--xk?;?-.xk?;?-dk ?;”:dk ’?;’--dk’;’-.dk’]; 


lineMarkers = cellstr(lineMarkers); 


fig2 = figure(2); 
hold on 
for m = 1:number_of_genes 
plot(all history generation collection(number of generation,:,m), 
char (lineMarkers (m))); 
end 
grid on 
^4 axis([1 number of population floor(min(min(min( 
all history generation collection(number of generation,:,:))))) 
ceil(max(max(max((all history generation collection( 
number of generation,:,:))))))D; 
axis([1 number of population O ceil(max(max(max(( 
all history generation collection(number of generation,:,:))))))]D; 
xlabel(?Individual number / Solutions?) 
ylabel(’Gene value?) 
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title(?Last generation?) 

legend(’1st gene’,’2nd gene’,’3rd gene’,’4th gene’,’5th gene’,’6th gene 
?,?7th gene’,’8th gene’,’9th gene’,’10th gene’,’11th gene’,’12th 
gene?) 


image filename = sprintf(”Genetic algorithm result % 
d layers with var L genes Ad theta.png!,number of layers,t); 
% saveas(gcf,image filename); 


maximum fitness - max(max(all history fitness collection)); 
display text = sprintf(’Maximum fitness found: %f’,maximum fitness); 
disp(display text) 


display text = sprintf(?Best fitness record: %f’ ‚best fitness record); 
disp(display text) 
display text = sprintf(?Best fitness gene: \n’); 
disp(display text) 
display text = ??; 
add gene = sprintf(’%f’ ,best fitness gene(1)); 
for n = 2:number_of_genes 
display text = [display text add gene ’\t’]; 
add gene = sprintf (*%f? ,best fitness gene(n)); 
end 
display text = [display text add gene ’\t’]; 
display text = sprintf(display text); 
disp(display text) 


double check = S1i1Fitness with var L var theta(best fitness gene,t); 

display text = sprintf(’Double checking the best fitness gene: %f”, 
double check); 

disp(display text) 


display text = sprintf(’Solution found at generation tt: %d’, 
generation number record); 
disp(display text) 


% Calling in the function to graph the result with best gene: 

[ fitness,theta in Hz,gamma | = Si1Fitness with var L var theta( 
best fitness gene,t ); 

SiiFitnessGraph with var L( theta in Hz, gamma ) 

image filename = sprintf('Genetic algorithm result % 
d. layers with var L genes $11, plot "^d theta.png?,number of layers,t 
); 


% saveas(gcf,image filename); 
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toc; 


total_search_time = toc; 

effective_search_time = toc*generation_number_record/ 
number_of_generation; 

display text = sprintf(’Estimated search time: JE", 
effective search, time); 

disp(display text) 


end 
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A.5 Supplemental MATLAB code of for GA 


function [ fitness,theta in Hz,gamma ] = S$11Fitness_with_var_L_var_theta( 
individual,t ) 

% fitness = $11Fitness(individual, numOfGenes) 

^ SiiFitness calculates the relative bandwidth of S11 at dB level from 
individual solutions 

^4 individual must be in normalized G 

Á Converting between L and x : 

^4 L \lambda = 0.25 * x Ntheta 


ZS = 1; 
ZO = 1; 
fO = 2.50e2; 


theta = [0:pi/8000:t*pi]; 
theta in Hz = 2*theta*f0/pi; 


numOfGenes = round(length(individual)/2); 
dB_level = -10; 


switch num0fGenes 
case 1 
G1 = individual (1); 
x1 = individual (2) ; 
gamma = mag2db(abs(-((ZS.*cos(x1.*theta) + 1i.*Z0.*(-1 + G1.*ZS).* 
sin(x1.*theta))./(ZS.*cos(x1.*theta) + 1i.*Z0.*(1 + G1.*ZS).*sin 
(x1.*theta))))); 
case 2 
G1 = individual (1); 
G2 = individual (2); 
x1 = individual (3); 
x2 = individual (4) ; 
gamma = mag2db(abs((cos(x1.*(theta)).*(ZS.*cos(x2.*(theta)) + 1i.* 
ZO.*(-1 + G2.*ZS) .*sin(x2.*(theta))) + sin(x1.*(theta)).*(1i.*Z0 
.*(-1 + G1.*ZS + G2.*Z8).*cos(x2.*(theta)) - (ZS + G1.*Z072.*(-1 
+ G2.*ZS)).*sin(x2.*(theta))))./ ((-cos(x1.*(theta))).*(ZS.*cos 
(x2.*(theta)) + 11.*Z0.*(1 + G2.*Z8).*sin(x2.*(theta))) + sin(x1 
.*(theta)) .*((-1i).*Z0.*(1 + G1.*ZS + G2.*ZS) .*cos(x2.* (theta) ) 
+ (ZS + G1.*Z072.*(1 + G2.*ZS)) .*sin(x2.*(theta)))))); 
case 3 
G1 individual(1); 
G2 = individual(2); 
G3 = individual(3); 
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x1 
x2 
x3 


individual (4); 
individual(5); 
= individual(6); 


gamma = mag2db(abs((cos(x1.*(theta)).*(cos(x2.*(theta)).*((-11).*ZS 


case 4 
G1 
G2 
G3 


G4 = 


x1 
x2 
x3 
x4 


.*cos(x3.*(theta)) + Z0.*(-1 + G3.*ZS).*sin(x3.*(theta))) + Sint 
x2.*(theta)).*(Z0.*(-1 + G2.*ZS + G3.*ZS).*cos(x3.*(theta)) + 1i 
.* (ZS + G2.*Z0°2.*(-1 + G3.*ZS)).*sin(x3.*(theta)))) + sin(x1.x*( 
theta)).*(cos(x2.*(theta)).*(ZO.*(-1 + G1.*ZS + G2.*ZS + G3.*ZS) 
.*cos(x3.*(theta)) + 1i.*(ZS + G1.*Z0°2.*(-1 + G3.*ZS) + G2.*Z0 
72.*(-1 + G3.*ZS)).*sin(x3.*(theta))) + sin(x2.*(theta)).*(1i.*( 
ZS + Gi.*Z0°2.*(-1 + G2.x*ZS + G3.*ZS)).*cos(x3.*(theta)) - ZO 
.*(-1 + G3.4ZS + G1.*(ZS + G2.*Z0^2.*(-1 + G3.*ZS))) .*sin(x3.*( 
theta)))))./ (1i.*cos(x1.* (theta) ).*(cos(x2.* (theta) ) .*(ZS.*cos ( 
x3.*(theta)) + 1i.*Z0.*(1 + G3.*ZS) .*sin(x3.*(theta))) - sin(x2 
.*(theta)).*((-1i).*Z0.*(1 + G2.*ZS + G3.*Z8) .*cos(x3.* (theta) ) 
+ (ZS + G2.*Z2072.*(1 + G3.*ZS)).*sin(x3.*(theta)))) + sin(x1.*( 
theta)).*((-1i).*cos(x2.*(theta)).*((-1i).*Z0.*(1 + G1.*ZS + G2 
.*Z8 + G3.*ZS).*cos(x3.*(theta)) + (ZS + G1.*Z2072.*(1 + G3.*ZS) 
+ G2.*Z0°2.*(1 + G3.*ZS)).*sin(x3.*(theta))) + sin(x2.*(theta)) 
.*((-1i).*(ZS + G1.*Z072.*(1 + G2.*ZS + G3.*ZS)).*cos(x3.*(theta 
)) + ZO.*(1 + G3.4ZS + G1.*(ZS + G2.*Z072.*(1 + G3.*ZS))).*sin( 
x3.*(theta))))))); 


= individual(1); 
= individual(2); 
= individual(3); 
individual(4); 
= individual (5); 
= individual (6) ; 
= individual (7); 
= individual (8) ; 


gamma = mag2db(abs((-ZS + (Z0.*(cos(x2.* (theta) ) .*(cos(x3.* (theta) ) 


.*((-1i).*cos(x4.*(theta)).*sin(x1.*(theta)) + ((-1i).*cos(x1.*( 
theta)) + (G1 + G2 + G3) .*Z0.*sin(x1.*(theta))).*sin(x4.*(theta) 
)) + sin(x3.*(theta)).*(cos(x1.*(theta)) .*((-1i).*cos(x4.*(theta 
)) + G3.*Z0.*sin(x4.*(theta))) + sin(x1.*(theta)).*((G1 + G2).x* 
ZO.*cos(x4.*(theta)) + 1i.*(1 + G1.*G3.*Z0°2 + G2.*G3.*Z072) .* 
sin(x4.*(theta))))) + sin(x2.*(theta)).*(cos(x1.*(theta)).*(cos( 
x3.*(theta)).*((-1i).*cos(x4.*(theta)) + (G2 + G3) .*Z0.*sin(x4 
.*(theta))) + 1i.*sin(x3.*(theta)).*((-1i) .*G2.*Z0.*cos (x4. *( 
theta)) + (1 + G2.*G3.*Z0^2).*sin(x4.*(theta)))) + sin(x1.x*( 
theta)).*((-sin(x3.*(theta))).*((-1i).*(1 + G1.*G2.*Z072) .*cos( 
x4.*(theta)) + ZO.*(G1 + G3 + G1.*G2.*G3.*Z072) .*sin(x4.*(theta) 
)) + 
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1i.*cos(x3.*(theta)).*((-1i).*G1.*Z0.*cos(x4.*(theta)) + (1 + G1.*(G2 + 
G3) .*Z0°2).*sin(x4.*(theta)))))))./ (cos(x1.*(theta)).*(cos(x2.*(theta 
)) .*C(-cos(x3.* (theta) )).*(cos(x4.*(theta)) + 11.*G4.*Z0.*sin(x4. zt 
theta))) + sin(x3.*(theta)).*((-1i).*(G3 + G4) .*Z0.*cos(x4.*(theta)) + 

(1 + G3.*G4.*Z2072) .*sin(x4.*(theta)))) + sin(x2.*(theta)).*(cos (x3. zt 
theta) ).*((-1i).*(G2 + G3 + G4).*Z0.*cos(x4.*(theta)) + (1 + G2.*G4.* 
ZO°2 + G3.*G4.*Z0^2) .*sin(x4.*(theta))) + sin(x3.*(theta)).*((1 + G2 
.*(G3 + G4) .*Z0°2).*cos(x4.*(theta)) + 11.*Z0.*(G2 + G4 + G2.*G3.*G4.* 
20^2).*sin(x4.*(theta))))) + sin(x1.*(theta)).*(cos(x2.*(theta)).*(cos 
(x3.*(theta)).*((-1i).*(G1 + G2 + G3 + G4).*Z0.*cos(x4.*(theta)) + (1 
+ G1.*G4.*Z0^2 + G2.*G4.*Z0°2 + G3.*G4.*Z072) .*sin(x4.*(theta))) + sin 
(x3.*(theta)).*((1 + G1.*(G3 + G4).*Z0°2 + G2.*(G3 + G4) .*Z072) .*cos( 
x4.*(theta)) + 1i.*Z0.*(G1 + G2 + G4 + 

G1.*G3.*G4.*Z072 + G2.*G3.*G4.+*Z0°2) .*sin(x4.*(theta)))) + sin(x2.*(theta 
)).*(cos(x3.*(theta)).*((1 + G1.*(G2 + G3 + G4).*Z0°2).*cos(x4.* (theta 
)) + 1i.*Z0.*(G1 + G4 + G1.*G2.*G4.*Z2072 + G1.*G3.*G4.*Z072) .*sin(x4 
.*(theta))) - sin(x3.*(theta)).*((-1i).*Z0.*(G3 + G4 + G1.*(1 + G2.*( 
G3 + G4) .*Z072)) .*cos(x4.*(theta)) + (1 + G3.*G4.*Z0"2 + G1.*Z072.*(G2 

+ G4 + G2.*G3.*G4.*Z072)) .*sin(x4.*(theta)))))))./ (ZS + (ZO.*(cos(x2 
.*(theta)).*(cos(x3.*(theta)).*((-1i).*cos(x4.*(theta)).*sin(x1.*( 
theta)) + ((-1i).*cos(x1.*(theta)) + (G1 + G2 + G3) .*Z0.*sin(x1.*( 
theta))).*sin(x4.*(theta))) + sin(x3.*(theta)).*(cos(x1.*(theta)) 
.*((-1i).*cos(x4.*(theta)) + G3.*Z0.*sin(x4.*(theta))) + sin(x1.x*( 
theta)).*((G1 + G2).*Z0.*cos(x4.*(theta)) + 1i.*(1 + G1.*G3.*Z072 + G2 
.*G3.*Z0°2).*sin(x4.*(theta))))) + 

sin(x2.*(theta)) .*(cos(x1.* (theta) ).*(cos(x3.*(theta)) .*((-1i) .*cos (x4. zt 
theta)) + (G2 + G3).*Z0.*sin(x4.*(theta))) + 1i.*sin(x3.*(theta)) 
.*((-1i).*G2.*Z0.*cos(x4.*(theta)) + (1 + G2.*G3.*Z0^2) .*sin(x4.*( 
theta)))) + sin(x1.*(theta)).*((-sin(x3.*(theta))).*((-1i).*(1 + G1.x* 
G2.*Z072) .*cos(x4.*(theta)) + ZO.*(G1 + G3 + G1.*G2.*G3.*Z20^2) .*sin(x4 
.*(theta))) + 1i.*cos(x3.*(theta)).*((-1i).*G1.*Z0.*cos(x4.*(theta)) + 

(1 + G1.*(G2 + G3).*Z0°2).*sin(x4.*(theta)))))))./ (cos(x1.*(theta) ) 
.*(cos(x2.* (theta) ) .*((-cos(x3.*(theta))).*(cos(x4.*(theta)) + 1i.*G4 
.*Z0.*sin(x4.*(theta))) + sin(x3.*(theta)).*((-1i).*(G3 + G4).*Z0.*cos 
(x4.*(theta)) + (1 + G3.*G4.*Z0^2) .*sin(x4.*(theta)))) + sin(x2.*( 
theta)).*(cos(x3.*(theta)).*((-1i).*(G2 + G3 + G4).*ZO.*cos(x4.*(theta 
)) + (1 + G2.*G4.*Z0^2 + G3.*G4.*Z072) .*sin(x4.*(theta))) + sin(x3.*( 
theta)).*((1 + G2.*(G3 + G4).*Z0°2).*cos(x4.*(theta)) + 

1i.*Z0.*(G2 + G4 + G2.*G3.*G4.*Z072) .*sin(x4.*(theta))))) + sin(x1.*( 

theta)).*(cos(x2.*(theta)).*(cos(x3.*(theta)).*((-1i).*(G1 + G2 + G3 + 

G4) .*Z0.*cos(x4.*(theta)) + (1 + G1.*G4.*Z0^2 + G2.*G4.*Z0^2 + G3.*G4 
.*Z072) .*sin(x4.*(theta))) + sin(x3.*(theta)).*((1 + G1.*(G3 + G4).*ZO 
72 + G2.*(G3 + G4) .*Z0^2) .*cos(x4.*(theta)) + 11.*Z0.*(G1 + G2 + G4 + 
G1.*G3.*G4.*Z072 + G2.*G3.*G4.*Z072) .*sin(x4.*(theta)))) + sin(x2.*( 

theta)).*(cos(x3.*(theta)).*((1 + G1.*(G2 + G3 + G4) .*Z0^2) .*cos (x4.*( 
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theta)) + 11.*Z0.*(G1 + G4 + G1.*G2.*G4.*Z0^2 + G1.*G3.*G4.*20^2).*sin 
(x4.*(theta))) - sin(x3.*(theta)).*((-1i).*Z0.*(G3 + G4 + G1.*(1 + G2 
.*(G3 + G4) .*Z0°2)).*cos(x4.*(theta)) + (1 + G3.*G4.*Z0^2 + G1.*Z0 
72.*(G2 + G4 + G2.*G3.*G4.*Z0°2)) .*sin(x4.*(theta))))))))); 


case 5 


G1 
G2 
G3 


x1 
x2 
x3 
x4 
x5 


= individual (1); 
= individual (2); 
= individual (3) ; 
G4 = 
G5 = 
= individual (6) ; 
= individual (7); 
= individual (8) ; 
= individual (9); 
= individual (10); 


individual (4); 
individual (5) ; 


gamma = mag2db(abs((-ZS + (1i.*Z0.*(cos(x1.*theta).*(cos(x2.*theta) 


.* (cos (x4. *theta) .*(cos(x5.*theta) .*sin(x3.*theta) + (cos (x3.* 
theta) + 1i.*(G3 + G4) .*Z0.*sin(x3.*theta)).*sin(x5.*theta)) + 
sin(x4.*theta) .*(cos(x3.*theta) .*(cos(x5.*theta) + 11.*G4.*Z0.* 
sin(x5.*theta)) - sin(x3.*theta).*((-1i).*G3.*Z0.*cos(x5.*theta) 
+ (1 + G3.*G4.*Z0°2).*sin(x5.*theta)))) + sin(x2.*theta).*(cos( 
x3.*theta).*(cos(x4.*theta).*(cos(xb.*theta) + 1i.*(G2 + G3 + G4 
).*Z0.*sin(x5.*theta)) - sin(x4.*theta).*((-1i).*(G2 + G3).*Z0.* 
cos(x5.*theta) + (1 + G2.*G4.*Z0^2 + G3.*G4.*Z072) .*sin(x5.* 
theta))) + sin(x3.*theta).*((-sin(x4.*theta)).*((1 + G2.*G3.*Z0 
72) .*cos(x5.*theta) + 1i.*Z0.*(G2 + G4 + G2.*G3.*G4.*20^2) .*sin( 
xb.*theta)) - cos(x4.*theta) .*((-1i) .*G2.*Z0.*cos(x5.*theta) + 
(1 + G2.*(G3 + G4).*Z0°2).*sin(x5.*theta))))) + sin(x1.*theta) 
.*((-sin(x2.*theta)).*(cos(x4.*theta).*(cos(x5.*theta).*sin(x3.* 
theta) + (cos(x3.*theta) + 


1i.*(G3 + G4) .*Z0.*sin(x3.*theta)).*sin(x5.*theta)) + sin(x4.*theta) .*( 
cos(x3.*theta).*(cos(x5.*theta) + 1i.*G4.*Z0.*sin(x5.*theta)) - sin(x3 
.*theta) .*((-1i).*G3.*Z0.*cos(x5.*theta) + (1 + G3.*G4.*Z0^2) .*sin(x5 
.*theta)))) + cos(x2.*theta).*(cos(x3.*theta).*(cos(x4.*theta).*(cos( 
xb.*theta) + 1i.*(G2 + G3 + G4) .*Z0.*sin(x5.*theta)) - sin(x4.*theta) 
*((-1i) .*(G2 + G3).*Z0.*cos(x5.*theta) + (1 + G2.*xG4.*Z0"2 + G3.*G4.* 
20^2).*sin(xb.*theta))) + sin(x3.*theta).*((-sin(x4.*theta)).*((1 + G2 
.*G3.*Z072).*cos(x5.*theta) + 1i.*Z0.*(G2 + G4 + G2.*G3.*G4.*Z072) .* 
sin(x5.*theta)) - cos(x4.*theta).*((-1i).*G2.*Z0.*cos(x5.*theta) + (1 
+ G2.*(G3 + G4).*Z0°2).*sin(x5.*theta)))) + 1i.*G1.*Z0.*(cos(x2.*theta 
).* (cos (x4.*theta) .*(cos(x5.*theta).*sin(x3.*theta) + (cos(x3.*theta) 
+ 11.*(G3 + G4).*Z0.*sin(x3.*theta)).*sin(x5.*theta)) + sin(x4.*theta) 
.*(cos(x3.*theta).*(cos(xb.*theta) + 11.*G4.*Z0.*sin(x5.*theta)) - sin 
(x3.*theta) .*((-1i) .*G3.*Z0.*cos(x5.*theta) + (1 + 
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G3.*G4.*Z0°2).*sin(x5.*theta)))) + sin(x2.*theta) .*(cos(x3.*theta) .*(cos( 
x4.*theta) .*(cos(x5.*theta) + 11.*(G2 + G3 + G4).*Z0.*sin(x5.*theta)) 
- sin(x4.*theta) .*((-1i).*(G2 + G3).*Z0.*cos(x5.*theta) + (1 + G2.*G4 
.*Z0^2 + G3.*G4.*Z072) .*sin(x5.*theta))) + sin(x3.*theta).*((-sin(x4.* 
theta)).*((1 + G2.*G3.*Z072) .*cos(x5.*theta) + 1i.*Z0.*(G2 + G4 + G2.* 
G3.*G4.*Z0^2) .*sin(x5.*theta)) - cos(x4.*theta) .*((-1i) .*G2.+*Z0.*cos( 
x5.*theta) + (1 + G2.*(G3 + G4) .*Z072) .*sin(x5.*theta))))))))./(cos(x1 
.*theta) .*(cos(x2.*theta) .*(cos(x3.*theta) .*(cos(x4.*theta) .*(cos(x5.* 
theta) + 11.*G5.*Z0.*sin(x5.*theta)) - sin(x4.*theta).*((-1i).*(G4 + 
GB) .*Z0.*cos(x5.*theta) + (1 + G4.*G5.*Z0°2).*sin(x5.*theta))) + 1i.* 
sin(x3.*theta) .*(cos(x4.*theta) .*((G3 + G4 + G5).*Z0.*cos(x5.*theta) + 

1i.*(1 + G3.*G5.*Z0^2 + G4.*G5.*Z0^2) .*sin(xb.*theta)) + 1i.*sin(x4.* 
theta) .*((1 + G3.*(G4 + G5).*Z072).*cos(x5.*theta) + 1i.*Z0.*(G3 + G5 
+ 

G3.*G4.*G5.*Z0°2) .*sin(x5.*theta)))) + 1i.*sin(x2.*theta) .*(cos(x3.*theta 
).*(cos(x4.*theta).*((G2 + G3 + G4 + G5).*Z0.*cos(x5.*theta) + 1i.*(1 
+ G2.*G5.*Z0°2 + G3.*G5.*Z0°2 + G4.*G5.*Z072) .*sin(x5.*theta)) + 1i.* 
sin(x4.*theta).*((1 + G2.*(G4 + G5).*Z0°2 + G3.*(G4 + G5) .*Z072) .*cos( 
x5.*theta) + 11.*Z0.*(G2 + G3 + G5 + G2.*G4.*G5.*2072 + G3.*G4.*G5.*Z0 
72) .*sin(x5.*theta))) + 1i.*sin(x3.*theta).*(cos(x4.*theta).*((1 + G2 
.*(G3 + G4 + G5).*Z0°2).*cos(x5.*theta) + 1i.*Z0.*(G2 + Gb + G2.*G3.x* 
G5.*Z0^2 + G2.*G4.*G5.*Z072) .*sin(x5.*theta)) - sin(x4.*theta).*((-1i) 
.*Z0.*(G4 + G5 + G2.*(1 + G3.*(G4 + G5).*Z072)).*cos(x5.*theta) + (1 + 
G4.*G5.*Z0°2 + G2.*Z0°2.*(G3 + Gb + G3.*G4.*G5.*Z0°2)).*sin(x5.*theta 
))))) + 1i.*sin(x1.*theta) .*(cos(x2.*theta) .*(cos(x3.*theta) .*(cos(x4 
.*theta).*((G1 + G2 + G3 + G4 + 

G5) .*Z0.*cos(x5.*theta) + 11.*(1 + G1.*G5.*Z0°2 + G2.*G5.*Z0°2 + G3.*G5.* 
Z0^2 + G4.*G5.*Z0°2).*sin(x5.*theta)) + 1i.*sin(x4.*theta).*((1 + G3.* 
G4.*Z0^2 + G3.*G5.*Z0°2 + G1.*(G4 + G5).*Z0°2 + G2.*(G4 + G5).*Z072).x* 
cos(x5.*theta) + 1i.*Z0.*(G1 + G2 + G3 + G5 + G1.*G4.*G5.*Z20^2 + G2.* 
G4.*G5.*Z0°2 + G3.*G4.*G5.*Z072) .*sin(x5.*theta))) + 1i.*sin(x3.*theta 
).*(cos(x4.*theta).*((1 + G1.*(G3 + G4 + G5) .*Z072 + G2.*(G3 + G4 + G5 
).*Z072) .*cos(x5.*theta) + 11.*Z0.*(G1 + G2 + G5 + G1.*G3.*G5.*Z072 + 
G2.*G3.*G5.*Z0^2 + G1.*G4.*G5.*Z0°2 + G2.*G4.*G5.*Z072).*sin(xb.*theta 
)) - sin(x4.*theta).*((-1i).*Z0.*(G2 + G4 + G5 + G2.*G3.*G4.*Z072 + G2 
.*G3.*G5.*Z0^2 + G1.*(1 + G3.*(G4 + G5) .*Z0°2)).*cos(x5.*theta) + (1 + 
G4.*G5.*Z0°2 + G1.*Z0°2.*(G3 + Gb + G3.*G4.+*G5.*Z0°2) + G2.*Z0°2.*(G3 
+ G5 + G3.*G4.*G5.*Z0°2)).*sin(x5.*theta)))) - sin(x2.*theta) .*(cos( 
x3.*theta).*((-1i).*cos(x4.*theta).*((1 + G1.*(G2 + G3 + G4 + 

GB) .*Z072) .*cos(x5.*theta) + 1i.*Z0.*(G5 + G1.*(1 + G2.*G5.*Z0°2 + G3.*G5 
.*Z072 + G4.*G5.*Z072)) .*sin(x5.*theta)) + sin(x4.*theta).*(Z0.*(G4 + 
GB + G1.*(1 + G2.*(G4 + G5).*Z0"2 + G3.*(G4 + G5) .*Z072)) .*cos(xb.* 
theta) + 1i.*(1 + G4.*G5.*Z0"2 + G1.*Z0°2.*(G2 + G3 + G5 + G2.*G4.*G5 
.*Z0^2 + G3.*G4.*G5.*Z072)).*sin(x5.*theta))) + sin(x3.*theta) .*(cos( 
x4.*theta) .*(Z0.*(G3 + G4 + G5 + Gi.*(1 + G2.*(G3 + G4 + G5).*Z0^2)) .* 
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cos(x5.*theta) + 1i.*(1 + G3.*G5.*Z072 + G4.*G5.*Z0"2 + G1.*Z072.*(G2 
+ G5 + G2.*G3.*G5.*Z0°2 + G2.*G4.*G5.*Z072)) .*sin(x5.*theta)) + 1i.* 
sin(x4.*theta).*((1 + G3.*(G4 + G5).*Z0°2 + G1.*Z0°2.*(G4 + G5 + G2 
(1 + G3.*(G4 + G5).*Z072))).*cos(x5.*theta) + 1i.*Z0.*(G3 + G5 + G3 
.*G4.*G5.*Z072 + G1.*(1 + G4.*G5.*Z0°2 + G2.*Z0^2.*(G3 + G5 + G3.*G4.* 
G5.*Z0°2))).*sin(x5.*theta)))))))./(ZS + (1i.*Z0.*(cos(x1.*theta).*( 
cos(x2.*theta) .*(cos(x4.*theta) .*(cos(x5.*theta).*sin(x3.*theta) + 

(cos(x3.*theta) + 11.*(G3 + G4).*Z0.*sin(x3.*theta)).*sin(x5.*theta)) + 
sin(x4.*theta) .*(cos(x3.*theta) .*(cos(x5.*theta) + 11.*G4.*Z0.*sin(x5 
.*theta)) - sin(x3.*theta).*((-1i).*G3.*Z0.*cos(x5.*theta) + (1 + G3.* 
G4.*Z2072) .*sin(x5.*theta)))) + sin(x2.*theta) .*(cos(x3.*theta) .*(cos( 
x4.*theta).*(cos(x5.*theta) + 11.*(G2 + G3 + G4).*Z0.*sin(xb.*theta)) 
- sin(x4.*theta).*((-1i).*(G2 + G3).*Z0.*cos(x5.*theta) + (1 + G2.*G4 
.*Z0^2 + G3.*G4.*Z072) .*sin(x5.*theta))) + sin(x3.*theta).*((-sin(x4.* 
theta)).*((1 + G2.*G3.*Z072) .*cos(x5.*theta) + 1i.*Z0.*(G2 + G4 + G2.x* 
G3.*G4.*Z0^2) .*sin(x5.*theta)) - cos(x4.*theta) .*((-1i) .*G2.+*Z0.*cos( 
x5.*theta) + (1 + G2.*(G3 + G4).*Z0°2).*sin(x5.*theta))))) + sin(x1.x* 
theta) .*((-sin(x2.*theta)).*(cos(x4.*theta) .*(cos(x5.*theta) .*sin(x3.* 
theta) + (cos(x3.*theta) + 1i.*(G3 + G4).*Z0.*sin(x3.*theta)).*sin(x5 
.*theta)) + sin(x4.*theta) .*(cos(x3.*theta) .*(cos(x5.*theta) + 

1i.*G4.*Z0.*sin(x5.*theta)) - sin(x3.*theta).*((-1i).*G3.*Z0.*cos(x5.* 
theta) + (1 + G3.*G4.*Z072) .*sin(x5.*theta)))) + cos(x2.*theta).*(cos( 
x3.*theta).*(cos(x4.*theta).*(cos(x5.*theta) + 1i.*(G2 + G3 + G4) .*Z0 
.*sin(x5.*theta)) - sin(x4.*theta).*((-1i).*(G2 + G3) .*Z0.*cos(x5.* 
theta) + (1 + G2.*G4.*Z0°2 + G3.*G4.*Z20^2).*sin(xb.*theta))) + sin(x3 
.*theta).*((-sin(x4.*theta)).*((1 + G2.*G3.*Z0°2).*cos(x5.*theta) + 1i 
.*Z0.*(G2 + G4 + G2.*G3.*G4.*Z0°2).*sin(x5.*theta)) - cos(x4.*theta) 
.*((-1i).*G2.*Z0.*cos(xb.*theta) + (1 + G2.*(G3 + G4).*Z0^2) .*sin(xb.* 
theta)))) + 11.*G1.*Z0.*(cos(x2.*theta) .*(cos(x4.*theta) .*(cos (x5. * 
theta) .*sin(x3.*theta) + (cos(x3.*theta) + 1i.*(G3 + G4).*Z0.*sin(x3.x* 
theta)).*sin(x5.*theta)) + sin(x4.*theta) .*(cos(x3.*theta) .*(cos(x5.* 
theta) + 1i.*G4.*Z0.*sin(x5.*theta)) - sin(x3.*theta).*((-1i).*G3.*Z0 
.*cos(x5.*theta) + (1 + G3.*G4.*20^2).*sin(xb.*theta)))) + sin(x2.x* 
theta) .*(cos(x3.*theta) .*(cos(x4.*theta) .*(cos(x5.*theta) + 

1i.*(G2 + G3 + G4).*Z0.*sin(x5.*theta)) - sin(x4.*theta).*((-1i).*(G2 + 
G3) .*Z0.*cos(x5.*theta) + (1 + G2.*G4.*Z0°2 + G3.*G4.*Z0°2) .*sin(xb.* 
theta))) + sin(x3.*theta) .*((-sin(x4.*theta)).*((1 + G2.*G3.*Z072) .* 
cos(x5.*theta) + 11.*Z0.*(G2 + G4 + G2.+*G3.+*G4.*Z0°2) .*sin(x5.*theta) ) 
- cos(x4.*theta) .*((-1i) .*G2.*Z0.*cos(x5.*theta) + (1 + G2.*(G3 + G4) 
.*Z0^2) .*sin(x5.*theta))))))))./(cos(x1.*theta).*(cos(x2.*theta).*(cos 
(x3.*theta) .*(cos(x4.*theta) .*(cos(x5.*theta) + 11.*G5.*Z0.*sin(x5.* 
theta)) - sin(x4.*theta).*((-1i).*(G4 + G5).*Z0.*cos(x5.*theta) + (1 + 
G4.*G5.*Z0°2) .*sin(x5.*theta))) + 1i.*sin(x3.*theta).*(cos(x4.*theta) 
.*((G3 + G4 + G5).*Z0.*cos(x5.*theta) + 11.*(1 + G3.*G5.*Z0°2 + G4.*G5 
.*Z072) .*sin(x5.*theta)) + 1i.*sin(x4.*theta).*((1 + G3.*(G4 + G5) .*Z0 
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72) .x*cos(x5.*theta) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z0°2) .*sin(xb.* 
theta)))) + 1i.*sin(x2.*theta).*(cos(x3.*theta).*(cos(x4.*theta).*((G2 
+ G3 + G4 + G5).*Z0.*cos(x5.*theta) + 1i.*(1 + G2.*G5.*Z0°2 + G3.*G5 

.*Z072 + 
G4.*G5.*Z0°2).*sin(x5.*theta)) + 1i.*sin(x4.*theta).*((1 + G2.*(G4 + G5) 
.*Z0^2 + G3.*(G4 + G5) .*Z0°2).*cos(x5.*theta) + 1i.*Z0.*(G2 + G3 + G5 
+ G2.*G4.*G5.*Z0°2 + G3.*G4.*G5.*Z0^2) .*sin(x5.*theta))) + 1i.*sin(x3 
.*theta).*(cos(x4.*theta).*((1 + G2.*(G3 + G4 + G5) .*Z072) .*cos (x5.* 
theta) + 11.*Z0.*(G2 + G5 + G2.*G3.*G5.*Z0°2 + G2.*G4.*G5.*Z072) .*sin( 
xb.*theta)) - sin(x4.*theta).*((-1i).*Z0.*(G4 + G5 + G2.*(1 + G3.*(G4 
+ G5).*Z0°2)).*cos(x5.*theta) + (1 + G4.*G5.*Z0°2 + G2.*Z0^2.*(G3 + G5 
+ G3.*G4.*G5.*Z0°2)).*sin(x5.*theta))))) + 1i.*sin(x1.*theta).*(cos( 
x2.*theta).*(cos(x3.*theta).*(cos(x4.*theta).*((G1 + G2 + G3 + G4 + G5 
).*Z0.*cos(xb.*theta) + 1i.*(1 + G1.*G5.*Z0^2 + G2.*G5.*Z072 + G3.*G5 
.*Z0^2 + G4.*G5.*Z0°2) .*sin(xb.*theta)) + 1i.*sin(x4.*theta).*((1 + G3 
.*G4.*Z2072 + G3.*G5.*Z0"2 + G1.*(G4 + 
G5) .*Z0°2 + G2.*(G4 + G5).*Z0^2).*cos(xb.*theta) + 1i.*Z0.*(G1 + G2 + G3 
+ G5 + G1.*G4.*G5.*Z2072 + G2.*G4.*G5.*Z0°2 + G3.*G4.*G5.*Z0°2) .*sin(x5 
.*theta))) + 1i.*sin(x3.*theta) .*(cos(x4.*theta).*((1 + G1.*(G3 + G4 + 
GB) .*2072 + G2.*(G3 + G4 + G5).*Z0°2).*cos(x5.*theta) + 1i.*Z0.*(G1 + 
G2 + Gb + G1.*G3.*G5.*Z20^2 + G2.*G3.*G5.*Z0°2 + G1.*G4.*G5.*Z072 + G2 
.*G4.*G5.*Z2072) .*sin(x5.*theta)) - sin(x4.*theta) .*((-1i).*Z0.*(G2 + 
G4 + G5 + G2.*G3.*G4.*Z0°2 + G2.*G3.*G5.*Z0°2 + G1.*(1 + G3.*(G4 + G5) 
.*Z0°2)).*cos(x5.*theta) + (1 + G4.*G5.*Z0°2 + G1.*Z0^2.*(G3 + G5 + G3 
.*G4.*G5.*Z0°2) + G2.*Z0°2.*(G3 + G5 + G3.*G4.*G5.*Z0°2)).*sin(x5.* 
theta)))) - sin(x2.*theta).*(cos(x3.*theta).*((-1i).*cos(x4.*theta) 
.*((1 + G1.*(G2 + G3 + G4 + G5).*Z0°2).*cos(x5.*theta) + 1i.*Z0.*(G5 + 
G1.*(1 + G2.*G5.*Z0°2 + G3.+*G5.*Z0°2 + G4.*G5.*Z072)).*sin(x5.*theta) 
) + sin(x4.*theta).*(Z0.*(G4 + G5 + Gi.*(1 + G2.*(G4 + G5) .*Z072 + 
G3.*(G4 + G5).*Z0°2)).*cos(x5.*theta) + 1i.*(1 + G4.*G5.*Z0"2 + G1.*Z0 
72.*(G2 + G3 + G5 + G2.*G4.*G5.*Z0°2 + G3.*G4.*G5.*Z072)) .*sin(x5.* 
theta))) + sin(x3.*theta).*(cos(x4.*theta).*(Z0.*(G3 + G4 + G5 + G1 
.*(1 + G2.*(G3 + G4 + G5).*Z072)).*cos(x5.*theta) + 1i.*(1 + G3.*Gb.x* 
Z0"2 + G4.*G5.*Z0°2 + G1.*Z0^2.*(G2 + G5 + G2.*G3.*G5.*Z0"2 + G2.*G4.* 
G5.*Z0°2)).*sin(x5.*theta)) + 1li.*sin(x4.*theta).*((1 + G3.*(G4 + G5) 
.*Z0^2 + G1.*Z0^2.*(G4 + G5 + G2.*(1 + G3.*(G4 + G5) .*Z0^2))) .*cos (x5 
.*theta) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z0°2 + G1.*(1 + G4.*G5.*Z0^2 
+ G2.*Z0°2.*(G3 + G5 + G3.+*G4.+*G5.*Z0°2))).*sin(x5.*theta))))))))); 
case 6 
G1 = individual (1); 
G2 = individual (2); 
G3 = individual (3); 
G4 = individual (4); 
G5 = individual (5); 
G6 = individual (6) ; 
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x1 = individual(7); 
x2 = individual(8); 
x3 = individual(9); 
x4 = individual (10); 
xb = individual (11); 
x6 = individual (12) ; 
gamma = mag2db(abs((-ZS + (11.*Z0.*(cos(x1.* (theta) ) .*(cos (x2. zt 
theta) ).*(cos(x3.* (theta) ).*(cos(x5.* (theta) ).*(cos(x6.* (theta) ) 
.*sin(x4.*(theta)) + (cos(x4.*(theta)) + 11.*(G4 + G5) .*Z0.*sin( 
x4.*(theta))).*sin(x6.*(theta))) + sin(x5.*(theta)).*(cos (x4. zt 
theta) ).*(cos(x6.*(theta)) + 1i.*G5.*Z0.*sin(x6.*(theta))) - sin 
(x4.* (theta) ).*((-1i) .*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.+*G5.* 
Z20^2).*sin(x6.*(theta))))) + sin(x3.*(theta)).*(cos(x4.*(theta)) 
.*(cos(x5.*(theta)).*(cos(x6.*(theta)) + 1i.*(G3 + G4 + G5) .*Z0 
.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4).*Z0.* 
cos(x6.*(theta)) + (1 + G3.*G5.*Z0^2 + G4.*G5.*Z072) .*sin(x6.*( 
theta)))) + sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + G3.*G4 
.*Z0^2).*cos(x6.*(theta)) + 11.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) 
.*sin(x6.*(theta))) - 
cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5) .* 
20^2).*sin(x6.*(theta)))))) + sin(x2.*(theta)).*((-sin(x3.*(theta))) 
.*(cos(x5.*(theta)).*(cos(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*( 
theta)) + 1i.*(G4 + G5).*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + 
sin(x5.*(theta)).*(cos(x4.*(theta)).*(cos(x6.*(theta)) + 11.*G5.*Z0.* 
sin(x6.*(theta))) - sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) 
+ (1 + G4.*G5.*Z0°2) .*sin(x6.*(theta))))) + cos(x3.*(theta)).*(cos(x4 
.*(theta)).*(cos(xb.*(theta)).*(cos(x6.*(theta)) + 1i.*(G3 + G4 + G5) 
.*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4).*Z0.*cos 
(x6.*(theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z0°2) .*sin(x6.*(theta)))) 
+ sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + G3.*G4.*Z0^2) .*cos (x6 
.*(theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z0°2) .*sin(x6.*(theta))) - 
cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + 
G5).*Z0^2) .*sin(x6.*(theta))))) + 1i.*G2.*Z0.*(cos(x3.*(theta)).*(cos( 
xb.*(theta)).*(cos(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) 
+ 11.*(G4 + G5).*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + 
sin(x5.*(theta)).*(cos(x4.*(theta)).*(cos(x6.*(theta)) + 11.*G5.*Z0.*sin( 
x6.*(theta))) - sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + 
(1 + G4.*G5.*Z20^2).*sin(x6.*(theta))))) + sin(x3.*(theta)) .*(cos (x4. zt 
theta) ) .*(cos(x5.*(theta)).*(cos(x6.*(theta)) + 1i.*(G3 + G4 + G5) .*Z0 
.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4) .*Z0.*cos(x6 
.*(theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z072) .*sin(x6.*(theta)))) + 
sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + G3.*G4.*Z0°2) .*cos (x6. *( 
theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) .*sin(x6.*(theta))) - 
cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5 
).*Z072) .*sin(x6.*(theta)))))))) + sin(x1.*(theta)).*((-sin(x2.* (theta 
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))) .* (cos (x3.* (theta) ) .*(cos(x5.* (theta) ).*(cos(x6.*(theta)).*sin(x4 
.*(theta)) + (cos(x4.*(theta)) + 11.*(G4 + G5) .*Z0.*sin(x4.*(theta))) 
.*sin(x6.*(theta))) + 

sin(x5.*(theta)) .*(cos(x4.* (theta) ).*(cos(x6.*(theta)) + 11.*G5.*Z0.*sin( 
x6.*(theta))) - sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + 
(1 + G4.*G5.*Z0°2).*sin(x6.*(theta))))) + sin(x3.*(theta)) .*(cos (x4. zt 
theta) ) .*(cos(x5.*(theta)).*(cos(x6.*(theta)) + 1i.*(G3 + G4 + G5) .*Z0 
.*sin(x6.*(theta))) - sin(xb.*(theta)).*((-1i).*(G3 + G4) .*Z0.*cos(x6 
.*(theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z072) .*sin(x6.*(theta)))) + 
sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + G3.*G4.*Z0°2) .*cos (x6. *( 
theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) .*sin(x6.*(theta))) - 
cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5 
).*Z072) .*sin(x6.*(theta)))))) + cos(x2.*(theta)).*((-sin(x3.* (theta) ) 
).* (cos (x5.* (theta) ) .*(cos(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*( 
theta)) + 1i.*(G4 + G5).*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + 
sin(x5.*(theta)).*(cos(x4.*(theta)).*(cos(x6.*(theta)) + 11.*G5.*Z0.* 
sin(x6.*(theta))) - 

sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5.*Z072) .* 
sin(x6.*(theta))))) + cos(x3.*(theta)).*(cos(x4.* (theta) ).*(cos (x5. zt 
theta)).*(cos(x6.*(theta)) + 11.*(G3 + G4 + G5) .*Z0.*sin(x6.* (theta) )) 
- sin(x5.*(theta)).*((-1i).*(G3 + G4) .*Z0.*cos(x6.*(theta)) + (1 + G3 
.*G5.*Z0°2 + G4.*G5.*Z0°2).*sin(x6.*(theta)))) + sin(x4.*(theta)).*((- 
sin(x5.*(theta))).*((1 + G3.*G4.*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 
+ G5 + G3.*G4.*G5.*Z0°2) .*sin(x6.*(theta))) - cos(x5.*(theta)).*((-1i 
).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5) .*Z072) .*sin(x6.*( 
theta))))) + 11.*G2.*Z0.*(cos(x3.* (theta) ) .*(cos(x5.* (theta) ) .*(cos(x6 
.*(theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) + 1i.*(G4 + G5).*Z0.* 
sin(x4.*(theta))).*sin(x6.*(theta))) + sin(x5.*(theta)).*(cos (x4. zt 
theta)).*(cos(x6.*(theta)) + 1i.*G5.*Z0.*sin(x6.*(theta))) - sin(x4.*( 
theta)) .*((-1i) .*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5.*Z072) .*sin(x6 
.*(theta))))) + 

sin(x3.*(theta)) .*(cos(x4.* (theta) ).*(cos(x5.* (theta) ) .*(cos(x6.* (theta) ) 
+ 11.*(G3 + G4 + G5).*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i 
).*(G3 + G4) .*Z0.*cos(x6.*(theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z072) 
.*sin(x6.*(theta)))) + sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + 
G3.*G4.*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) 
.*sin(x6.*(theta))) - cos(x5.*(theta)).*((-1i) .*G3.*Z0.*cos(x6.* (theta 
)) + (1 + G3.*(G4 + G5) .*Z072) .*sin(x6.*(theta))))))) + 1i.*G1.*Z0.*( 
cos (x2. * (theta) ) .*(cos(x3.*(theta)) .*(cos(x5.* (theta) ) .*(cos (x6. *( 
theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) + 1i.*(G4 + G5).*Z0.*sin 
(x4.* (theta) )).*sin(x6.*(theta))) + sin(x5.*(theta)).*(cos(x4.*(theta) 
).*(cos(x6.*(theta)) + 1i.*G5.*Z0.*sin(x6.*(theta))) - sin(x4.*(theta) 
).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5.*Z072) .*sin(x6.*( 
theta))))) + sin(x3.*(theta)).*(cos(x4.* (theta) ).*(cos(x5.* (theta) ) .*( 
cos(x6.*(theta)) + 1i.*(G3 + G4 + 
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GB) .*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4).*Z0.*cos 
(x6.*(theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z072) .*sin(x6.* (theta) ))) 
+ sin(x4.*(theta)).*((-sin(xb.*(theta))).*((1 + G3.*G4.*Z0^2) .*cos (x6 
.*(theta)) + 11.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) .*sin(x6.*(theta))) - 

cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + 
GB) .*Z0°2) .*sin(x6.*(theta)))))) + sin(x2.*(theta)).*((-sin(x3.*(theta 
))) .* (cos (x5.* (theta) ) .*(cos(x6.* (theta) ).*sin(x4.*(theta)) + (cos(x4 
.*(theta)) + 1i.*(G4 + G5).*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + 

sin(x5.*(theta)).*(cos(x4.*(theta)).*(cos(x6.*(theta)) + 11.*G5.*Z0.* 
sin(x6.*(theta))) - sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) 

+ (1 + G4.*G5.*Z0°2) .*sin(x6.*(theta))))) + cos(x3.*(theta)).*(cos(x4 
.*(theta)).*(cos(x5.*(theta)).*(cos(x6.*(theta)) + 1i.*(G3 + G4 + G5) 
.*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4).*Z0.*cos 
(x6.*(theta)) + (1 + 

G3.*G5.*Z0^2 + G4.*G5.*Z072) .*sin(x6.*(theta)))) + sin(x4.*(theta)).*((- 
sin(x5.*(theta))).*((1 + G3.*G4.*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 
+ G5 + G3.*G4.*G5.*Z0°2) .*sin(x6.*(theta))) - cos(x5.*(theta)).*((-1i 
).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5) .*Z072) .*sin(x6.*( 
theta))))) + 11.*G2.*Z0.*(cos(x3.* (theta) ) .*(cos(x5.* (theta) ) .*(cos (x6 
.*(theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) + 1i.*(G4 + G5).*Z0.* 
sin(x4.*(theta))).*sin(x6.*(theta))) + sin(x5.*(theta)).*(cos (x4. zt 
theta)) .*(cos(x6.*(theta)) + 1i.*G5.*Z0.*sin(x6.*(theta))) - sin(x4.*( 
theta)) .*((-1i) .*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5.*Z072) .*sin(x6 
.*(theta))))) + sin(x3.*(theta)).*(cos(x4.*(theta)).*(cos(x5.*(theta)) 
.*(cos(x6.*(theta)) + 1i.*(G3 + G4 + G5).*Z0.*sin(x6.*(theta))) - Sint 
xb.*(theta)).*((-1i).*(G3 + G4) .*Z0.*cos(x6.*(theta)) + (1 + G3.*G5.* 
Z0"2 + G4.*G5.*Z072) .*sin(x6.*(theta)))) + sin(x4.*(theta)).*((-sin(x5 
.*(theta))).*((1 + 

G3.*G4.*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) .* 
sin(x6.*(theta))) - cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) 
+ (1 + G3.*(G4 + G5).*Z072) .*sin(x6.*(theta)))))))))))./(cos(x1.*( 
theta)).*(cos(x2.*(theta)).*(cos(x3.*(theta)).*(cos(x4.*(theta)).*(cos 
(x5.*(theta)).*(cos(x6.*(theta)) + 1i.*G6.*Z0.*sin(x6.*(theta))) - sin 
(x5.*(theta)).*((-1i).*(G5 + G6).*Z0.*cos(x6.*(theta)) + (1 + G5.*G6.* 
20^2).*sin(x6.*(theta)))) + 1i.*sin(x4.*(theta)).*(cos(x5.*(theta)) 
.*((G4 + G5 + G6).*Z0.*cos(x6.*(theta)) + 1i.*(1 + G4.*G6.*Z072 + G5.* 
G6.*Z2072) .*sin(x6.*(theta))) + 1i.*sin(x5.*(theta)).*((1 + G4.*(G5 + 
G6) .*Z0^2) .*cos(x6.*(theta)) + 11.*Z0.*(G4 + G6 + G4.*G5.*G6.*Z072).x* 
sin(x6.*(theta))))) + 1i.*sin(x3.*(theta)).*(cos(x4.*(theta)).*(cos(x5 
.*(theta)) .*((G3 + G4 + G5 + G6).*Z0.*cos(x6.*(theta)) + 1i.*(1 + G3.x* 
G6.*Z0°2 + G4.*G6.*Z0°2 + G5.*G6.*Z072) .*sin(x6.*(theta))) + 1i.*sin( 
x5.*(theta)).*((1 + 

G3.*(G5 + G6) .*Z072 + G4.*(G5 + G6) .*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*( 
G3 + G4 + G6 + G3.*G5.*G6.*Z0^2 + G4.*G5.*G6.*Z0°2) .*sin(x6.* (theta) ) ) 
) + 1i.*sin(x4.* (theta) ).*(cos(x5.*(theta)).*((1 + G3.*(G4 + G5 + G6) 
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.*Z0^2).*cos(x6.*(theta)) + 1i.*Z0.*(G3 + G6 + G3.*G4.*G6.*20^2 + G3. 
G5.*G6.*Z072).*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*Z0.*(G5 
G6 + G3.*(1 + G4.*(G5 + G6).*Z0°2)).*cos(x6.*(theta)) + (1 + G5.*G6. 
ZO°2 + G3.*Z07"2.*(G4 + G6 + G4.*G5.*G6.*Z072)).*sin(x6.*(theta)))))) 
1i.*sin(x2.*(theta)).*(cos(x3.*(theta)).*(cos(x4.*(theta)).*(cos(x5 
.*(theta)).*((G2 + G3 + G4 + G5 + G6).*Z0.*cos(x6.*(theta)) + 1i.*(1 + 
G2.*G6.*Z0°2 + G3.*G6.*Z0°2 + G4.*G6.*Z0"2 + G5.*G6.*Z072) .*sin(x6.*( 
theta))) + 1i.*sin(x5.*(theta)).*((1 + G4.*G5.*Z0°2 + G4.*G6.*20^2 + 
G2.*(G5 + G6) .*Z072 + G3.*(G5 + G6) .*Z072) .*cos(x6.*(theta)) + 1i.*Z0 
.*(G2 + G3 + G4 + G6 + 
G2.*G5.*G6.*Z072 + G3.*G5.*G6.*Z072 + G4.*G5.*G6.*Z072).*sin(x6.*(theta)) 
)) + 1i.*sin(x4.*(theta)) .*(cos(x5.*(theta)).*((1 + G2.*(G4 + Gb + G6) 
.*Z0^2 + G3.*(G4 + G5 + G6).*Z0°2).*cos(x6.*(theta)) + 1i.*Z0.*(G2 + 
G3 + G6 + G2.*G4.*G6.*Z20^2 + G3.*G4.*G6.*Z0°2 + G2.*G5.*G6.*Z072 + G3 
.*G5.*G6.*Z072) .*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*Z0.*(G3 
+ G5 + G6 + G3.*G4.*G5.*Z0°2 + G3.*G4.*G6.*Z072 + G2.*(1 + G4.*(G5 + 
G6) .*Z0°2)) .*cos(x6.*(theta)) + (1 + G5.*G6.*Z0°2 + G2.*Z0^2.*(G4 + G6 
+ G4.*G5.*G6.*Z0°2) + G3.*Z0^2.*(G4 + G6 + G4.*G5.*G6.*Z072)).*sin(x6 
.*(theta))))) - sin(x3.*(theta)).*(cos(x4.* (theta) ) .*((-1i) .*cos (x5. zt 
theta)).*((1 + G2.*(G3 + G4 + G5 + G6) .*Z072) .*cos(x6.*(theta)) + 1i.* 
Z0.*(G6 + G2.*(1 + G3.*G6.*Z0"2 + G4.*G6.*Z0°2 + G5.*G6.*Z072)) .*sin( 
x6.*(theta))) + sin(x5.*(theta)).*(Z0.*(G5 + G6 + G2.*(1 + G3.*(G5 + 
G6) .*Z072 + 
G4.*(G5 + G6).*Z20^2)) .*cos(x6.*(theta)) + 1i.*(1 + G5.*G6.*Z072 + G2.*Z0 
72.*(G3 + G4 + G6 + G3.*G5.*G6.*Z072 + G4.+*G5.*G6.*Z0°2)) .*sin(x6.*( 
theta)))) + sin(x4.*(theta)).*(cos(x5.*(theta)).*(Z0.*(G4 + G5 + G6 + 
G2.*(1 + G3.*(G4 + G5 + G6).*Z0^2)).*cos(x6.*(theta)) + 1i.*(1 + G4.x* 
G6.*Z20^2 + G5.*G6.*Z20^2 + G2.*Z0°2.*(G3 + G6 + G3.*G4.*G6.*Z0°2 + G3.x* 
G5.*G6.*Z072)).*sin(x6.*(theta))) + 1i.*sin(x5.*(theta)).*((1 + G4.*( 
G5 + G6).*Z0°2 + G2.*Z0°2.*(G5 + G6 + G3.*(1 + G4.*(G5 + G6) .*Z072))) 
.*cos(x6.*(theta)) + 1i.*Z0.*(G4 + G6 + G4.*G5.*G6.*Z072 + G2.*(1 + G5 
.*G6.*Z0^2 + G3.*Z0°2.*(G4 + G6 + G4.*G5.*G6.*Z0°2))) .*sin(x6.* (theta) 
11333) + 1i.*Z0.*sin(x1.*(theta)).*((1./Z0).*(sin(x2.*(theta)).*(cos( 
x3.*(theta)).*(1i.*cos(x4.*(theta)).*(cos(x5.*(theta)).*(cos(x6.*( 
theta)) + 1i.*G6.*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*( 
G5 + G6).*Z0.*cos(x6.*(theta)) + (1 + 
G5.*G6.*Z072).*sin(x6.*(theta)))) + sin(x4.*(theta)).*((-1i).*cos(xb.*( 
theta)).*((-1i).*(G4 + G5 + G6).*Z0.*cos(x6.*(theta)) + (1 + G4.*G6.* 
Z0"2 + G5.*G6.*Z072) .*sin(x6.*(theta))) + sin(x5.*(theta)).*((-1i).*(1 
+ G4.*(G5 + G6) .*Z072) .*cos(x6.*(theta)) + Z0.*(G4 + G6 + G4.*G5.*G6 
.*Z0^2).*sin(x6.*(theta))))) + sin(x3.*(theta)) .*((-1i) .*cos (x4. *( 
theta) ).*(cos(x5.* (theta) ).*((-1i).*(G3 + G4 + G5 + G6).*Z0.*cos (x6. *( 
theta)) + (1 + G3.*G6.*Z0°2 + G4.*G6.*Z0°2 + G5.*G6.*Z072) .*sin(x6.*( 
theta))) + sin(x5.*(theta)).*((1 + G3.*(G5 + G6).*Z0°2 + G4.*(G5 + G6) 
.*Z0^2).*cos(x6.*(theta)) + 11.*Z0.*(G3 + G4 + G6 + G3.*Gb.*G6.*Z0"2 + 


+ * + %* 


82 


G4.*G5.*G6.*Z072) .*sin(x6.*(theta)))) + sin(x4.*(theta)).*(cos(xb.*( 
theta)).*((-1i).*(1 + G3.*(G4 + GB + G6).*Z0^2) .*cos(x6.*(theta)) + ZO 
.*(G3 + G6 + G3.*G4.*G6.*Z0°2 + G3.*G5.*G6.*Z072) .*sin(x6.*(theta))) + 

sin(x5.*(theta)).*(ZO.*(G5 + G6 + G3.*(1 + G4.*(G5 + G6).*Z0°2)).*cos 

(x6.*(theta)) + 1i.*(1 + G5.*G6.*Z072 + G3.*Z0^2.*(G4 + G6 + 

G4.*G5.*G6.*Z072)) .*sin(x6.*(theta))))))) + (1./Z0) .*(cos(x2.* (theta) ) .*( 
cos(x3.* (theta) ) .*(cos(x4.* (theta) ) .*(cos(x5.*(theta)).*((G2 + G3 + G4 

+ G5 + G6).*Z0.*cos(x6.*(theta)) + 1i.*(1 + G2.*G6.*Z0°2 + G3.*G6.*Z0 
^2 + G4.*G6.*Z0^2 + G5.*G6.*Z072) .*sin(x6.*(theta))) + 1i.*sin(x5.*( 
theta)).*((1 + G4.*G5.*Z0"2 + G4.*G6.*Z0°2 + G2.*(G5 + G6).*Z0°2 + G3 
.*(G5 + G6) .*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G2 + G3 + G4 + G6 + G2 
.*G5.*G6.*Z0°2 + G3.*G5.*G6.*Z0"2 + G4.*G5.*G6.*Z072) .*sin(x6.* (theta) 
))) + 1i.*sin(x4.* (theta) ).*(cos(x5.*(theta)).*((1 + G2.*(G4 + G5 + Gë 
).*Z0^2 + G3.*(G4 + G5 + G6).*Z0^2).*cos(x6.*(theta)) + 1i.*Z0.*(G2 + 
G3 + G6 + G2.*G4.*G6.*Z072 + G3.*G4.*G6.*Z0°2 + G2.*G5.*G6.*Z072 + G3 
.*G5.*G6.*Z072) .*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*Z0.*(G3 
+ G5 + G6 + G3.*G4.*G5.*Z0°2 + G3.*G4.*G6.*Z072 + G2.*(1 + G4.*(G5 + 
G6) .*Z0^2)) .*cos(x6.*(theta)) + (1 + 

G5.*G6.*Z072 + G2.*Z0°2.*(G4 + G6 + G4.*G5.*G6.*Z072) + G3.*Z0"2.*(G4 + 
G6 + G4.*G5.*G6.*Z072)).*sin(x6.*(theta))))) - sin(x3.*(theta)) .*(cos( 
x4.*(theta)).*((-1i).*cos(x5.*(theta)).*((1 + G2.*(G3 + G4 + G5 + G6) 

.*Z0^2).*cos(x6.*(theta)) + 1i.*Z0.*(G6 + G2.*(1 + G3.*G6.*Z0"2 + G4.* 
G6.*Z0°2 + G5.*G6.*Z072)).*sin(x6.*(theta))) + sin(x5.*(theta)).*(Z0 
.*(G5 + G6 + G2.*(1 + G3.*(G5 + G6).*Z0°2 + G4.*(G5 + G6).*Z072)).*cos 
(x6.*(theta)) + 1i.*(1 + G5.*G6.*Z072 + G2.*Z0°2.*(G3 + G4 + G6 + G3.* 
G5.*G6.*Z072 + G4.*G5.*G6.*Z0°2)).*sin(x6.*(theta)))) + sin(x4.*(theta 
)) .*(cos(x5.*(theta)).*(Z0.*(G4 + G5 + G6 + G2.*(1 + G3.*(G4 + G5 + Gë 
).*Z072)).*cos(x6.*(theta)) + 1i.*(1 + G4.*G6.*Z0°2 + G5.*G6.*Z072 + 
G2.*Z0°2.*(G3 + G6 + G3.*G4.*G6.*Z20^2 + G3.*G5.*G6.*Z072)) .*sin(x6.*( 
theta))) + 1i.*sin(x5.*(theta)).*((1 + G4.*(G5 + 

G6) .*Z0°2 + G2.*Z0°2.*(G5 + G6 + G3.*(1 + G4.*(G5 + G6).*Z072))).*cos(x6 

.*(theta)) + 1i.*Z0.*(G4 + G6 + G4.*G5.*G6.*Z0°2 + G2.*(1 + G5.*G6.*Z0 
72 + G3.*Z20^2.*(G4 + G6 + G4.*G5.*G6.*Z0°2))).*sin(x6.*(theta))))))) + 
G1.*(cos(x2.* (theta) ) .*(cos(x3.* (theta) ).*(cos(x4.* (theta) ) .*(cos (x5 
.* (theta) ).*(cos(x6.*(theta)) + 1i.*G6.*Z0.*sin(x6.*(theta))) - sin(x5 
.*(theta)).*((-1i).*(G5 + G6) .*Z0.*cos(x6.*(theta)) + (1 + G5.*G6.*Z0 
72) .*sin(x6.*(theta)))) + 1i.*sin(x4.*(theta)).*(cos(x5.*(theta)).*(( 
G4 + GB + G6).*Z0.*cos(x6.*(theta)) + 1i.*(1 + G4.*G6.*Z0°2 + G5.*G6.* 
20^2).*sin(x6.*(theta))) + 1i.*sin(x5.*(theta)).*((1 + G4.*(G5 + G6).* 
20^2).*cos(x6.*(theta)) + 1i.*Z0.*(G4 + G6 + G4.*G5.*G6.*Z072) .*sin(x6 
.*(theta))))) + 1i.*sin(x3.*(theta)).*(cos(x4.* (theta) ) .*(cos(x5. zt 
theta)).*((G3 + G4 + G5 + G6).*Z0.*cos(x6.*(theta)) + 

1i.*(1 + G3.*G6.*Z0"2 + G4.*G6.*Z20^2 + G5.*G6.*Z0°2) .*sin(x6.*(theta))) + 
1i.*sin(x5.*(theta)).*((1 + G3.*(G5 + G6).*Z0°2 + G4.*(G5 + G6).*Z0 
72) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 + G4 + G6 + G3.*G5.*G6.*Z0"2 + G4 
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.*G5.*G6.*Z0°2) .*sin(x6.*(theta)))) + 1i.*sin(x4.*(theta)).*(cos(x5.*( 
theta)).*((1 + G3.*(G4 + G5 + G6) .*Z072) .*cos(x6.*(theta)) + 1i.*ZO.*( 
G3 + G6 + G3.*G4.*G6.*Z072 + G3.*G5.*G6.*Z072) .*sin(x6.*(theta))) - 
sin(x5.*(theta)).*((-1i).*Z0.*(G5 + G6 + G3.*(1 + G4.*(G5 + G6).*Z072) 
).*cos(x6.*(theta)) + (1 + G5.*G6.*Z0^2 + G3.*Z0^2.*(G4 + G6 + G4.*G5 
.*G6.*Z0°2)).*sin(x6.*(theta)))))) + 1i.*sin(x2.*(theta)).*(cos(x3.*( 
theta)).*(cos(x4.*(theta)).*(cos(x5.*(theta)).*((G2 + G3 + G4 + G5 + 
G6) .*Z0.*cos(x6.*(theta)) + 1i.*(1 + G2.*G6.*Z0°2 + G3.*G6.*Z0"2 + G4 
.*G6.*Z0^2 + G5.*G6.*Z20^2) .*sin(x6.*(theta))) + 1i.*sin(x5.*(theta)) 
A(( + G4.*G5.*Z0°2 + G4.*G6.*Z072 + G2.*(G5 + G6).*Z0°2 + G3.*(G5 + 
G6) .*Z072) .*cos(x6.*(theta)) + 11.*Z0.*(G2 + G3 + G4 + G6 + 
G2.*G5.*G6.*Z072 + G3.*G5.*G6.*Z072 + G4.*G5.*G6.*Z0°2) .*sin(x6.* (theta) ) 
)) + 1i.*sin(x4.*(theta)).*(cos(x5.*(theta)).*((1 + G2.*(G4 + G5 + G6) 
.*Z072 + G3.*(G4 + GB + G6).*Z072).*cos(x6.*(theta)) + 11.*Z0.*(G2 + 
G3 + G6 + G2.*G4.*G6.*Z072 + G3.*G4.*G6.*Z0°2 + G2.*G5.*G6.*Z072 + G3 
.*G5.*G6.*Z0°2) .*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*Z0.*(G3 
+ G5 + G6 + G3.*G4.*G5.*Z0°2 + G3.*G4.*G6.*Z072 + G2.*(1 + G4.*(G5 + 
G6) .*Z072)) .*cos(x6.*(theta)) + (1 + G5.*G6.*Z072 + G2.*Z0^2.*(G4 + G6 
+ G4.*G5.*G6.*Z0°2) + G3.*Z0°2.*(G4 + G6 + G4.*G5.*G6.*Z072)) .*sin(x6 
.*(theta))))) - sin(x3.*(theta)).*(cos(x4.*(theta)) .*((-1i).*cos(x5.*( 
theta)).*((1 + G2.*(G3 + G4 + G5 + G6).*Z2072).*cos(x6.*(theta)) + 1i.* 
Z0.*(G6 + G2.*(1 + G3.*G6.*Z0^2 + G4.*G6.*Z0^2 + G5.*G6.*Z072)).* sin( 
x6.*(theta))) + sin(x5.*(theta)).*(Z0.*(G5 + G6 + G2.*(1 + G3.*(G5 + 
G6) .*Z072 + 
G4.*(G5 + G6) .*Z0°2)) .*cos(x6.*(theta)) + 1i.*(1 + G5.*G6.*Z072 + G2.*Z0 
72.*(G3 + G4 + G6 + G3.*G5.*G6.*Z0"2 + G4.*G5.*G6.*Z072)) .*sin(x6.*( 
theta)))) + sin(x4.*(theta)).*(cos(x5.*(theta)).*(Z0.*(G4 + G5 + G6 + 
G2.*(1 + G3.*(G4 + GB + G6).*Z072)).* cos(x6.*(theta)) + 1i.*(1 + G4.x* 
G6.*Z20^2 + G5.*G6.*Z20^2 + G2.*Z0°2.*(G3 + G6 + G3.*G4.*G6.*Z0°2 + G3.x* 
G5.*G6.*Z0°2)) .*sin(x6.*(theta))) + 1i.*sin(x5.*(theta)).* ((1 + G4.*( 
G5 + G6).*Z0°2 + G2.*Z0°2.*(G5 + G6 + G3.*(1 + G4.*(G5 + G6) .*Z072))) 
.*cos(x6.*(theta)) + 1i.*Z0.*(G4 + G6 + G4.*G5.*G6.*Z072 + G2.*(1 + G5 
.*G6.*Z0°2 + G3.*Z0^2.*(G4 + G6 + G4.*G5.*G6.*Z072))) .*sin(x6.* (theta) 
)))))))))./ (ZS + (1i.*Z0.* (cos (x1.* (theta) ) .*(cos(x2.* (theta) ).*(cos( 
x3.*(theta)).*(cos(x5.*(theta)).*(cos(x6.*(theta)).*sin(x4.*(theta)) + 
(cos(x4.*(theta)) + 11.*(G4 + G5) .*Z0.*sin(x4.*(theta))).*sin(x6.*( 
theta))) + 
sin(x5.*(theta)) .*(cos(x4.* (theta) ).*(cos(x6.*(theta)) + 11.*G5.*Z0.*sin( 
x6.*(theta))) - sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + 
(1 + G4.*G5.*Z072) .*sin(x6.*(theta))))) + sin(x3.*(theta)) .*(cos (x4. zt 
theta) ).*(cos(x5.*(theta)).*(cos(x6.*(theta)) + 1i.*(G3 + G4 + G5) .*Z0 
.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4) .*Z0.*cos(x6 
.*(theta)) + (1 + G3.*G5.*Z0"2 + G4.*G5.*Z072) .* sin(x6.*(theta)))) + 
sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + G3.*G4.*Z0^2) .*cos (x6. *( 
theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) .*sin(x6.*(theta))) - 
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cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5 
).*Z072) .*sin(x6.*(theta)))))) + sin(x2.*(theta)).*((-sin(x3.*(theta)) 
).*(cos(x5.*(theta)).*(cos(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*( 
theta)) + 1i.*(G4 + G5).*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + 
sin(x5.*(theta)).*(cos(x4.*(theta)).*(cos(x6.*(theta)) + 11.*G5.*Z0.* 
sin(x6.*(theta))) - 

sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5.*Z072) .* 
sin(x6.*(theta))))) + cos(x3.*(theta)).*(cos(x4.*(theta)).*(cos(x5.*( 
theta)).*(cos(x6.*(theta)) + 11.*(G3 + G4 + G5).*ZO.*sin(x6.*(theta))) 
- sin(x5.*(theta)).*((-1i).*(G3 + G4).*Z0.*cos(x6.*(theta)) + (1 + G3 
.*G5.*Z0^2 + G4.*G5.*Z07"2).* sin(x6.*(theta)))) + sin(x4.*(theta)) 
.*((-sin(x5.*(theta))).*((1 + G3.*G4.*Z0°2) .*cos(x6.*(theta)) + 1i.*Z0 
.* (G3 + G5 + G3.*G4.*G5.*Z2072) .*sin(x6.*(theta))) - cos(x5.*(theta)) 
(Eli) .*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5) .*Z072) .*sin(x6 
.*(theta))))) + 11.*G2.*Z0.*(cos(x3.* (theta) ).*(cos(x5.*(theta)).*(cos 
(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) + 1i.*(G4 + G5).x* 
Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + sin(x5.*(theta)).*(cos(x4 
.*(theta)).*(cos(x6.*(theta)) + 11.*G5.*Z0.*sin(x6.*(theta))) - sin(x4 
.* (theta) ) .*((-1i) .*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5.*Z072) .*sin 
(x6.*(theta))))) + 

sin(x3.*(theta)) .*(cos(x4.* (theta) ).*(cos(x5.* (theta) ) .*(cos(x6.* (theta) ) 
+ 11.*(G3 + G4 + G5).*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i 
).*(G3 + G4) .*Z0.*cos(x6.*(theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z072) 
.*sin(x6.*(theta)))) + sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + 
G3.*G4.*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) 
.*sin(x6.*(theta))) - cos(x5.*(theta)).*((-1i) .*G3.*Z0.*cos(x6.* (theta 
)) + (1 + G3.*(G4 + G5) .*Z072) .*sin(x6.*(theta)))))))) + sin(x1.*( 
theta)) .*((-sin(x2.*(theta))).*(cos(x3.* (theta) ) .*(cos(x5.* (theta) ) .*( 
cos(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) + 1i.*(G4 + G5) 
.*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + sin(x5.*(theta)).*(cos(x4 
.*(theta)).*(cos(x6.*(theta)) + 11.*G5.*Z0.*sin(x6.*(theta))) - sin(x4 
.*(theta)) .*((-1i) .*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5.*Z072) .*sin 
(x6.*(theta))))) + 

sin(x3.*(theta)) .*(cos(x4.* (theta) ).*(cos(x5.* (theta) ) .*(cos(x6.* (theta) ) 
+ 11.*(G3 + G4 + G5).*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i 
).*(G3 + G4) .*Z0.*cos(x6.*(theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z072) 
.* sin(x6.*(theta)))) + sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + 
G3.*G4.*Z072) .*cos(x6.*(theta)) + 11.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) 
.*sin(x6.*(theta))) - cos(x5.*(theta)).*((-1i) .*G3.*Z0.*cos(x6.* (theta 
)) + (1 + G3.*(G4 + G5).*Z0^2) .*sin(x6.*(theta)))))) + cos(x2.*(theta) 
).*((-sin(x3.*(theta))).*(cos(x5.*(theta)).*(cos(x6.*(theta)).*sin(x4 
.*(theta)) + (cos(x4.*(theta)) + 11.*(G4 + G5).*Z0.*sin(x4.*(theta))) 
.*sin(x6.*(theta))) + sin(xb.*(theta)).*(cos(x4.*(theta)) .*(cos(x6.*( 
theta)) + 11.*G5.*Z0.*sin(x6.*(theta))) - sin(x4.*(theta)).*((-1i).*G4 
.*Z0.*cos(x6.*(theta)) + (1 + G4.+*G5.*Z0°2).*sin(x6.*(theta))))) + cos 
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G5) 


(x3.*(theta)).*(cos(x4.*(theta)).*(cos(x5.*(theta)).*(cos(x6.*(theta)) 
+ 1i.*(G3 + G4 + 
.*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4).*Z0.*cos 
(x6.*(theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z0°2).* sin(x6.*(theta)))) 
+ sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + G3.*G4.*Z072) .*cos(x6 
.*(theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z0^2) .*sin(x6.*(theta))) - 
cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + 
G5).*Z0^2) .*sin(x6.*(theta))))) + 1i.*G2.*Z0.*(cos(x3.*(theta)).*(cos( 
xb.*(theta)).*(cos(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) 
+ 11.*(G4 + G5).*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + sin(x5.*( 
theta)).*(cos(x4.*(theta)).*(cos(x6.*(theta)) + 1i.*G5.*Z0.*sin(x6.*( 
theta))) - sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + (1 + 
G4.*G5.*Z072).*sin(x6.*(theta))))) + sin(x3.*(theta)) .*(cos(x4.*(theta 
)).*(cos(x5.*(theta)).*(cos(x6.*(theta)) + 11.*(G3 + G4 + G5).*Z0.*sin 
(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4).*ZO.*cos(x6.*( 
theta)) + (1 + 


G3.*G5.*Z0°2 + G4.*G5.*Z0^2).*sin(x6.*(theta)))) + sin(x4.*(theta)).*((- 


sin(x5.*(theta))).*((1 + G3.*G4.*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 
+ G5 + G3.*G4.*G5.*Z0°2) .*sin(x6.*(theta))) - cos(x5.*(theta)).*((-1i 
) .*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5) .*Z072) .*sin(x6.*( 
theta))))))) + 11.*G1.*Z0.*(cos(x2.* (theta) ) .*(cos(x3.* (theta) ) .*(cos( 
x5.* (theta) ).*(cos(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) 
+ 11.*(G4 + G5).*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + sin(x5.*( 
theta) ) .*(cos(x4.* (theta) ).*(cos(x6.*(theta)) + 1i.*G5.*Z0.*sin(x6.*( 
theta))) - sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + (1 + 
G4.*G5.*Z0^2) .*sin(x6.*(theta))))) + sin(x3.*(theta)).*(cos(x4.* (theta 
)).*(cos(x5.*(theta)).*(cos(x6.*(theta)) + 11.*(G3 + G4 + G5).*Z0.*sin 
(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G3 + G4) .*Z0.*cos (x6. *( 
theta)) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z0"2).*sin(x6.*(theta)))) + sin( 
x4.*(theta)).*((-sin(x5.*(theta))).*((1 + 


G3.*G4.*20^2) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) .* 


sin(x6.*(theta))) - cos(x5.*(theta)).*((-1i) .*G3.*Z0.*cos(x6.* (theta) ) 
+ (1 + G3.*(G4 + G5).*Z072) .*sin(x6.*(theta)))))) + sin(x2.*(theta)) 
.*((-sin(x3.* (theta) )).*(cos(x5.* (theta) ).*(cos(x6.*(theta)) .*sin(x4 
.*(theta)) + (cos(x4.*(theta)) + 11.*(G4 + G5).*Z0.*sin(x4.* (theta) )) 
.*sin(x6.*(theta))) + sin(xb.*(theta)).*(cos(x4.*(theta)) .*(cos(x6.*( 
theta)) + 11.*G5.*Z0.*sin(x6.*(theta))) - sin(x4.*(theta)).*((-1i).*G4 
.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5.*Z0^2) .*sin(x6.*(theta))))) + cos 
(x3.*(theta)).*(cos(x4.*(theta)).*(cos(x5.*(theta)).*(cos(x6.*(theta)) 
+ 11.*(G3 + G4 + G5).*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i 
).*(G3 + G4).*Z0.*cos(x6.*(theta)) + (1 + G3.*G5.*Z0^2 + G4.*G5.*Z072) 
.*sin(x6.*(theta)))) + sin(x4.*(theta)).*((-sin(x5.*(theta))).*((1 + 
G3.*G4.*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) 
.*sin(x6.*(theta))) - 
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cos(x5.*(theta)).*((-1i).*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + G5) .* 
20^2).*sin(x6.*(theta))))) + 1i.*G2.*Z0.*(cos(x3.*(theta)).*(cos(x5.*( 
theta) ) .*(cos(x6.*(theta)).*sin(x4.*(theta)) + (cos(x4.*(theta)) + 1i 
.*(G4 + G5).*Z0.*sin(x4.*(theta))).*sin(x6.*(theta))) + sin(x5.*(theta 
)).* (cos (x4.* (theta) ).*(cos(x6.*(theta)) + 11.*G5.*Z0.*sin(x6.* (theta 
))) - sin(x4.*(theta)).*((-1i).*G4.*Z0.*cos(x6.*(theta)) + (1 + G4.*G5 
.*Z0^2).*sin(x6.*(theta))))) + sin(x3.*(theta)) .*(cos(x4.*(theta)) .*( 
cos(x5.*(theta)).*(cos(x6.*(theta)) + 1i.*(G3 + G4 + G5).*Z0.*sin(x6 
.*(theta))) - sin(xb.*(theta)).*((-1i).*(G3 + G4) .*Z0.*cos(x6.* (theta) 
) + (1 + G3.*G5.*Z0°2 + G4.*G5.*Z072) .*sin(x6.*(theta)))) + sin(x4.*( 
theta)) .*((-sin(x5.*(theta))).*((1 + G3.*G4.*Z0°2) .*cos(x6.*(theta)) + 
1i.*Z0.*(G3 + G5 + G3.*G4.*G5.*Z072) .* sin(x6.*(theta))) - cos(x5.*( 
theta) ) .*((-1i) .*G3.*Z0.*cos(x6.*(theta)) + (1 + G3.*(G4 + 
GB) .*Z0°2) .*sin(x6.*(theta)))))))))))./ (cos(x1.*(theta)) .*(cos (x2. zt 
theta) ) .*(cos(x3.*(theta)).*(cos(x4.* (theta) ).*(cos(x5.*(theta)).*(cos 
(x6.*(theta)) + 11.*G6.*Z0.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1 
i).*(G5 + G6).*Z0.*cos(x6.*(theta)) + (1 + G5.*G6.*Z072) .* sin(x6.*( 
theta)))) + 1i.*sin(x4.*(theta)).*(cos(x5.*(theta)).*((G4 + G5 + G6).* 
ZO.*cos(x6.*(theta)) + 1i.*(1 + G4.*G6.*Z0°2 + G5.*G6.*Z072) .*sin(x6 
.*(theta))) + 1i.*sin(x5.*(theta)).*((1 + G4.*(G5 + G6).*Z0^2) .*cos(x6 
.*(theta)) + 1i.*Z0.*(G4 + G6 + G4.*G5.*G6.+*Z0°2) .*sin(x6.*(theta))))) 
+ 1i.*sin(x3.*(theta)) .*(cos(x4.* (theta) ) .*(cos(x5.*(theta)).*((G3 + 
G4 + G5 + G6).*Z0.*cos(x6.*(theta)) + 1i.*(1 + G3.*G6.*Z20^2 + G4.*G6.x* 
ZO°2 + G5.*G6.*Z072) .*sin(x6.*(theta))) + 1i.*sin(x5.*(theta)).*((1 + 
G3.*(G5 + G6) .*Z072 + G4.*(G5 + G6) .*Z072) .*cos(x6.*(theta)) + 1i.*Z0 
.*(G3 + G4 + G6 + 
G3.*G5.*G6.*Z0°2 + G4.*G5.*G6.*Z0°2) .*sin(x6.*(theta)))) + 1i.*sin(x4.*( 
theta)).*(cos(x5.*(theta)).*((1 + G3.*(G4 + G5 + G6) .*Z072) .*cos(x6.*( 
theta)) + 1i.*Z0.*(G3 + G6 + G3.*G4.*G6.*Z072 + G3.*G5.*G6.*Z072) .*sin 
(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*Z0.*(G5 + G6 + G3.*(1 + G4 
.*(G5 + G6) .*Z0°2)).*cos(x6.*(theta)) + (1 + G5.*G6.*Z0^2 + G3.*Z0 
72.*(G4 + G6 + G4.*G5.*G6.*Z0°2)) .*sin(x6.*(theta)))))) + 1i.*sin(x2 
.* (theta) ).*(cos(x3.* (theta) ) .*(cos(x4.* (theta) ).*(cos(x5.* (theta) ) 
.*((G2 + G3 + G4 + GB + G6).*Z0.*cos(x6.*(theta)) + 1i.*(1 + G2.*G6.* 
ZO°2 + G3.*G6.*Z0"2 + G4.*G6.*Z0°2 + G5.*G6.*Z0°2).* sin(x6.*(theta))) 
+ 1i.*sin(x5.*(theta)).*((1 + G4.*G5.*Z0°2 + G4.*G6.*Z0°2 + G2.*(G5 + 
G6) .*Z072 + G3.*(G5 + G6) .*Z0°2) .*cos(x6.*(theta)) + 1i.*Z0.*(G2 + G3 
+ G4 + G6 + G2.*G5.*G6.*Z0°2 + G3.*G5.*G6.*Z072 + G4.*G5.*G6.*Z072) .* 
sin(x6.*(theta)))) + 
1i.*sin(x4.*(theta)).*(cos(x5.*(theta)).*((1 + G2.*(G4 + G5 + G6).*Z0°2 + 
G3.*(G4 + G5 + G6).*Z0°2).* cos(x6.*(theta)) + 1i.*Z0.*(G2 + G3 + Gë 
+ G2.*G4.*G6.*Z0^2 + G3.*G4.*G6.*Z0°2 + G2.*G5.*G6.*Z0^2 + G3.*G5.*G6 
.*Z0^2).*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*Z0.*(G3 + G5 + 
G6 + G3.*G4.*G5.*Z0°2 + G3.*G4.*G6.*Z0"2 + G2.*(1 + G4.*(G5 + G6).*ZO 
72)).*cos(x6.*(theta)) + (1 + G5.*G6.*Z072 + G2.*Z0°2.*(G4 + G6 + G4.x* 
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G5.*G6.*Z0°2) + G3.*Z0°2.*(G4 + G6 + G4.*G5.*G6.*Z072)) .*sin(x6.*( 
theta))))) - sin(x3.*(theta)).*(cos(x4.*(theta)) .*((-1i) .*cos (x5. *( 
theta)).*((1 + G2.*(G3 + G4 + G5 + G6) .*Z072) .*cos(x6.*(theta)) + li.* 
Z0.*(G6 + G2.*(1 + G3.*G6.*Z0"2 + G4.*G6.*Z0°2 + G5.*G6.*Z20^2)).* Sint 
x6.*(theta))) + sin(x5.*(theta)).*(Z0.*(G5 + G6 + G2.*(1 + G3.*(G5 + 
G6) .*Z0°2 + G4.*(G5 + G6) .*Z0°2)).*cos(x6.*(theta)) + 1i.*(1 + G5.*G6 
.*Z0^2 + G2.*Z0°2.*(G3 + G4 + G6 + 

G3.*G5.*G6.*Z0°2 + G4.*G5.*G6.*Z0°2)).*sin(x6.*(theta)))) + sin(x4.*( 
theta)).*(cos(x5.*(theta)).*(Z0.*(G4 + G5 + G6 + G2.*(1 + G3.*(G4 + G5 
+ G6).*Z0°2)).* cos(x6.*(theta)) + 1i.*(1 + G4.*G6.*Z0"2 + G5.*G6.*Z0 
^2 + G2.*Z0°2.*(G3 + G6 + G3.*G4.*G6.*Z20^2 + G3.*G5.*G6.*Z072)) .*sin( 
x6.*(theta))) + 1i.*sin(x5.*(theta)).*((1 + G4.*(G5 + G6).*Z0^2 + G2.* 
Z0°*2.*(G5 + G6 + G3.*(1 + G4.*(G5 + G6).*Z072))).*cos(x6.*(theta)) + 1 
i.*Z0.*(G4 + G6 + G4.*G5.*G6.*Z0"2 + G2.*(1 + G5.*G6.*Z072 + G3.*Z0 
72.*(G4 + G6 + G4.*G5.*G6.*Z0°2))).*sin(x6.*(theta))))))) + 1i.x*Z0.x* 
sin(x1.*(theta)).*((1./Z0) .*(sin(x2.* (theta) ).*(cos(x3.*(theta)) .*(1i 
.*cos(x4.* (theta) ) .*(cos(x5.* (theta) ).*(cos(x6.*(theta)) + 1i.*G6.*Z0 
.*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*(G5 + G6) .*Z0.*cos(x6 
.*(theta)) + (1 + G5.*G6.*Z0°2) .*sin(x6.*(theta)))) + sin(x4.*(theta) ) 
.*((-1i).*cos(x5.*(theta)).*((-1i).*(G4 + Gb + 

G6) .*Z0.*cos(x6.*(theta)) + (1 + G4.*G6.*Z0°2 + G5.*G6.*Z072) .*sin(x6.*( 
theta))) + sin(x5.*(theta)).*((-1i).*(1 + G4.*(G5 + G6) .*Z0°2) .*cos (x6 
.*(theta)) + Z0.*(G4 + G6 + G4.+*G5.*G6.*Z0°2) .*sin(x6.*(theta))))) + 
sin(x3.*(theta)).*((-1i).*cos(x4.*(theta)).*(cos(x5.*(theta)).*((-1i) 
.*(G3 + G4 + G5 + G6).*Z0.*cos(x6.*(theta)) + (1 + G3.*G6.*Z0"2 + G4.* 
G6.*Z072 + G5.*G6.*Z072).*sin(x6.*(theta))) + sin(x5.*(theta)).*((1 + 
G3.*(G5 + G6) .*Z072 + G4.*(G5 + G6) .*Z072) .*cos(x6.*(theta)) + 1i.*Z0 
.*(G3 + G4 + G6 + G3.*G5.*G6.*Z072 + G4.*G5.*G6.*Z072) .*sin(x6.*(theta 
)))) + sin(x4.*(theta)).*(cos(x5.*(theta)).*((-1i).*(1 + G3.*(G4 + G5 
+ G6).*Z072).*cos(x6.*(theta)) + ZO.*(G3 + G6 + G3.*G4.*G6.*Z072 + G3 
.*G5.*G6.*Z0°2) .*sin(x6.*(theta))) + sin(xb.*(theta)).*(Z0.*(G5 + G6 + 
G3.*(1 + G4.*(G5 + G6) .*Z0°2)).*cos(x6.*(theta)) + 1i.*(1 + Gb.*G6.x* 
Z0"2 + G3.*Z0°2.*(G4 + G6 + 

G4.*G5.*G6.*Z0°2)) .*sin(x6.*(theta))))))) + (1./Z0) .*(cos(x2.* (theta) ) .*( 
cos (x3.* (theta) ) .*(cos(x4.* (theta) ) .*(cos(x5.*(theta)).*((G2 + G3 + G4 
+ G5 + G6).*Z0.*cos(x6.*(theta)) + 1i.*(1 + G2.*G6.*Z0°2 + G3.*G6.*Z0 
72 + G4.*G6.*Z0^2 + G5.*G6.*Z20^2).* sin(x6.*(theta))) + 1i.*sin(x5.*( 
theta)).*((1 + G4.*G5.*Z072 + G4.*G6.*Z0°2 + G2.*(G5 + G6).*Z0°2 + G3 
.*(G5 + G6) .*Z072) .*cos(x6.*(theta)) + 1i.*Z0.*(G2 + G3 + G4 + G6 + G2 
.*G5.*G6.*Z0°2 + G3.*G5.*G6.*Z0°2 + G4.*G5.*G6.*Z072) .*sin(x6.* (theta) 
))) + 1i.*sin(x4.* (theta) ).*(cos(x5.*(theta)).*((1 + G2.*(G4 + G5 + Gë 
).*Z0^2 + G3.*(G4 + G5 + G6).*Z072).*cos(x6.*(theta)) + 1i.*Z0.*(G2 + 
G3 + G6 + G2.*G4.*G6.*Z0"2 + G3.*G4.*G6.*Z0°2 + G2.*G5.*G6.*Z072 + G3 
.*G5.*G6.*Z072) .*sin(x6.*(theta))) - sin(x5.*(theta)).*((-1i).*Z0.*(G3 
+ G5 + G6 + G3.*G4.*G5.*Z0°2 + G3.*G4.*G6.*Z072 + G2.*(1 + G4.*(G5 + 
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G6) .*Z072)).* cos(x6.*(theta)) + (1 + 

G5.*G6.*Z0°2 + G2.*Z0°2.*(G4 + G6 + G4.*G5.*G6.*Z072) + G3.*Z072.*(G4 + 
G6 + G4.*G5.*G6.*Z072)).*sin(x6.*(theta))))) - sin(x3.*(theta)) .*(cos( 
x4.*(theta)).*((-1i).*cos(xb.*(theta)).*((1 + G2.*(G3 + G4 + G5 + G6) 
.*Z0^2).*cos(x6.*(theta)) + 1i.*Z0.*(G6 + G2.*(1 + G3.*G6.*Z0"2 + G4.x* 
G6.*Z072 + G5.*G6.*Z072)).* sin(x6.*(theta))) + sin(xb.*(theta)).*(ZO 
.*(G5 + G6 + G2.*(1 + G3.*(G5 + G6).*Z0°2 + G4.*(G5 + G6).*Z072)).*cos 
(x6.*(theta)) + 1i.*(1 + G5.*G6.*Z072 + G2.*Z0°2.*(G3 + G4 + G6 + G3.* 
G5.*G6.*Z072 + G4.*G5.*G6.*Z0°2)) .*sin(x6.*(theta)))) + sin(x4.*(theta 
)) .* (cos (x5.*(theta)).*(Z0.*(G4 + Gb + G6 + G2.*(1 + G3.*(G4 + G5 + G6 
).*Z072)).* cos(x6.*(theta)) + li.*(1 + G4.*G6.*Z0°2 + G5.*G6.*Z072 + 
G2.*Z0°2.*(G3 + G6 + G3.*G4.+*G6.*Z0°2 + G3.*G5.*G6.*Z072)) .*sin(x6.*( 
theta))) + 1i.*sin(x5.*(theta)).*((1 + G4.*(G5 + G6).*Z0°2 + 

G2.*Z0°2.*(G5 + G6 + G3.*(1 + G4.*(G5 + G6).*Z0°2))).*cos(x6.*(theta)) + 
1i.*Z0.*(G4 + G6 + G4.*G5.*G6.*Z0"2 + G2.*(1 + G5.*G6.*Z0^2 + G3.*Z0 
72.*(G4 + G6 + G4.*G5.*G6.*Z0°2))).*sin(x6.*(theta))))))) + G1.*(cos( 
x2.* (theta) ).*(cos(x3.* (theta) ) .*(cos(x4.* (theta) ).*(cos(x5.* (theta) ) 
.*(cos(x6.*(theta)) + 1i.*G6.*Z0.*sin(x6.*(theta))) - sin(x5.*(theta) ) 
.*((-1i).*(G5 + G6) .*Z0.*cos(x6.*(theta)) + (1 + G5.*G6.*Z072) .* sin( 
x6.*(theta)))) + 1i.*sin(x4.*(theta)).*(cos(x5.*(theta)).*((G4 + G5 + 
G6) .*Z0.*cos(x6.*(theta)) + li.*(1 + G4.*G6.*Z0°2 + G5.*G6.*20^2).*sin 
(x6.*(theta))) + 1i.*sin(x5.*(theta)).* ((1 + G4.*(G5 + G6) .*Z072) .* 
cos(x6.*(theta)) + 11.*Z0.*(G4 + G6 + G4.+*G5.*G6.*Z0°2) .*sin(x6. zt 
theta))))) + 1i.*sin(x3.*(theta)).*(cos(x4.* (theta) ).*(cos(x5.* (theta) 
).*((G3 + G4 + G5 + G6).*Z0.*cos(x6.*(theta)) + 

1i.*(1 + G3.*G6.*Z0"2 + G4.*G6.*Z0"2 + G5.*G6.*Z0°2) .*sin(x6.*(theta))) + 
1i.*sin(x5.*(theta)).* ((1 + G3.*(G5 + G6).*Z0°2 + G4.*(G5 + G6).*Z0 
72) .*cos(x6.*(theta)) + 11.*Z0.*(G3 + G4 + G6 + G3.*G5.*G6.*Z0"2 + G4 
.*G5.*G6.*Z072) .*sin(x6.*(theta)))) + 1i.*sin(x4.*(theta)).*(cos(x5.*( 
theta)).*((1 + G3.*(G4 + G5 + G6) .*Z072) .*cos(x6.*(theta)) + 11.*Z0.*( 
G3 + G6 + G3.*G4.*G6.*Z072 + G3.*G5.*G6.*Z20^2) .*sin(x6.*(theta))) - 
sin(x5.*(theta)).* ((-1i).*Z0.*(G5 + G6 + G3.*(1 + G4.*(G5 + G6) .*Z0 
72)).*cos(x6.*(theta)) + (1 + G5.*G6.*Z072 + G3.*Z0°2.*(G4 + G6 + G4.* 
G5.*G6.*Z072)) .*sin(x6.*(theta)))))) + 1i.*sin(x2.* (theta) ) .*(cos(x3 
.* (theta) ).*(cos(x4.* (theta) ) .*(cos(x5.*(theta)).*((G2 + G3 + G4 + G5 
+ G6) .*Z0.*cos(x6.*(theta)) + 1i.*(1 + G2.*G6.*Z0°2 + G3.*G6.*20^2 + 
G4.*G6.*Z0°2 + G5.*G6.*Z0°2).* sin(x6.*(theta))) + 11.*sin(x5.* (theta) 
).*((1 + G4.*G5.*Z0"2 + G4.*G6.*Z072 + 

G2.*(G5 + G6) .*Z072 + G3.*(G5 + G6).*Z0°2).*cos(x6.*(theta)) + 1i.*ZO.*( 
G2 + G3 + G4 + G6 + G2.*G5.*G6.*Z0°2 + G3.*G5.*G6.*Z072 + G4.*G5.*G6.* 
20^2).*sin(x6.*(theta)))) + 1i.*sin(x4.*(theta)).*(cos(x5.* (theta) ) 
.*((1 + G2.*(G4 + G5 + G6).*Z0°2 + G3.*(G4 + G5 + G6) .*Z0°2) .*cos(x6 
.*(theta)) + 1i.*Z0.*(G2 + G3 + G6 + G2.*G4.*G6.*Z0°2 + G3.*G4.*G6.*Z0 
^2 + G2.*G5.*G6.*Z0°2 + G3.*G5.*G6.*Z0°2).*sin(x6.*(theta))) - sin(x5 
.*(theta)).*((-1i).*Z0.*(G3 + G5 + G6 + G3.*G4.*G5.*Z0°2 + G3.*G4.*G6 
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G5 


end 


.*Z0^2 + G2.*(1 + G4.*(G5 + G6).*Z0°2)).* cos(x6.*(theta)) + (1 + G5.* 
G6.*Z072 + G2.*Z0°2.*(G4 + G6 + G4.*G5.*G6.*Z072) + G3.*Z0^2.*(G4 + G6 
+ G4.*G5.*G6.*Z0°2)).*sin(x6.*(theta))))) - sin(x3.*(theta)).*(cos(x4 
.*(theta)).*((-1i).*cos(x5.*(theta)).*((1 + G2.*(G3 + G4 + G5 + G6).* 

20^2).*cos(x6.*(theta)) + 1i.*Z0.*(G6 + G2.*(1 + G3.*G6.*Z0^2 + G4.*G6 
.*Z0^2 + 


.*G6.*20^2)).* sin(x6.*(theta))) + sin(x5.*(theta)).*(Z0.*(G5 + G6 + G2 


.*(1 + G3.*(G5 + G6).*Z0^2 + G4.*(G5 + G6).*Z072)).*cos(x6.*(theta)) + 
1i.*(1 + G5.*G6.*Z0°2 + G2.*20^2.*(G3 + G4 + G6 + G3.*G5.*G6.*Z072 + 
G4.*G5.*G6.*Z072)) .*sin(x6.*(theta)))) + sin(x4.*(theta)).*(cos(xb.*( 
theta)).*(Z0.*(G4 + G5 + G6 + G2.*(1 + G3.*(G4 + G5 + G6) .*Z072)) .* 
cos(x6.*(theta)) + 1i.*(1 + G4.*G6.*Z0°2 + G5.*G6.*Z20^2 + G2.*Z0"2.x*( 
G3 + G6 + G3.*G4.*G6.*Z0°2 + G3.*G5.*G6.*Z0°2)) .*sin(x6.*(theta))) + 1 
i.*sin(xb.*(theta)).* ((1 + G4.*(G5 + G6).*Z0^2 + G2.*Z0°2.*(G5 + G6 + 
G3.*(1 + G4.*(G5 + G6) .*Z0°2))).*cos(x6.*(theta)) + 1i.*Z0.*(G4 + G6 
+ G4.*G5.*G6.*Z0°2 + G2.*(1 + G5.*G6.*Z072 + G3.*Z0"2.*(G4 + G6 + G4.x* 
G5.*G6.*Z072))) .*sin(x6.*(theta)))))))))))); 


record theta = zeros(1,length(theta)); % Initialization for record theta 
record theta(1,:) = ((gamma(1,:)) < dB level).*theta in Hz(1,:); % Take 


values of theta when gamma < -10 dB 


flag = 0; % Flag for helping calculate bandwidth 
bandwidth low end = 0; % Initialization 
bandwidth high end = 0; % Initialization 


^ Note: only records the first bandwidth that goes below -20 dB 


for 


end 


n = 1:length(theta) 


if ((record theta(1,n) “= 0) && (flag == 0)) 
bandwidth low end = theta in Hz(1,n); 
flag = 1; 

end 


if ((record_theta(1,n) == 0) && (flag == 1)) 
bandwidth high end = theta in Hz(1,n); 
flag = 2; 

end 


if(flag == 1) % If the bandwidth high end isn't detected due to the 


crossing point too close to the end, then set the high end at the end 
of the bandwidth 
bandwidth_high_end = theta_in_Hz(1,n); 
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end 
bandwidth = bandwidth_high_end - bandwidth_low_end; 
fitness = bandwidth/f0; 


end 
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A.6 MATLAB code for plotting convenience 


function S11FitnessGraph_with_var_L( theta in Hz, gamma ) 
% An add-on function to Si1Fitness with var L 
% Feed it theta in Hz and gamma and it outputs a graph of $11 


dB. level = -20; 


record theta = zeros(1,length(theta in Hz)); Z Initialization for 
record theta 

record theta(1,:) = ((gamma(1,:)) < dB level).*theta in Hz(1,:); % Take 
values of theta when gamma « dB level 

flag = 0; % Flag for helping calculate bandwidth 

bandwidth low end = 0; % Initialization 

bandwidth, high end = 0; % Initialization 


^ Note: only records the first bandwidth that goes below the dB level 
for n = 1:length(theta in Hz) 


if ((record theta(1,n) “= 0) && (flag == 0)) 
bandwidth low end = theta, in Hz(1,n); 
flag = 1; 

end 


if ((record_theta(1,n) == 0) && (flag == 1)) 
bandwidth_high_end = theta_in_Hz(1,n); 
flag = 2; 

end 


end 


if(flag == 1) % If the bandwidth high end isn’t detected due to the 
crossing point too close to the end, then set the high end at the end 
of the bandwidth 
bandwidth_high_end = theta_in_Hz(1,n); 

end 


bandwidth = bandwidth_high_end - bandwidth_low_end; 

^ display text = sprintf(’Bandwidth: %f - %f’ ,bandwidth_low_end, 
bandwidth high end); 

^ disp(display text); 


fO - 2.50e2; 


k low = bandwidth low. end/f0; 
k high = bandwidth high end/f0; 
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display text = sprintf ("Relative bandwidth: %f fO - %f f0”,k_low,k_high); 
disp(display_text) ; 
disp(?Factors:?); 


display text = sprintf (Af Nt 4f^,k low,k high); 
disp(display text) 


% Line markers 
lineMarkers [?=k ?3?:k ?;?--k ?;?-.K ?;?-sk ?;”:sk ’;’?--sk’;’-.sk’]; 
lineMarkers = cellstr(lineMarkers); 


figure 

hold on 

plot(theta in Hz,gamma,char(lineMarkers(1))); 
axis([min(theta in Hz) max(theta in Hz) -50 01) 
grid on 

title(?S (11) vs Frequency?) 

xlabel(’Frequency (MHz)?) 

ylabel(”S (11) (dB)’) 

end 
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A.7 General form of reflection coefficients with variable transmission 
line length for 1, 2, 3, 4, 5, and 6 layers case 


The following reflection coefficients are obtained through the same method described 
in chapter 2 except that each transmission line sections are treated as a unique 
variable. In other words, the transmission line length in the first, second, third, 


fourth, fifth, and sixth section are 01, 02, 05, 04, 05, and Og respectively. 


PM (Zs cos (110) + jZo(—1 + G1Zs) sin (210)) 
Y= (Zs cos (210) + jZo(1 + G1Zs) sin (x10)) (A.17) 








T —  (cos(x10)(Zg cos (z20)4-j Zo(—14-G2 Zs) sin (120))+sin (710) (jZ9(—14+G1 Zs +G2Zs) cos (220) -(Z5+G1Z3(—1+G2Z3)) sin (20))) 
2 ((— cos (z10))(Zs cos (1220) +j Zo (1+-G2Zs) sin (x20))+sin (210) (3 Zo (1--G1 Zs +G2Zg) cos (120) +(Z3+G1Z3 (14G2Z3)) sin (z20))) 


(A.18) 





























Ls = 

(A.19) 
Ez 

(A.20) 
Ds = 

(A.21) 


(Place holder for six layers case. Error because the expression is too long.) 
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A.8 Cut-and-try Algorithm Results 
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FIGURE A.1: Cut-and-try result for 2 sections case 
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FIGURE A.2: Cut-and-try result for 3 sections case 
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Sa vs Frequency 


S,, (dB) 
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FIGURE A.3: Cut-and-try result for 4 sections case 


Si vs Frequency 


S,, (dB) 
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FIGURE A.4: Cut-and-try result for 5 sections case 


A.9 Genetic Algorithm Expanded Results Tables 


Note: the standard result is equivalent to expanded result with theta limit set to 1 

This chapter lists all the GA expanded result. Column ”G„” is the normalized 
conductance and " X,,” is the separation between the element n and (n — 1) expressed 
in the multiples of 1/4 (the full detail is available in Section 2.3). The Oth element 
is always a short circuit. The ”RBW”, "Start", and "End" columns indicate the 


relative bandwidth, the start, and the end of the bandwidth in the multiples of fo 
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respectively. For example, if the bandwidth spans from 0.5 fo to 1.7 fo, then the entry 
is listed as "1.2" in the "RBW" column, "0.5" in the "Start" column, and ”1.7” in 
the "End" column. Finally, the "T" columns stands for "Theta limit”. The theta 
limit defines the search boundary in frequency in terms of fo between 0 to 2T fy. For 
example, if the theta limit is set to "77, then the algorithm searches for the best 
solution—considering other limitations as described in Section 3.4)—between 0 to 


14 fo. 
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A.9.1 Threshold level: -10 dB 


One section 


Table A.1: GA expanded result for 1 section with I" threshold -10 dB. 

































































Gi Xi RBW | Start End | T 
1.210 | 0.695 | 1.122 | 0.878 | 2.000 1 
1.220 | 0.347 | 2.244 1.756 | 4.000 | 2 
1.232 | 0.232 | 3.364 | 2.632 | 5.995 | 3 
1.222 | 0.174 | 4.486 | 3.509 | 7.995 | 4 
1.224 | 0.139 | 5.608 | 4.386 | 9.994 | 5 
1.224 | 0.116 | 6.730 | 5.271 | 12.000 | 6 
1.221 | 0.099 | 7.853 | 6.147 | 14.000 | 7 
1.216 | 0.087 | 8.975 | 7.025 | 16.000 | 8 
1.220 | 0.077 | 10.097 | 7.903 | 18.000 | 9 
1.220 | 0.070 | 11.220 | 8.780 | 20.000 | 10 
1.225 | 0.063 | 12.342 | 9.654 | 21.996 | 11 
1.223 | 0.058 | 13.465 | 10.536 | 24.000 | 12 
1.223 | 0.054 | 14.586 | 11.409 | 25.995 | 13 
1.225 | 0.050 | 15.707 | 12.293 | 28.000 | 14 
1.226 | 0.046 | 16.811 | 13.189 | 30.000 | 15 
1.222 | 0.044 | 17.939 | 14.031 | 31.970 | 16 
1.225 | 0.041 | 19.072 | 14.929 | 34.000 | 17 
1.219 | 0.039 | 20.180 | 15.821 | 36.000 | 18 
1.223 | 0.037 | 21.312 | 16.688 | 38.000 | 19 
1.219 | 0.035 | 22.400 | 17.521 | 39.922 | 20 
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Two sections 


Table A.2: GA expanded result for 2 sections with I' threshold -10 dB. 

































































Gy Go Xi X» RBW Start End T 
1.274 | 0.360 | 0.843 | 1.654 | 1.646 | 0.354 | 2.000 | 1 
0.979 | 0.901 | 0.419 | 0.419 | 3.229 | 0.771 | 4.000 | 2 
0.943 | 0.860 | 0.282 | 0.278 | 4.844 | 1.154 | 5.997 | 3 
1.674 | 0.657 | 0.426 | 0.216 | 6.478 | 1.400 | 7.878 | 4 
1.543 | 0.657 | 0.354 | 0.175 | 8.251 | 1.560 | 9.812 | 5 
1.569 | 0.599 | 0.256 | 0.130 | 9.886 | 2.114 | 12.000 | 6 
1.589 | 0.577 | 0.211 | 0.106 | 11.400 | 2.600 | 14.000 | 7 
1.364 | 0.389 | 0.105 | 0.220 | 13.337 | 2.643 | 15.980 | 8 
1.370 | 0.465 | 0.093 | 0.186 | 15.180 | 2.820 | 18.000 | 9 
0.960 | 0.879 | 0.082 | 0.084 | 16.134 | 3.866 | 20.000 | 10 
0.941 | 0.825 | 0.078 | 0.075 | 17.749 | 4.249 | 21.997 | 11 
1.573 | 0.639 | 0.145 | 0.073 | 19.902 | 3.827 | 23.728 | 12 
1.541 | 0.676 | 0.121 | 0.060 | 21.415 | 4.585 | 26.000 | 13 
1.656 | 0.664 | 0.114 | 0.056 | 22.666 | 5.334 | 28.000 | 14 
0.984 | 0.860 | 0.058 | 0.054 | 24.180 | 5.820 | 30.000 | 15 
1.028 | 0.929 | 0.052 | 0.053 | 25.789 | 6.212 | 32.000 | 16 
1.270 | 0.392 | 0.051 | 0.101 | 28.153 | 5.577 | 33.730 | 17 
0.939 | 0.851 | 0.046 | 0.047 | 29.058 | 6.942 | 36.000 | 18 
0.969 | 0.861 | 0.045 | 0.044 | 30.648 | 7.325 | 37.973 | 19 
1.308 | 0.396 | 0.041 | 0.079 | 32.998 | 7.003 | 40.000 | 20 
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Three sections 


Table A.3: GA expanded result for 3 sections with T threshold -10 dB. 































































































G4 Go G3 Xi X» Xa RBW Start End T 
1.600 | 0.524 | 0.769 | 0.928 | 0.916 | 1.796 | 1.817 | 0.183 | 2.000 | 1 
1.487 | 0.497 | 0.699 | 0.431 | 0.445 | 0.933 | 3.628 | 0.373 | 4.000 | 2 
1.411 | 0.568 | 0.644 | 0.286 | 0.277 | 0.533 | 5.373 | 0.627 | 6.000 | 3 
1.159 | 0.628 | 0.537 | 0.208 | 0.198 | 0.618 | 7.281 | 0.719 | 8.000 | 4 
1.256 | 0.782 | 0.611 | 0.336 | 0.387 | 0.166 | 9.175 | 0.825 | 10.000 | 5 
1.191 | 0.611 | 0.470 | 0.142 | 0.159 | 0.468 | 10.985 | 1.016 | 12.000 | 6 
1.184 | 0.656 | 0.570 | 0.122 | 0.122 | 0.374 | 12.834 | 1.167 | 14.000 | 7 
1.442 | 0.660 | 0.490 | 0.106 | 0.108 | 0.328 | 14.590 | 1.410 | 16.000 | 8 
1.314 | 0.609 | 0.522 | 0.097 | 0.096 | 0.291 | 16.460 | 1.541 | 18.000 | 9 
1.236 | 0.628 | 0.553 | 0.089 | 0.096 | 0.280 | 18.388 | 1.568 | 19.956 | 10 
1.243 | 0.707 | 0.543 | 0.077 | 0.083 | 0.237 | 20.160 | 1.840 | 22.000 | 11 
1.088 | 0.616 | 0.521 | 0.068 | 0.080 | 0.225 | 21.994 | 2.006 | 24.000 | 12 
1.188 | 0.776 | 0.507 | 0.067 | 0.062 | 0.197 | 23.696 | 2.304 | 26.000 | 13 
1.283 | 0.606 | 0.381 | 0.057 | 0.179 | 0.121 | 25.947 | 2.053 | 28.000 | 14 
1.169 | 0.604 | 0.423 | 0.051 | 0.144 | 0.100 | 27.574 | 2.427 | 30.000 | 15 
1.194 | 0.709 | 0.533 | 0.055 | 0.056 | 0.163 | 29.322 | 2.679 | 32.000 | 16 
1.233 | 0.668 | 0.580 | 0.055 | 0.052 | 0.160 | 31.317 | 2.684 | 34.000 | 17 
1.225 | 0.683 | 0.545 | 0.046 | 0.043 | 0.142 | 32.839 | 3.162 | 36.000 | 18 
1.230 | 0.681 | 0.566 | 0.044 | 0.042 | 0.128 | 34.632 | 3.369 | 38.000 | 19 
1.363 | 0.425 | 0.619 | 0.044 | 0.038 | 0.099 | 36.147 | 3.853 | 40.000 | 20 
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Four sections 
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OG | 00070» | I9Z | 682'9£ | 6600 | eut | 4800 | 600 | IPO | 140 | 9660 | PSOE 
GI | 000'8£ | ELEC | LZ9°ZE | €80'0 | 0700 | 7EO'O | OTTO | FSO | G£90 | PIG'E | PEST 
ST | 000:98 | 6979 | TPG'6Z | 8800 | YF0'O | 1080 | 6000 | STFO | 060T | 9879 | TILO 
LT | 000'P£ | GEO | SP6'0€ | EPTO | 000 | 9700 08400 | PIFO | 290 | OCC T | OFF 
OT | 000'c8 | 7E0'9 | 446 `9% | £900 | 6F0'0 | 60010 | 67£'0 | 9780 | C980 | €€0`¿ | T69T 
ST | 000'08 | £F T | 27987 | 96T0 | 0920 | 8900 | 62010 | £960 | 69F0 | EFRO | GCTT 
PI | 000'8G | ETT | 998°Sz | 0900 | 0900 | 9810 | $400 | 4670 | 970 | SZOT | CZ26T 
ET | 000796 | TLE | 6cc c6 | 8900 | ceTO | 88010 | 6070 | PISO | OICT | 2199 | 29918 
ZL | 00076 | 9071 | 96/6 | 8910 | 0200 | PIEO | 000 | CT£€'O | 4780 | 1360 | 6SF'T 
TT | 000.37 | TOE'T | 66903 | EET'O | 400 9080 | 2800 | 62€'0 | 9970 | OITT | SLET 
OT | 000:0 | £80'T | ST6'8T | 0600 | 6200 | ETEO | 9⁄0 | 2880 | ¿870 | 9280 | 2871 
6 | 000'8T | £80'T | ¿T6'9T | 06T'0 | 0600 | TOSO | I800 | 2070 | BTO | 070 | ZZ0'T 
8 | 000'9T | £480 | 9ZT'ST | 9870 | 0370 | P9€'0 | 8800 | Z8TO | PLEO | 890 | GIGTT 
L [00071 | 8741 | CAC CT | 0800 | PITO | 6810 | OPT | 9FS'0 | 980 | 206% | 860% 
9 | 000'GT | 2890 | ETE IT | 8870 | STO | 3670 | OTTO | 0280 | 8970 | 962'0 | 9890 
G | 2166 | 68€ T | Z898 | ELO | FAT'O | C0 | 69€0 | 6870 | ZCO'T | ZE6'G | C66 C 
P | 98674 | TOO | ESE | 8910 0670 | 960 | POSO | 2600 | 2370 | 0820 | GET 
€ | 0009 | 3070 | 6699 | 4680 | PESO | 0270 | 080 | E880 | 2980 | £070 | FI9'T 
6 | 0007 | €970 | rie 2980 | 1660 | 9680 | 6190 | 29€0 | 81910 | TRIO | LEET 
I | 0007 | SFT O | S&T | 0980 | 9960 | IIST | £960 | 6290 | T870 | C480 | SKOTI 
L| Pua | wes | MAH | ÓX X X CR TÐ Ð Sei !ð 








"OD OT- PIOYS9IUI J YIM SUOIDAS p 10} 3[nso1 popuedxo VI :p y 


QLL 
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Five sections 
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OG | 00070» | 966'p | POO'SE | OPO'O | 8800 SS0'0 | 9490 | 6110 | 9770 | C920 | TAPT | VOET | 6590 
61 | 00088 | SPE T | 299 € | 9700 | 0600 | ZPO'0 | LIZO | AF€'O | EZ70 | 6790 | cS T | COV G | 416€ 
ST | 000'98 TPT | 68C T£ | 6600 | £500 | T£00 | ZOTO | 2220 | SFO | 80810 | SLE | GATT | PPT O 
LT | 000'P£ | Z287 | SET T€ SOTO | 8700 | 8800 | ITO | 29€0 | TEVO | 3890 | AIT T | 96c v | 992 1 
OT | 000'7£ | TET | 98/66 | 8:00 | 9900 | CAT O | ITTO | F£0'0 | Pero | 0970 | FLO | 8997 | ZOET 
ST | 000'08 | 379% | 85€ LT | ALTO | 8800 | 4900 | £900 1400 | GLTO | 8690 | TZLO | OLTT | ECVT 
PI | 000'8G | GETZ | I98'SC | 2900 | T90'0 | 80C'0 | 1900 | 9600 | TREO | ZEE'O | 2960 | FEST | CGT T 
ET | 00093 | L9E°% | FE9'EZ | 6000 | 2900 | 88010 | $760 | ITO | 90F'0 | 929'0 | 0071 | TELE IOC 
ZL | 000 v6 | I90' v | 6£6'6T | 7900 | 7200 4800 | €TO'O | TSVO | C90 | C080 | 967€ | GET T | LETZ 
IL|000c|TYCLT | 92c0c | 120'0 | 9800 | 6370 | C00 | 6100 | EFO | TEEO | SPOTT | STE'E | T09T 
OT | 000'0G | F86'0 | 9IO'GT | I800 | $960 | 8910 | 19010 | ESTO | 06€0 | £660 | 6970 | 6880 | ZIT 
6 | 000'8T 008 | 00291 | 96T'0 | 9600 | 22010 6600 | 2870 | 9I€0 | ST8'O | OTZ'G | €68 T | 9380 
8 | 000'9T | CIV T | 88€ FT | €1£'0 | SIT 0 | 9010 | SET'O | 80€0 | PITO | 9820 | 0020 | £06'£ | LOGT 
L [00071 | 0790 | I9S€ $T 9190 | 660 | 8660 | IZZO | ETTO | 0£€0 | 8970 | 2870 | TPLO | C28 T 
9 | 000'cI 0190 | 06€ IT | 8930 | PETO | 2360 | OCTO | 2870 | 8200 | 1670 | 8950 | SET | 4991 
Q | 00001,8690 | ZAF6 | ATO | 06T0 | 9680 | 8600 | PSE'O | 99Z0 | 9080 | F89°0 | 3290 | £560 
P | 0008 | 0680 | IIFA | PETO | C60 | 6870 | STO | TOST | 89€0 | 0280 | 6990 | 980T | 9907 
£ | 0009 | 6£6'0 | 190°S | PEZO | 9870 | 90F'G | £261 | 2910 | 8270 | 62€0 | OI9T | BILT | 2290 
c | 0007 0980 | OP9'€ | LISO | 6920 | £480 | ECHO | TESP | 2270 | 8180 | $890 | 6960 | 886 C 
I | 0007 |6ITO| Z88'T | 2097 | FOOT | 9180 | 2460 | 66870 | 6970 | T9£'0 | 9770 | 6990 | S9Z'I 
L| Pua | wes | MAH | Dy YX X X. X Ð IÐ tÐ ð 19 
“AP OT- PIOUYS9IYY J YALA suorgjoes G 10] Ynsor popuedxo yO "CV 9[qer, 
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Six sections 
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OG | 00070» | E7T'9 | 29888 | 0900 | 6800 | 9PT'O | 8800 FCO | 8700 | 8010 | 8860 | F99 € | 6ETO | Y9T C | GOTT 
61 | 000'8£ | LE0'F | £96'££ | 3700 | 6600 | £800 | 9770 | P6F T | 6290 | 7090 | 6E9'0 | 996 1 | FP € | TOSI | F69'1 
ST | 000'98 | 476% | £90'££ | 08070 | 9700 | 6PT'O | 180'0 | 6200 | £400 | zeen | £TSO | £240 | I8T'G | £960 | TGS T 
LT | 000'P£ | T6/'€ | 60Z'0E | PPO'O | 6600 | 8900 9100 | 9870 | PIO'O | OPO | 222'0 | 1791 | CAF T | €66T | 9091 
OT | 000€ | SFO | 99626 | 60T'O | £900 6800 | LOFO | SAGT | C990 | TIGO | 8880 | OTE | FOE'E | 16870 | SAT 
ST | 00008 | 846 T | ZZ0'8z | 0900 | P6T'O | 8010 | £900 | 6080 | PITT | 90€0 | 2970 | ZLO'T | Z6VZ | 49€ | 86€ 0 
PI | 000'8G | GOTE | 26277 | 6710 | 8900 | 98010 | I800 | S08 T | ESZT | SPTO | 4990 | EET C | 860€ | SOOT | OFF C 
ET | 00093 | 8871 | CIS TC | SET'O | 8900 | TZC'0 | 0/00 | 99T0 | £900 | 99810 | FOP O | OFS'O | 1691 | YOC T | 0860 
ZL | 00076 |089'€ |0c£'0c | 9900 | PET'O | 9800 | STC'O | 6860 | 9870 | FFF'0 | 86810 | 9676 | 9€9T | 9qZ' | 9040 
IL|000'cc | €90°% | LE6°6T | £660 | 6200 | 9900 0400 | OTS'G | £6€0 | €66€0 | ZZS'0 | 6860 | GOTT | 2790 | £010 
OT | 92661 GETT 929'8T | 800 | 3670 | 810 | 6600 | FO T | POTT | ZOSO | SPFO 9860 | LENT | 9781 | 610 1 
6 | 000'8T | POTT | 968'9T | FIC'O | ten | 6600 Leen | 6600 | SIKT | F€C'O | 2290 | 96690 | POEZ | OIAT | 1607 
8 | 000'9T | 8800 49681 | 6TC'O | 9010 | Z80°0 4660 | LETT | SLOT | TIGO | 8660 | ZE6Z | Z86'% | 809 T | 9060 
A | 000'PT | EIT 992 21 | VETO | 850 | 2660 | 69080 | PITO | TOSO 68F0 | SPSO | 040 | 228€ | SOFT | PISZ 
9 | 000'cT | £060 | 260 IT | SCTO | 24€ '0 | 997'0 9060 | Y89'T | Z6CT | 167:0 | 9970 | EFRO | GET | 6067 | COST 
G | 000'0T | E40 | 84€6 | SETO | 16£'0 | ESZO | 8870 | 1670 | IS00 | PETO | 2080 | STGO | 0280 | €680 | TOT 
y | 0008 | PLO 9977 | 0000 | A80 | 8870 | LEGT | S9T'0 | OCTO | TOPO | 88€0 | GIOT | HES Z | PST | SPST 
€ | 0009 4490 | EEG | 0900 | €79'0 | 0890 | 99€% | S6E'0 | 90690 | 62€0 | 9T90 | 9180 | SETT | IC40 | ZS0'T 
c | 0007 | 9290 | eege | PESO | £0 | 6001 | 49/0 | PLZT | 8€0 | FEO | GIGO | PESO | TICT | 9961 | 8980 
I | 0007 | 1400 | 6:61 | 040 | 0607 | PET'S | STAT | 8£7T | GC | 97€0 | 8810 | Z970 | ETFO | £990 | 866 T 
L| Pua | wes | MAH | °X SX YX d X N Ð iÐ ES? Ð Ð Ð 








"OD OT- PIOYS9IUI J YIM suor)oəs Y 10} 3[nso1 popuedxo VI ON 
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A.9.2 Threshold level: -20 dB 


One section 


Table A.7: GA expanded result for 1 section with I" threshold -20 dB. 





























G4 X4 | RBW | Start End | T 
1.028 | 0.564 | 0.450 | 1.549 | 1.999 | 1 
1.014 | 0.282 | 0.901 | 3.100 | 4.000 | 2 
1.026 | 0.188 | 1.349 | 4.644 | 5.993 | 3 
1.009 | 0.141 | 1.798 | 6.196 | 7.994 | 4 
1.017 | 0.113 | 2.250 | 7.742 | 9.992 | 5 
1.021 | 0.094 | 2.694 | 9.268 | 11.961 | 6 
1.022 | 0.081 | 3.135 | 10.784 | 13.919 | 7 
1.032 | 0.071 | 3.579 | 12.337 | 15.915 | 8 





1.026 | 0.063 | 4.029 | 13.864 | 17.893 | 9 
1.037 | 0.056 | 4.488 | 15.496 | 19.984 | 10 
1.022 | 0.051 | 4.950 | 17.029 | 21.978 | 11 
1.000 | 0.047 | 5.379 | 18.606 | 23.985 | 12 
1.022 | 0.043 | 5.852 | 20.132 | 25.984 | 13 
1.007 | 0.040 | 6.295 | 21.703 | 27.998 | 14 
1.018 | 0.038 | 6.734 | 23.167 | 29.901 | 15 
1.017 | 0.035 | 7.201 | 24.775 | 31.976 | 16 
1.005 | 0.084 | 7.555 | 26.071 | 33.626 | 17 
1.040 | 0.031 | 8.074 | 27.920 | 35.993 | 18 
1.016 | 0.030 | 8.439 | 29.036 | 37.475 | 19 
1.022 | 0.028 | 8.968 | 30.851 | 39.819 | 20 
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Two sections 


Table A.8: GA expanded result for 2 sections with I' threshold -20 dB. 

























































































Gi Go Xi X» RBW Start End T 
1.274 | 0.149 | 0.604 | 1.768 | 0.499 | 1.439 | 1.937 | 1 
2.120 | 0.566 | 0.324 | 0.318 | 1.409 | 2.455 | 3.863 | 2 
1.286 | 0.404 | 0.251 | 0.261 | 3.286 | 2.246 | 5.531 | 3 
1.699 | 0.480 | 0.178 | 0.183 | 3.485 | 3.691 | 7.176 | 4 
1.543 | 0.499 | 0.172 | 0.175 | 4.215 | 3.571 | 7.787 | 5 
2.017 | 0.718 | 0.340 | 0.111 | 4.547 | 6.597 | 11.144 | 6 
1.630 | 0.405 | 0.114 | 0.103 | 6.024 | 6.581 | 12.605 | 7 
1.364 | 0.389 | 0.105 | 0.114 | 7.101 | 5.419 | 12.521 | 8 
1.588 | 0.465 | 0.093 | 0.095 | 7.509 | 6.763 | 14.272 | 9 
1.635 | 0.495 | 0.075 | 0.077 | 8.629 | 8.426 | 17.054 | 10 
1.720 | 0.498 | 0.063 | 0.062 | 10.107 | 11.013 | 21.120 | 11 
1.573 | 0.438 | 0.058 | 0.063 | 10.876 | 10.304 | 21.180 | 12 
1.421 | 0.433 | 0.058 | 0.060 | 13.435 | 10.038 | 23.473 | 13 
1.337 | 0.381 | 0.049 | 0.056 | 14.480 | 11.375 | 25.855 | 14 
1.368 | 0.439 | 0.052 | 0.053 | 15.923 | 11.153 | 27.076 | 15 
1.542 | 0.489 | 0.040 | 0.042 | 17.033 | 14.830 | 31.863 | 16 
1.427 | 0.372 | 0.043 | 0.040 | 18.071 | 15.166 | 33.237 | 17 
1.612 | 0.457 | 0.036 | 0.037 | 18.454 | 17.498 | 35.952 | 18 
1.456 | 0.521 | 0.041 | 0.040 | 19.902 | 15.111 | 35.012 | 19 
1.424 | 0.489 | 0.041 | 0.039 | 20.802 | 15.250 | 36.052 | 20 
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Three sections 


Table A.9: GA expanded result for 3 sections with I' threshold -20 dB. 































































































Gi Go G3 Xi X» Xa RBW Start End T 
1.377 | 0.533 | 0.321 | 0.784 | 0.791 | 0.813 | 1.506 | 0.495 | 2.000 | 1 
1.365 | 0.508 | 0.313 | 0.402 | 0.405 | 0.400 | 3.015 | 0.976 | 3.991 | 2 
1.410 | 0.576 | 0.317 | 0.253 | 0.265 | 0.262 | 4.467 | 1.533 | 6.000 | 3 
1.230 | 0.541 | 0.273 | 0.201 | 0.198 | 0.192 | 5.995 | 2.005 | 8.000 | 4 
1.376 | 0.573 | 0.273 | 0.155 | 0.151 | 0.164 | 7.359 | 2.641 | 10.000 | 5 
1.385 | 0.480 | 0.300 | 0.131 | 0.129 | 0.125 | 8.907 | 3.093 | 12.000 | 6 
1.522 | 0.489 | 0.261 | 0.107 | 0.114 | 0.109 | 10.199 | 3.802 | 14.000 | 7 
1.462 | 0.507 | 0.309 | 0.102 | 0.100 | 0.095 | 11.935 | 4.065 | 16.000 | 8 
1.181 | 0.518 | 0.252 | 0.088 | 0.087 | 0.087 | 13.444 | 4.556 | 18.000 | 9 
1.263 | 0.518 | 0.284 | 0.077 | 0.076 | 0.076 | 14.836 | 5.164 | 20.000 | 10 
1.280 | 0.496 | 0.273 | 0.065 | 0.066 | 0.064 | 15.906 | 6.095 | 22.000 | 11 
1.270 | 0.499 | 0.274 | 0.068 | 0.067 | 0.064 | 18.030 | 5.970 | 24.000 | 12 
1.281 | 0.523 | 0.273 | 0.061 | 0.062 | 0.063 | 19.485 | 6.476 | 25.961 | 13 
1.428 | 0.561 | 0.334 | 0.057 | 0.057 | 0.057 | 20.942 | 7.059 | 28.000 | 14 
1.313 | 0.513 | 0.261 | 0.051 | 0.052 | 0.057 | 22.272 | 7.690 | 29.962 | 15 
1.445 | 0.497 | 0.266 | 0.049 | 0.047 | 0.053 | 23.475 | 8.525 | 32.000 | 16 
1.357 | 0.487 | 0.270 | 0.045 | 0.047 | 0.043 | 25.058 | 8.942 | 34.000 | 17 
1.410 | 0.548 | 0.341 | 0.042 | 0.042 | 0.043 | 26.718 | 9.283 | 36.000 | 18 
3.654 | 0.964 | 0.380 | 0.014 | 0.034 | 0.038 | 23.548 | 14.452 | 38.000 | 19 
0.541 | 1.229 | 0.417 | 0.007 | 0.028 | 0.036 | 23.780 | 16.006 | 39.786 | 20 
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Four sections 
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OG | 0000? | 8FS°ZZ | CST LL | TEO'O | 600 | G00'0 | T60°0 | 6990 | LIET | 1480 | 0880 
GI | 000'88 | ¿z€ °T | 6E9'E | 28010 | £00 | 9000 | O£TO | C4€'0 | FOUT | F€ | ISZ 2 
ST | 000'98 | 0£€ FI | 049 IC | 6£0'0 | 07070 | EPPO | ICTO | POFO | GSTT | CITO | IRT 
LT | 0007€ | 969 TI | POE'6T | 9900 | 9£0'0 24100 | ITO | S8€'0 | 6960 | Z8T8€ | OGE'£ 
OT | 0003€ | ISTST | GFS'ST | 7700 | C900 | ZOTO 9000 | 68€'0 | 98TT | 66€ 7. | SLET 
ST | 000'08 | 9172 | GRE ZT | 600 | 1700 | 6700 9700 | OFTO | 9080 | 02770 | POST 
PI | 98622 | 90401 | 2C AI | TS0'0 | 1900 | 0700 | OLT'O 0480 | 3560 | TELE | OLST 
ET | £96'G6 | FEL OI | 6IS'ST | 88010 | 800 | STO'O | TREO | 96£'0 | 160 T | LICG | 662€ 
ZL | 96886 | 916 | 626'ST | 190'0 | 9900 | ZE0'0 | 0100 | 90£'0 | 2090 | 98/T | G88T 
IT | 00077 | 0089 | 089'ST | €900 | 120'0 | 0900 | 4000 | 9870 | ¿190 | 90F T | EEFT 
OT | 00003 | 7292 | 97€'GT | 120'0 | 8900 | a70'0 | BTO | 0010 | 0660 | 06247 | TITT 
6 | F86'2T | EETL | 69801 | 2200 | 8900 | 60070 | 8000 | 08570 | TOC T | 2560 692 
8 | 000191 | ZIT? | F88 IT | 6800 | I800 | 800 | 8800 | TIEO | 9670 | 06€0 | F8OT 
A | 000'pT | FORS | 9698 | TOTO | ZOTO | 1800 | €T€'O | 5680 | $660 | $08 € | 238 T 
9 100071 | 8ZOF | 7262 | 90TO | ITTO | OCTO | 8000 | POZO | YC9'O | TEOT | 6680 
Q | 000'0T | S968 | SFO'9 | 9PTO rt | 090 | 2880 | 2880 | 9901 | TEGE | 099T 
P | 0008 | Z8? | SI€€ | TOTO | 9270 | 3270 ZI00 | 8890 | SEKT | ST6'T | 198°0 
€ | 0009 | P6YT | 409v | eco | SECO | OTTO | SOT'O | 20€0 | 1270 | SZOT | OSE G 
c | 0007 | 0780 | ORTE Iren LE7'O | 89810 | 0370 | 0070 | ZZE'0 | ECCO | ZOE T 
T | 000% | 0470 | OET | 940 | 99/0 | EELO | 00F0 | 2400 | 2870 | 6680 | L68°% 
L| puq | weg | MGU | ÓX od X "X ro ð eÐ IÐ 
“AP 0c- PIOUYS9IYY J YHA SUOT}OeS p 10J Ynsor papuedxa VI :0IV ALL 
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Five sections 


112 

















































































































OG | ST€'6€ | 19Z°9% | LGO'ST | 08010 | ZETO | TOTO | ZS0°0 | ZRZ'O | 1690 | GEIT | 9Z9'T | 1660 | OLTT 
61 | 967'££ | IZOTZ | SAV'GT | 28010 | OTTO | 6070 | £660 | 60810 | £790 | IPFT | TEOS | TETT | 9690 
ST | $9066 | IS6'8T | £80'0T | Z70'0 | LZ1'0 | OTTO | F8C0| 4090 | C4//0 | 25€ 6 | 0621 | 6T | OFTT 
LT | 6787€ | ECT TC | 229 TT | 9800 | 9600 | £460 | 04/0 | LISO | £260 | PE9T | 8070 | 6T0'€ | 90877 
OT | 12616 | S€ 81 | 95€ 6 | 7700 | ITO | 8600 | £900 | £870 | 8890 | VERT | SLOT | OSO'T | LOE T 
GT | 86120 | 67921 | 69€6 | SFO'O | 6710 | 8790 | 88T0 | 4640 | 8890 | 089 | LEVT | ZETT | LEGT 
PI | OU Z | 660'CI | 900'€T | TSO'0 | 2800 | 00:0 | 9270 | £0€0 | c8€0 | 00€ T | L6T'% | EET | €67T 
ET | O€7PG | 29991 | PPSS | 02010 | ECTO | 0900 | 99T0 | 96810 | £170 | 6777 | 664T | cG0€ | GIST 
ZL | SSTOT | 8072 | 82/8 | 9800 | I800 | C940 | 80 | IG T | PIDO | GITT | 92€ € | SETT | 697% 
TT 123607 | I6V'ST orl | 6900 | OLT'O | 6870 | 8800 | 8880 | €29'0 | 860€ | GIST | 8900 | SLS 
OT | OTZ'6T | 00731 | 018 | 6900 | 06T'0 | IST'O | 8620 | 89910 | 0/0 | TELT | 99 T | PEF'T | 2810 
6 | 6q€ ¿1 | 190 1 | c6c'8 | 8900 | £070 | 960 | ETFO | 8601 | 6920 | 986 | STOG | 68TT | Un 
8 |8P8'£T | 90 | org | TOTO | OTTO | OTTO | SOCI | 0680 | 2880 | 6:80 0961 | 2907 | 2891 
A | $166 | 90r Y | 10€ 8 | SET'O | 98TO | 8770 | 16010 | 69T'0 | 0370 | Z90'T | 891 | LOUT | TT C 
9 | 6Z8'0T | 6969 | 016 € | TITO | 67£'0 | IT | TPSO | ZETO | ZE9'0 | 889 T | ISKT | 698 T | 68T 0 
Q | 7996 | 9609 | ZIE | TO | C£TO | 9870 | ICTO | 6800 8650 | LEIT | STAT | 2841 | 2060 
P | 2404 | 084€ | 26€ € | ISTO | 99T0 | 6690 | 99TO | Z6Z'T | 89810 | I9CT | 6c€ 6 | OLT'T Leen 
€ | POTS | 094% | FORT | TOZO | $970 | OFTE T | SOFT | 9601 | 8010 | €6T'T | OSOT | BOLT | Leen 
c | LOTS | 4S/T | TOGT | 96810 | I9€0 | I9TO | CO | SST | FTO | SE€T | FTO'G | 8880 | 0690 
I | 6941 | 2090 | GOTT | 0160 | 6640 0690 | CEO | FIOT | 0970 | 8690 | GICT | C28 T | STET 
L| puq | wes | MAH | °X YX EX ex på Ð Ð tÐ e) t9 
“HP 0c- PIOYSOIYY J YIM suonoos G 10} Haat popuedxo WO “ITV AQLL 
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Six sections 
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OG | 9:998 | Z6£'ST | Sec I6 | ¿€0'0 | 60:0 | 90T'O 0080 | 8770 | $00 | ISO | SITT | TELE | 6000 | SLET | ALRT 
GI | CZAS | TET 9T | OPS 'G6T | 200 | 0700 | 0800 | 260 | 95/0 | LEZO | EGO | SETT | LEET | 666 T | 1971 | 929 1 
ST | 68198 | 6€ TT | 46071 | S£0'0 | 98010 | 6900 | OTT | 9/60 | ITTO | 2880 | ZZF'T | STO'G | 8681 | 8000 | SOFT 
LT | 92T'8£ | ¿090% | 04€ 71 | 98010 | ££0'0 | SPTO | 8800 | £970 | 280 | TECO | SPOT | FEET | 600€ | 9STT | ZE9T 
OT | SZYTE | 9ITOG | STS'TT | 28010 | OCTO | €90 0 | TOGO | OIZT | C660 | 8690 | 9€S'T 6892 | SPAT | 90€ | TEST 
GT | 98087 | GEL AT | 6€ 0T | £000 | ZET'O | 2680 | CIT | £860 | OCT TE | 8490 | 296 | I8 | TOF'G | TOT | ZETT 
PI | 92610 | 888 II | 680701 | 6900 | 68010 | 9270 | GEL T | 992'0 | PPO'O | £870 | 06€ | OFG'T | EET | €180 | PTT 
ET | 70887 | GLE LT | 8818 | A900 | £2€0 | 1900 | 6070 | 0280 | TOT | PEDO | GIT C | LIST | 666 | AT | ECVT 
ZL | £9P'£G | OPO'ST | €782 | 0900 | 68TO | £280 ¿ESO | SPOT | T870 | 1040 | 687% | 009% | 6273 | OCT | 9TOT 
IL|A4c90c|06c8T1 | 2874 | 8800 | 6/10 | 6000 | Z£9'T | 2992 | TE9O | £690 | 0900 | 2966 | SATZ | CVLO | 7460 
OT | 90861 | ESTE | ZZT'L | 89010 | 2810 | ¿PEO | LOZO | STO | LLT | 9690 | BEAT | LOST | CEST | 4871 | 08010 
6 | V6C PT | 60c8 | 9809 | 800 | 6.0 0 | 9:000 | TOTO | S&z'0 | OF9 T | F9€0 | SPRL | 96€ T | 096T | 8887 | one" 
8 | 88€'ST | Z0TOT | ISC'8 | 8200 | 2330 | 618'0 | OGTG | 8610 | 8600 | 998910 | LIET | 9980 | 96V | 66v c | 9890 
A | SETET | 1298 | STOP | 1600 | 6600 | 8890 | C880 | PCT | GAT'T | £690 | 688 T | OET | £46 T | EET | CFLS 
9 | 260 IT | 1849 | STEP | Z0T0 | 2880 | TREO | 8490 | 96TT | 9671 | TIGO | OET | LOUT | ICGT | TOGT | SSE T 
Q | 66896 | I£V9 | SOF'E | ICTO | 89810 | I9c0 | 89/0 | IGTO | 0290 | 0690 | TILT | G6FT | 86€T | 7020 | 80910 
P | $664 | 7667 | 676% | VETO | 280 | 0980 | SETT | 8010 | 98T0 | 990 | I6/T | 969 T | GOET | 2690 | COCO 
£ | P888 | ZEGE | LOGT | 0070 | CF9'0 | PEZO GEI | 8810 | £790 | 7990 | LPL T | OOF'T | 7980 | Dron | 6980 
T | £998 | L07 | 68ST | 8880 | 8280 | 9970 | OSET | 800 | 80810 | FZE'O | SETT | SOLT | TTGT | OPE T | ShZ'O 
T | 49/1 | 9490 | I60T | OG80 | 669'0 | 660 | 6£0'% | LEO T | TYKT | 6880 | 09910 | SIT T | 6880 | 7060 | ZFZ Š 
L| Pua | es | MAH | OV Ý YX EX eg D4 Ð Ð ES Eb) Ð t9 
‘AP 0c- Prousaxq3 J YIM suorjoes 9 10} 3jnso1 papuedxa YO ZTV e[qer 
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A.10 Transmission Line Design and Calculation 


In this work, the first step in transmission line design is to determine the width- 
to-height (%) ratio for parallel plate line on an FR-4 board dielectric. Here we use 
subscript "ms" to denote "microstrip" and "pp" to denote parallel plate. We use 


the following formula for microstrip [40, 41]: 


Ww 8e^ 
COM = 22 (A.22) 
where: 














Er T 


(7). = = |B-1-mn@B-1)+ $= (n(B-1)+0.39 - 52) 


(A.24) 


3771 
B= ——— A.25 
2Z0/€r ( ) 


Where e, = €g = 4.8 for FR-4 board and Zo = 25 Q (for microstrip). The height 


of the board is fixed at i inch for parallel plate (5 inch for microstrip). Plugging 
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in the numbers, we get: 


(=) = 5.1086 (not used) 
msA 


h 
w 
— — 4.9614 
(>) a 
1. 
hms = 35 inch = 0.0794cm > w = 0.3938 cm 


Then we stack two microstrips together in series like from Fig. A.5 into A.6 to 


form 50 €) parallel plate line. The 7 ratio for parallel plate case is: om = 2.4807. 


Copper tape 


VAN Yes 


FR-4 board W hf 


LREBBLELEBLEBLELLELELEBLEBLELELELELEBLELEBLELEBLEESE 





Ground plane 
FIGURE A.5: Microstrip line 


Copper tape 
4 







E-field 


WD 





FIGURE A.6: Parallel plate line 


Finally, we can figure out the effective permittivity of the board using Eq. A.26 


below: 





(A.26) 
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Epp = 3.9276 eo (A.27) 


In conclusion, we know the width, height, and effective permittivity of the trans- 
mission line. The length of the board can be easily calculated after we decide on 
operating frequency (which gives us the wavelength and separation between resistor 
elements). 


In section 5.1, the dimensions are: 


Wg = 10.16 cm (Width of board) 
Lg = 15.24 cm (Length of board) 
wo = 0.3937 cm (Width of copper tape) 
tg = 0.15875 cm (Thickness of FR-4 Board) 
eu = 4.8 (Relative permittivity of FR-4 Board /Epoxy-glass) 


(A.28) 
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A.11 Details of Location Measurements 


A.11.1 CIEMAS 


Table A.13: Details of measurement locations in CIEMAS. 
































Location Category 
Outside 1411 near router 
Outside DiVE virtual reality lab corner 
Outside 3467 hall, near router 
Outside 3478A room 
3rd floor S wing near stairs NE corner near router, corner 
3rd floor S wing near stairs S corner near router, corner, near window 
3rd floor skyway hall, near window 
2nd floor conference room above Twinnies room 











A.11.2 Teer 


Table A.14: Details of measurement locations in Teer. 
































Location Category 
2nd floor near information board corner 
near 213D corner, near window 
outside basement lab near router 
106 room 
114 room 
outside 118 near window 
outside 423A near router 
breakroom 4rd floor room, near router 
outside P04 near router 
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A.11.3 Physics 


Table A.15: Details of measurement locations in Physics. 





Location 


Category 





Outside 151 


hall 





Outside 161 


near router, corner 





Outside 179 


near router, hall 





101 commons 


room 





Outside 029 


near router, corner 





Outside 014 


near router, hall 





Outside 06 


hall, near router, near window 





Outside 035 


near router 





NW basement exit 


near router, corner 





Outside 094 


near router, corner 





Outside 285 


hall, near router 








Outside 0010 





near router, hall 
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A.11.4 French 


Table A.16: Details of measurement locations in French. 





Location 


Category 





Outside 2225 


hall 





Outside 2336 


hall, near window 





Outside 2313 


near router, corner 





Conf room near 2313 


room, near router 





Outside 2303 


near router, hall 





Outside 2117 


hall 





Outside 2135 


corner, near router, near window 





Outside 1124 


corner 





Outside 1125 


corner 





Biosci hall in front of 1207 


hall, near router 





Outside 0069 


near router, hall 


























1243 room 
Basement corner, near router 
Outside 0012 hall 
Outside 3135B near window, corner 
Outside 3232 hall 
Outside 4223 corner 
5316 room 
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